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PREFACE 


Part X, Sections 1 through 10, -and Part XI, Section 1, contain a large 
body of background and explanatory material. Part of this material is 
necessary to an understanding of the analysis of the MFASMR system, which is 
presented in Part IT., Section 2, in response to Task 1 of the work statement. 

The result of the analysis in Fart II, Section 3, indicates that the 
MFASMR system does not exhibit high performance because much of the available 
thermal power is not u#cd in the construction of the image and because the image 
which can be formed has a resolution of only ten lines. Because of this 
unfavorable result, Task 2, which called for a Comparison with filled aper- 
ture systems in terms of weight, complexity, reliability, deployment, storage, 
performance, cost and development risk, became unnecessary. Therefore, the 
study effort was redirected to emphasize the analysis of image reconstruction 
and an understanding of aperture synthesis. 

Task 3, image reconstruction, or, more formally, the image quality com- 
puter study, was performed in three ways: 1) two-dimensional inverse Fast 

Fourier Transformation, 2) Wiener filtering plus the inverse Fast Fourier 
Transformation and 3) maximum entropy reconstruction. The’ results of these 
computer analyses are summarized in Part II, Section 2, and are presented in 
detail in Volume II, Appendix. In addition, the computer programs are 
delivered as part of the contract effort. 

Task 4, comparison with conventional aperture synthesis systems, has been 
completed. This effort is reported in Part I, Section 8, and in Volume II, 
Appendices. The appendices in Volume II also contain discussions in more 
detail of the theory of radiometry by correlation and aperture synthesis. 

Because referencing is implicit in the nature of this material, a bibliog- 
raphy for Part I is in Section 11, Part I, and a bibliography for Part II is 
in Section 6, Part II. In the appendices in Volume II, the references are 
footnoted within each appendix; 
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1. INTRODUCTION 


1.1 RECAP OF THE STUDY OF THE MFASMR SYSTEM 

This report describes a study of a device for high resolution radiometric 
mapping. This Multifrequency Aperture-Synthesizing Microwave Radiometer Sys- 
tem (MFASMR) is an example of a class of imaging systems which form images by 
a two-step process. They offer an alternative imaging method which competes 
with more familiar imagers such as lenses and pencil-beam scanning antennas. 

Direct imagers sense the individual resolvable elements in the imaged 
scene (object space) and build up the image by recording the intensity of each 
resolvable element (pixel). Lenses display all the pixel intensities simul- 
t'.iiu'ouM |,y on ;i nitrfacv called the image plane. In photography ul.L Intensities 
are recorded simultaneously on light-sensitive film. In the microitave region 
the most successful direct imaging method has been the use rf scanning, 
pencil-beam antenna. The antenna beam looks at the scene pi.'_.t.s, one-nt-n- 
time. The antenna output power is measured by a microwave receiver. Then a 
display device, such as a CRT, forms image pixels with brightness proportional 
to the receiver output. They are ordered' to match the correspondingpixels 
in the original scene. 

Such image-making processes are familiar due to long, usage. The MFASMR 
and related devices form images in a manner which is not so familiar. In 
addition to the image plane, there is a second surface Ln optical systems 
called the frequency plane or Fourier plane. A set of measurements in this 
plane can be used as the input to a Fourier transform operation. The trans- 
form output is the image. MFASMR and similar systems make such measurements 
ln the Fourier piano in order to form images. 

In a little while it will be shown that the electromagnetic flux in the 
frequency plane is the two-dimensional spectrum of the image-plane field. 

This spectrum is a spatial frequency spectrum. Each spectrums are not as 
well known as more familiar temporal frequency spectrums. However both 
spectrum types possess many analogous properties. 
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There are several methods of making the frequency plane measurements. 

The most common method involves movement of an antenna over the Fourier plane 
while cross-correlating its output with the output of a stationary reference 
antenna. This process is used by many radio telescopes. It is called aper- 
ture synthesis. Alternately the antenna pair can remain stationary and the 
operating frequency varied in order to obtain the required samples. MFASMR 
uses this approach. 

In t ho case of aperture synthesis, the collection of samples is usually, 
divided into two steps. The first step consists of movement of the roving 
antenna along a track. The fixed and roving antennas are at the ends of a 
variable baseline. The result can be thought of as an interferometer with 
variable antenna spacing. After the roving antenna reaches a nmvimnm dis- 
tance from the fixed antenna, the baseline is rotated to a new orientation. 
The roving antenna repeats its excursion along the new baseline. Eventually, 
the baseline rotates 180°. If then coincides with its original orientation, 
but with baseline ends interchanged. Since the antennas are identical, the 
Interchange of baseline ends does not affect measurement. Therefore, the 
measurements during the first 180° rotation can be used for the measurement 
values over the next 180° rotation. Thus, samples over a circular disk in 
the Fourier plane are obtained. In common parlance, only the measurement 
process which occurs as the roving antenna moves along the baseline Is called 
aperture synthesis. The second step, consisting of rotation of the baseline, 
is called rotation synthesis or suporsynthesis. 

f 

Now MFASMR does not use rotation synthesis. Instead, it uses' antennas 
at the ends of two fixed baselines at right angles to each other. Also, 
MFASMR cannot operate at a set of frequencies dense enough to meet the 
Nyquist criterion. Only frequencies where there are no transmitters can be 
used. Thermal radiation is very weak compared to man-made sources, .. Such 
sources easily destroy the radiometer accuracy. Only bands reserved for 
radiometers, or bands where transmitters are not operating for other reasons, 
cun be used. Thu MFASMR design identifies five such spectrum windows at 2.4, 
4.75, 9.5, 19.0 and 38.0 GHz. Both the undersampling and- poor distribution 
of samples damage the data collected. 

Some of the damage done to the data can be repaired, if we are willing 
to use more general methods of forming the image from the spectrum samples 
than merely using the Fourier transform. These processes of forming images 
from measurement data urc called image reconstruction. 

Compensation for undersampling the spectrum depends on the fact that the 
image signal is redundant. Because of this redundancy, missing data can be 
rcpluccd with vulues predicted from knowledge of sample values around the 
missing value. The average difference between this predicted value and the 
true value can be made much less than the average difference to be expected 
when redundancy is lacking. 
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If che restoration process is linear, there are analytic methods called 
rate distortion theory for computing the rms image distortion as a function 
of the amount of data lost in the process of sampling and coding the image 
data* This rate distortion theory lias been used in this study to compute the 
percentage of image distortion to be expected in MFASMR systems with plausible 
parameters. It is found that about 30% image distortion is to be expected. 
This distortion is too great to render this approach attractive. 

However, this result does not settle the question as to whether MFASMR 
is a viable system or not. A significant restriction was placed on the image 
reconstruction processes when rate distortion theory was used to determine 
system performance. Namely it was postulated that image reconstruction was a 
linear process. In fact, reconstruction need not be linear. 

During the last decade, a group of nonlinear data processing techniques 
have been developing which are called by names such as "maximum entropy 
method" or "maximum likelihood method." 

These methods are recursive. Multiple cycles around the recursion loop 
maximize the output data entropy or other closely allied parameters. Signal 
til.'iL I «t I oh do not need to bo known u priori, aw in linear reconstruction 
methods. Data statistics are taken into account by their influence on the 
recursion process. 

An oxoo I (out ruv low ol’ llionu muLliodn In provided by uu IEEE i'russ Reprint 
Volume, Modern Spectrum Analysis, edited by Donald Childers Spectrum estima- 
tion, image reconstruction and adaptive arrays are treated as aspects of the 
maximum entropy method (MFM) and its relatives. 

MEM processing has been controversial, in spite of the excellent, results 
often obtained. Although it often outperforms linear methods, the MEM recur- 
sion cycle has turned out to he unstable in many cases. The causes of insta- 
bility iiuve not been well understood. In addition, the maximization of 
entropy, or its twin, information, has appeared to be a vary abstract objec- 
tive compared to such down-to-earth objectives as minimization of rms data 
distortion. 

However, Stephen Warncko (References 1, 2, 3, A) has developed a MEM 
algorithm for image reconstruction from frequency plane samples which has con- 
straints built into the recursion loop which ensure stability. It lias yielded 
excellent results compared to linear reconstruction methods. Dr. Wernecke has 
kindly given us a FORTRAN program using his MEM algorithm. The results 
obtained using this MEM algorithm to reconstruct images from simulated MFASMR 
data are given in Volume II of this study report. 
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The simulated MFASMR data is produced by use of microwave radiometer 
image signals provided by Robert Moore, Head of the Microwave Radiometer Group 
at MOTS, Inyokern. These images were generated by a Ka-Band pencil-beam, scan- 
ning radiometer. A Fourier-transform operation turns these images into spatial 
frequency spectrums. These spectrums can now be sampled in the same manner as 
MFASMR samples the actual spectrum. Those sample sets, an well as sample sets 
corresponding to other frequency-plane imagers, provide the data output for the 
computer study of MEM image reconstruction reported in Volume II. 


1.2 STUDY EXTENSION TO THE PROPERTIES OF FOURIER-PLANE IMAGERS 

There are several themes which will be developed in this study which 
show why frequency-plane imagers are superior to direct imagers for orbital 
mapping. At this point, these major themes will be sketched in broad out- 
line, permitting the reader to see the goals of the detailed system analysis 
as it develops. 

The first and most important theme is a direct extension of the MFASMR 
approach, namely, the conviction that frequency-plane imagers will be the 
dominant type of radiometric orbital mapper in the future. This conviction 
stems from this study of frequency-plane imagers. 

Later in this report, the nature of the frequency plane will be sharpened 
up. However, Llie most importune property of the frequency plane needs no 
further analysis. The frequency plane is the electromagnetic field, generated 
by the scene being viewed, at the surface of imaging lens or reflector, which 
acts as a decoder of this f rcqucncy-p Juno data . Decod- 1 «■„ transforms the dutu 
into a form we can understand. This form is an image the scene. 

Building better hardware for the purpose of decoding the frequency-plane 
data impinging on optical systems has preoccupied people since the time of 
Newton. 

In the last few decades a new idea has grown up. This new strategy for 
forming images involves the two-step process mentioned in Section 1.1. 

Instead of transforming the frequency plane data into an image by means of 
optical hardware, the frequency plane data is measured by means of small 
probes moving over this plane. Once this data has been recorded, the function 
of the imaging lens or reflector is carried out by a computer program. 

Neither the probes which measure the frequency-plane nor the computer 
are anywhere near as large as the hardware they supplant. This is especially 
true when the hardware supplanted is a very large scanning antenna constructed 
to do high resolution radiometric mapping from orbit. The replacement of the 
imaging antenna by a computer algorithm has many other advantages besides 
small sias. The computer program is highly flexible. The hardware is not. 
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We can construct programs for transformation from frequency-plane data 
to image with many different properties. We can transform in such a way as 
to make the image the most faithful replica of the scene. Or we can trans- 
form so as to obtain the best visual appearance. 

If wc please, we cun measure the properties of the image obtained and 
then change the transform computation so as to make the Image more and more 
like siome standard defining a "best" image. This is comparable to adaptive 
arrays. However, we do not need to make physical changes in a hardware 

Imager, such ns the sotting of phasu shifters. Our adaptive process is only 
a recursive loop in a computer program. 

Indeed, wo can transform the same data in many different ways in parallel 
computations in order to output images of the same scene optimized in many 
different ways simultaneously. 

This approach focuses attention on the transform processes which have 
replaced the imaging antenna. As a result, our theme of replacing hardware 
imagers by software imagers causes heavy emphasis to be placed on these 
transform processes. These processes are coming to be known as "image 
reconstruction. " 

The second theme is the observation that although frequency-plane 
imngcrs use small antennns compared to direct imagers of comparable resolu- 
tion, the Lr sensitivity and ability to resolve temperature variations in the 
mapped region is just as good as the direct imagers. 

Roughly speaking, the antennas of Fourier imagers arc sized to sec the 
untire field of view subtended by the region being mapped. This selection 
of antenna beamwidth is dictated by the consideration that these antennas 
measure periodicities in the brightness of the scene which extend throughout 
the entire mapped region. Thu mugniLudu and piiase of uii possible periodici- 
ties comprise the spatial frequency spectrum of the scene. It is this spec- 
trum which is Fourier-transformed to form the final image. Ths scanning , 
anteuna of the competing direct imager will be sized to resolve a single pixel 
In the scene. If the final Image lias M by M resolvable lines, tbc scanner 
will have an area greater than the area of the Fourier imager antennas. 

Now, it is an interesting fact that £he large scanner antenna and the 
little frequency-plane antennas have the same total capture cross section. 

That is, their capture cross sections integrated over spheres surrounding the 
antennas are the same. A fundamental antenna theorem states that if the 
antenna is matched to its transmission line feeder, this total cross section 
is always when the operating wavelength is X (Reference 5) . It is easy to 
see that, as a consequence of the above theorem, the sversge power received 
by the scanner antenna end the frequency sampling antenna is the same. 


T.-1-5 
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In order Co prove this assertion , let the thermal flux density from the 
scene be p (watts/nr-steradian) ; let the solid angle subtended by the f ield- 
of-vicw be ft. If there are M? pixels in Che scene, the solid angle subtended 
by the scanning beam of the direct imager. Aft, is equal to fl/M^. The fre- 
quency plane antenna is sized to reduce its cross section outside of its field 
of-viev (FOV) to zero. Therefore, its average cross section must satisfy the 
relation 


^FPA^ 


or 




(I-l-l) 


The thermal power received by the frequency plane antenna is 


p^.pA ■ (flux dunslty) (l f OV solid angle) (average cross section) 




watts 


iu the case of the seamier, ideally the antenna cross suction 1 outside 
the beam is reduced to zero. Therefore 


^eam 




( 1 - 1 - 2 ) 


so the power received by the beam is 


P BEAM 



( 1 - 1 - 3 ) 


2 

Therefore, the power received by the two antennas is the same, namely pX . 

If we concentrate, our attention on a single pixel, we see that it 
delivers a power of pX^/M^ to the frequency-plane antenna. The same pixel 
delivers power to the scanner in pulses. When the scanner beam falls on the 
pixel, it delivers power to the scanner equal to PbeAM or pV-. However, the 
antenna must scan all of the other M* pixels, so it spends only (l/M^) of the 
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total time looking at a single pixel. So the average power delivered by the 
pixel to the antenna is the duty factor times the power level during the 
pulse. This is (l/M 2 )(pX 2 ) - pX 2 /M 2 . So the average power delivered by a 
pixel to the sensor is the same in both cases. 

Since aperture synthesis uses two antennas, the fixed reference and the 
roving antenna, which see the entire FOV, the frequency-plane imager actually 
receives twice the power received by the equivalent direct imager. However, 
the number of frequency plunc samples required for input to the Fourier trans- 
form operation is twice the number of pixels produced. Usually, the sample 
doubling occurs because M 2 frequency plane samples are taken which are complex 
numbers. The noise power of the real and imaginary parts of the sample add 
to double the noise in the M 2 pixel-bin outputs from the transform. Since 
both noise power and signal power are doubled, the signal-to-noise ratio (SNR) 
of the frequency-plane imager is the same as the SNR of the direct imager. 

Another source of confusion lurks in the radio-telescope literature. 
Hewish (Reference 6) asserts that an aperture-synthesis device has a sensi- 
tivity which is reduced by the factor, (number of positions of the movable 
element)"*^/ 2 , as compared with a conventional antenna of the same resolving 
power. For proof he refers to an earlier paper by himself and Ryle 
(Reference 7). 

Without attempting to find the earlier obscure reference, it is easy to 
sec that the factor arises because hewish is considering a point source. In 
that case, the Fourier plane samples add incoherently (power-wise) ’ while the 
scanner antenna adds the point source signal inputs to the array coherently 
(amplitude-wise). This is true if the scene is only one point source, such 
us a radio star, however, if r.he scene is a surface which radiates non- 
coherently, the thermal radiation from all the surface elements will add 
power-wise in both the aperture-synthesis device and the direct imaging . 
antenna. This is true even when the scene is a set of point sources which 
are not coherent with respect to each other. For this case, Hewish's factor 
disappears. This observation is proven by the detailed derivations of the 
performance of various imager types given as appendices to this report. Thus 
Hewish's observation is true, but applies only to the trivial case of a scene 
consisting of a single point source. 

Scanning the literature will show that the correct treatment of coherency 
conditions for frequency-plane imagers caused much trouble before these condi- 
tions were rigorously defined. Hence, it is not surprising to find observa- 
tions like the one by Hewish in the early papers on Fourier imagers. This 
theme will be worked out by detailed sensitivity calculations. 

A third theme, which will be of importance in this paper, is the observa- 
tion that rotation synthesis is a technique which will be used independently 
from aperture synthesis. Usually rotation synthesis is thought of as only 
part of the total aperture-synthesis process, rotation of the interferometer 
baseline being needed in order to make measurements over an area. However, 
rotation synthesis is a technique which stands alone and is vary useful even 
when aperture synthesis is not used. 


1-1-7 
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As Kraus demonstrates (Reference 8) , the multiplying interferometer 
measures the visibility function, V 0 (Sx) e^ffS as the baseline length in 

wavelengths, S\, is varied from zero to some maximum value D. This function 
will be identified as the spatial frequency function of the scene at a later 
time. 


Let this function be measured over a set of baseline orientations, 6, 
extending over 180° in angular displacement. Let the set of orientations be 
dense enough to make the measurements of V 0 meet the Nyqulst condition. Then 
it can be shown that the two-dimensional Fourier transform of V Q is an image 
of the scene in the region selected by the interferometer antennas. 

In order to have a compact notation, we will indicate a Fourier- 
transform pair of functions like this: 


G(f) v===ip(x) OR THIS C(f) - F{P(x)} 


(1-1-4) 


where G(f) is a spectrum of some kind and P(x) is a function of a linear or 
angular displacement. Then the assertion in the last paragraph that the 
image is the Fourier transform of the two-dimensional visibility function 
can be written in the form 


2— d 

V(S X , 0) * * B(», 0) 


(1-1-5) 


where 4 is the nadir angle and B(<t>, 0) is the radiometric brightness ''of scene. 

A valuable feature of this process is the fact that all the baseline 
orientations can be measured in time succession yet no sensitivity is lost, 
if the interferometer antennas are correctly sized to receive power from a 
f ield-of-view defining the scene which is to be mapped. This is in strong 
contrast to systems where B(4>, 0) is mapped by a pencil-beam scanner. 

The scanning array or reflector covers an .aperture area as shown in 
Figure 1-1. The comparable frequency-plane probes are shown to the right of 
the scanner aperture. 

The total array or reflector must be present all the time. If the aper- 
ture is "thinned," flux from the scene streams t.hrough the gaps created by 
ChLiut Lug. if the apurturu is thiimud enough to produce a significant weight 
saving, sensitivity is also significantly reduced. The thinned array is of 
course just as big as the original antenna. Therefore, thinning is of no help 
in stowage and deployment. 
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FIGURE M. PENCIL-BEAM SCANNER CONTRASTED WITH ROTATION 
SYNTHESIS IMAGER HAVING THE SAME SENSITIVITY 


Up to now we have described the measurement of V(S. , 6) by aperture 
synthesis. 

However, Kraus (Reference 8) in his Eq. (6-83) gives the relation 


V(S X ) ^u(+) 


( 1 - 1 - 6 ) 


where g(4>) Is the projection of the source brightness on the axis defined by 
the baseline oriented at angle 6. That is B($) is the' scene brightness as a 
function of the semiapex angle of the conical antenna patterns of the linear 
.'irr.’iy. These patterns define hyperbolas on the surface defining the sccnu. 
For small FOV, Eq. (1-1-6) is obvious. 

Now B(«p) can be measured by any line-source antenna. ' Then V(S^) can be 
computed in a data processor and stored. This can be done for a set of line 
source orientations identical to the baseline orientations used in the 
aperture-synthesis case. This is the same V(S^, 8) generated by aperture 
synthesis. The two-dimensional transform indicated in Eq. (1-1-5) again pro- 
- duces the image. 

Thus we have the option of using a line source, instead of aperture 
synthesis, to measure V(S^, 8) and then form an image from" this function. 

This can be a valuable option. Aperture synthesis involves the move- 
ment of one interferometer antenna over an area in the 1 Fourier plane. This 
involves quite a bit of mechanical complexity. 


1-1-9 
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In addition, uperture-synthesis is limited in bandwidth by the requirement 
that the outputs of the two interferometer antennas must correlate. Some line 
sources do not have this restriction. This can be of great importance when 
high resolution and good sensitivity are desired at the same time, 

A final theme is emphasis on coherency. Frequency plane measurements 
are coherency measurements. Therefore, a correlating radiometer is needed to 
make these measurements. A correlating radiometer is not a radiometer mode 
but one of two possible radiometer types. The other radiometer type measures 
power, rather limn coherency. Either radiometer can operute in one of many 
possible modes, such as Dicke, total-power, phase-switching, etc. So the 
division of radiometers into coherent end noncoherent (power-measuring) types 
Is comparable to the classification of radars as either coherent or noncoherent. 

As in the radar case, the theory of coherent radiometers is more subtle. 
However, the capabilities of correlating radiometers are much greater than 
the capabilities of radiometers which only measure power. It should not be 
thought that because the theory is moze complicated, the equipment is also 
mora complicated. It is often mentioned in the literature that the correlator 
is more complicated because it requires two receiver channels. However, in 
most remote sensing applications data about both signal polarizations are 
required. Two noncoherent radiometers (two receiver channels) are needed to 
measure both polarizations. Now both polarizations can pass through a single 
correlating radiometer by injecting them into the correlator in phase quadra- 
ture. So again, two polarization measurements require two channels. 

** l 
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2. BASIC CORRELATING RADIOMETER CONFIGURATIONS 


2.1 COHERENCE MEASUREMENT 

Before looking at frequency-plane imagers in detail, it is necessary to 
review the characteristics of correlating radiometers. They differ from 
ordinary radiometers in that they measure coherence or correlation between 
two channels. The ordinary radiometer measures power in a single channel. 

In many ways, these radiometer types mirror the classification 'of radars into 
noncoherent (power-measuring) and coherent radars. 

Correlating radiometers operate in many of the same modes as noncoherent 
^ radiometers. For instance in both cases the simplest mode of operation is as 
a total power radiometer. The correlating version Is shown In FIriitc* 2-1. 

« 

(. 


TOTAL POWSft 



OUT 


FIGURE 2-1 . TOTAL POWER CORRELATING 
RADIOMETER 

t 

f < • 
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The input to this radiometer configuration is shown as two antennas. 

This is the way the antenna is used is aperture synthesis. One antenna is 
fixed at the end of the interferometer baseline. The other antenna is the 
roving antenna which moves along the baseline. Of course* when the maximum 
baseline length is long, the two channels are not fed directly from the 
antennas with transmission line. Instead, the signals are reduced to IF, 
amplified, and distributed. 

In the case of Fourier- transform spectrometers, very often there is only 
one antenna. This single signal is then divided in a hybrid. The hybrid 
outputs go to the mixers. However, a preferred configuration is still the 
two-antenna configuration. In this dual configuration there is ho portion of 
thu signal channels which is common to both channels. Therefore, none of the 
radiated thermal power associated with the channel losses appears in the out- 
put as part of the signal. Since the radiation comes from different lossy 
elements, it does not correlate and is rejected. 

This makes radiometer call! ration much simpler. The effect of internal 
radiation and of reflected radiation due to channel mismatches is very com- 
plicated. A paper by Otoshi (Reference 9) bears out this point. 

A similar problem, which is even harder to control, is due to mixer 
noise. Due to the passage of bias current through the mixer diodes, they act 
as noise sources. The effective source temperature is ib out 400°K. This 
radiation travels out through the microwave channels. £ portion is reflected 
at mismatches such as hybrids, circulator, bends, and the antenna itself. 
Although the radiation is noise, it is coherent with itself. Therefore, a 
standing wnvo Is formed. Tf a mismatch produces, total reflection, the peuk 
of the mixer-noise standing wave will be (1 + 1)^ x 400°K ■ 1600 # K. The mea- 
sured mixer-noise error is proportional to a voltage at the radiometer output. 
It can be seen that if 0.16°K is allotted for mixer-noise error, the reflected 
mixer noise must be reduced by 10^ (-80 db) . This assumes that the standing 
wave loop can be anywhere in the microwave channel. If its position can be 
stabilized in the channel, the errors produced can be calculated out. How- 
ever, this can be hard to do, especially in wideband systems. 

Figure 2-2 shows the effect of mixer noise on a National Radio Astronomy 
Observatory radio telescope. The long-period fluctuation is due to mismatch 
at the antenna input. The short-period fluctuation with about 10*K amplitude 
is due to reflection off the edge of the Cassegrain subref lector . 

The dual antenna setup makes it easier to obtain the required 80 db iso- 
lation. Mixer noise reflected back to the mixer making the noise finds 
nothing to be correlated with in the other channel and is eliminated by what- 
ever means the radiometer uses to eliminate other kinds of receiver noise. 

Only mixer noise which couples from one channel to the other will be treated 
like a signal and appear in the output. Careful design of the dual antennas 
can ensure that the coupling between them is very small. When isolators are 
also used in tha two channels, the required 80 db isolation is obtainable. 
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.FIGURE 2-2. EFFECTS OF MIXER NOISE ON MEASURED ANTENNA 
TEMPERATURE 


In Che case of a pencil-beam scanning antenna, no increase in thermal 
power results from using two antennas. Both antennas must produce beams 
which tiue the kciiic pixel, at the Name time. i*'ur instance, in Che cuhc of a 
reflector, the two antennas would be produced by using a dual faedhorn like 
that used for a monopulse feed. However, in the case of aperture synthesis, 
each antenna may see the total without robbing power from the other 
antenna. Therefore, a two-to-one gain in signal power is obtained for this 
case. 


Total power radiometers do not modulate the input signal. As a result, 
fluctuation of the receiver gain produces fluctuation in receiver noise which 
cannot be distinguished from the signal. For instance, let the system noise 
temperature be 2000*K, let the gain stability be one percent (AG/G ■ 10 -2 ), 
then the output receiver iioIno fluctuations will be 20*K rus. This is un 
order of magnitude greater than typical temperature resolutions. 

Engineers here at Hughes Space and Communications C?uup have achieved 
remarkable gain stabilities. As a consequence, they have been using the total 
power mode. However, the conventional attitude has been that it is too hard 
to get the required stability under field conditions. Therefore, the usual 
design approach has been to modulate the thermal noise before it becomes 
mixed with receiver noise. The modulation then permits separation of the 
signal from the receiver noise in the video portion of the radiometer. 
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The original radiometer of this type was the Dicke. It had an RF twitch 
which twitched the receiver input between the antenna and a load. The tem- 
perature difference between load and antenna then appears as a thermal signal 
which is square wave modulated. After square-law detection, the square wave 
becomes a voltage wave. Synchronous demodulation by the RF switch actuation 
waveform then picks the thermal difference signal out of the receiver noise. 

The corresponding correlation radiometer is the phase-switched radiometer 
shown in Figure 2-3 on the following page. 

* 

Both of the radiometer modes just discussed are not fully coherent. The 
output of the correlator Is at zero frequency. Let the input to Channel 1 be 
fA(t) + jB(t)]. The input to Channel 2 is [A(t + t) + JB(t + r)J. The cor- 
relator .multiplies these waveforms together and averages the products to give 
the correlation fun ction in the form [A(t) Alt + t) - B(t) B(t + t)!|i + 

2j[A(t) B(t + f) + A(t + E(t)]2. It d* the modulus of this final expression 

which emerges from the radiometers we have been discussing. We cannot get at 
the rcsl component of the correlation, [ nor the imaginary component, 

[ In the same manner, the in-phase and quadrature receiver noise powers 

add power-wise in the output. 






PHASE SWITCHED 



FIGURE 2-3. PHASE-SWITCHED 
CORRELATING RADIOMETER 
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This same series of operations on the signal occur in a homodyne radar, 
except that one channel is replaced by a sinusoidal reference. Both systems 
are degenerate in the quantum mechanical sense. Two separate signal modes, 
corresponding to two eigenvalues have assumed the same values and have become 
indistinguishable. 

In order to avoid this defect, the two radiometer channels can be offset 
in frequency by using two local oscillators. These LOs are frequency offset 
by means of a servo loop* or by feeding one LO signal through a phase-shifter. 

Tl» lt» shifter keeps Increasing the phase change at a steudy rate in order to* 
produce the frequency offset. 

The output of the correlator consists of DC terms sucli as ( t ) , which 
arc removed and the cross product terms given above. This signal is demodu- 
lated by two phase detectors. One phase detector is a switch whose switching 
waveform or reference is in-phase with the frequency offset. The other phase 
detector reference is in quadrature with the offset. The in-phase detector 
selects only [ Jj,. The quadrature detector sees only [ ]2* 

A block diagram for this frequency-offset mode is shown in Figure 2-4 
below. 

Most modern radio telescopes use this frequency offset technique (Refer- 
ences 10, 11). Older instruments use the phase-switching method. They recover 
thu in-pliuHu and quadrature components by shifting the signal from one end of 
the baseline 90° with respect to the signal from the other end of the baseline. 
Then they perform a second homodyne detection. These two correlator outputs 
are often called "promp" and "semi-promp." 

Wessellng (Reference 12) is one of the first to point out the substantial 
advantages of the fully-coherent radiometer over the earlier circuits. 

Figure 2-5 shows the microwave portion of a correlating radiometer. The 
signal enters through the hybrid at the center of the assembly. The two chan- 
nels are at the top right and bottom left. At the extreme right a noise diode 
mount can bo soon. Tin* diode* output is Lnjuctud into butli channels in u 
wervoed bridge arrangement for measuring the unknown signal temperature. 


2.2 TEMPERATURE MEASURING RADIOMETERS 

The earliest use of radiometers was the measurement of radiant tempera- 
ture. These early instruments were noncoherent Dicke radiometers. At first 
tho voltage out of the phase detector was measured directly. This voltage is 
at least ideally proportionate to the difference between the reference load 
and the antenna temperature. 
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FIGURE 2-4. FULLY COHERENT 
CORRELATING RADIOMETER 
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FIGURE 2-5. 1.5 cm CORRELATING RADIOMETER 
RF HEAD (PHOTO 79-62146) 


Ic was soon found chat amplifiers were noc stable enough to maintain an 
accurate ratio between antenna temperature and output voltage. The obvious 
answer is a servo. One approach is to servo the gain (Reference 13). The 
other possibility is to servo the reference load temperature to the antenna 
temperature by adding thermal noise power to the output from the cold load. 

The noise is generated by a hot standard and Injected through a directional 
coupler (Reference 14). 

TlK'tt'rvmd ul tenu.’ilor iii'.Ikim* urvur went anywhere. After uil, ;m error 
oi 0.1 db cor responds, to an error of 3*K or so. Radiometers are better at 
calibrating attenuators than vice versa. Neither has the servoed gain 
approach been too successful for reasons which are obscure. 

The most successful solution has been the Hach servo loop (Reference 15). 
In this loop, a calibrated source, square-wave modulated, is Injected into 
either the antenna channel or cold load channel as the case may be. Then the 
duty factor of the square wave is servoed to balance antenna and load chan- 
nels. Initial calibration is carried out by noting the noise source duty 
factor using standard temperature loads. This is the system used by the 
radiometer in figure 2-5. 

The unknown temperature is then found by calculations using the duty 
factor values. The duty factor measurements are time measurements. Time 
measurements can be made with very great precision. 


1-2-7 
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Typical results for a Hach loop radiometer are reported by Hardy et al 
(Reference 16) and by Love and Van Melle (Reference 17) . A block diagram of 
the Love/Van Melle radiometer is given in Figure 2-6. 

Figure 2-7 shows a measurement of the boiling point of liquid argon. It 
agrees to 0.1°K with the tabulated boiling point of 87.3*K given for the 
existing air pressure. As a matter of fact, this radiometer maintained an 
accuracy of 0.1°K throughout an entire flight test program without recallbra- 
tion of any kind. 

However, the stability of transistor amplifiers has become good enough to 
make temperature measurements with an unservoed radiometer, if accuracy specs 
nrc not too demanding. The radiometer is calibrated while its feed horn looks 
at loud near 290°K, and again while the horn looks at a cold load or the sky 
(Reference 18). As the title indicates, the radiometer whose performance was 
measured was a total-power radiometer. However, this calibration technique 
is applicable to any radiometer including correlators. 

The final type of temperature measurement of importance is the thermal- 
gradient mode. This is a mapping mode where the antenna scans a pencil beam 
fast enough to cause all the itarmonics of the mapped image to occur at fre- 
quencies higher than those where gain fluctuation and 1/f noise is important. 

The amplified square-law detector output is passed through a bandpass 
filter. This filter's upper cutoff Is set by the usual considerations when 
designing a smoothing filter. The lower cutoff is set to reject the 1/f 
noi6e. 

Tills simple radiometer can make excellent microwave radiometer images. 
Removal of the low frequencies improves the image appearance. However, loss 
of these frequencies prevents making absolute temperature measurements. 

Figure 2-8 shows a thermal gradient radiometer designed and tested by 
Wiley Electronics in 1960. The antenna is a spherical reflector illuminated 
by three rotating line source feeds which correct for spherical aberration. 
This Is Lite first mapping microwuvu radiometer. It made excellent images. 


2.3 SPECTROMETERS 

When spectrum measurement is mentioned, usually the first thought is the 
use of filters (Reference 19) . If speed is not too important, a single scan- 
ning filter will do. In order to get more speed a bank of channelizing micro- 
wave filters can be used. • This is the route taken by the writer in designing 
a "frequency scanning" mapping radiometer for JPL. This radiometer's ampli- 
fier traveling wave tubes (TWTs) and filter bank are shown in Figure 2-9. 
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FIGURE 2-6. DETAILED BLOCK DIAGRAM OF A PRACTICAL DESIGN FOR A RADIOMETER 
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SYSTEM RESPONSE TIME T • 1.0 nc 

MLEtHENCt TEMPERATURE T c • 303.3k 

ANTENNA AVIPAOI TEMPERATURE T p - 1B.0°C 

VERTICAL SCALE Ik PER INCH 

FIGURE 2-7. LIQUID ARGON TEST RECORD. 3-8-72 



FIGURE 2-8. MAPPING MICROWAVE 
RADIOMETER MOUNTED IN C-82 
AIRCRAFT 
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FIGURE 2-9. FIFTY -CHANNEL FILTER 
BANK DRIVE BY TWO LOW NOISE TWTS 


The channelizing filter banks can be seen at the top of the cabinet. 

Sunk* of Hie «J I Mod van Lnil«.m til' lit lit approach arc obvious from the photograph. 

The filters are bulky. They are also lossy; so considerable gain between the 
antenna and filter bank is needed. In this case the gain was provided by 
two TWTs. 

Another flaw appeared too late in the program to rectify. The sensi- 
tivity r the individual channels drifted with temperature changes, time, 
drive level, and other clauses which could not bu identified. 

Although the basic mapping method was proven out, this radiometer made 
poor images. The uneven channel reponse causes heavy streaking in the longl- 
t ml Inn I (flight path) direction. 

This radiometer was built over a decade ago. Modern components would 
make it possible to build a much better instrument. Feedback loops in each 
filter channel could eliminate uneven response. However, there is a more 
subtle deficiency which the competing Fourler-transform radiometers do not 
have. 

This deficiency is the fact that the performance parameters of the spec- 
trum estimation procass are set by the properties of the filter bank and 
associated equipment. This means chat these properties are "cast in concrete." 
There is no way of changing them short of rebuilding the radiometer. 
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The alternate spectrometer design is Che Fourler-transf orm spectrometer. 
This device measures the autocorrelation function of the signal. The Fourier- 
transform of the autocorrelation is then found. The Wiener-Khintchine theorem 
states that this transform is the power-density spectrum of the signal 

(Reference 20). 

To form the autocorrelation function the signal must be split into two 
components. A differential time delay between channels must be Introduced. 

The delayed signals are then multiplied together and averaged to form ♦ (t). 

Figure 2-10 shows a correlator which uses a variable ultrasonic delay line. 

Figure 2-11 shows a correlation function measured by this radiometer. The 

measurement is made at various noise levels in order to measure sensitivity. 

If a variable path length (trombone) is used to produce the differential 
delay, we have the equivalent of the scanning microwave filter. Such a device 
was designed and tested by Long and Butterworth in 1963 (Reference 21). They 
did not transform <)(?) to get die spectrum. Instead they drove the trombone 
at a steady rate. The trombone output is Doppler-shifted by the trombone 
movement. 

A little thought will show the Doppler spectrum appearing at baseband in 
a replica of the microwave spectrum. Long and Butterworth therefore measured 
the RF spectrum by measuring the video spectrum with a low frequency spectrum 
analyzer. Figure 2-12 shows a block diagram of such a variable delay spec- 
trometer. The spectrum analysis is carried out digitally in accordance with 
contemporary design practice. 



FIGURE 2-10. ELECTRONIC CORRELATOR WITH VARIABLE 
QUARTZ DELAY LINE TO PRODUCE DIFFERENTIAL DELAY 
BETWEEN TWO SIGNALS (PHOTO 700-64-3 A) 


1 - 2-12 





zimon 








Cl~tt€0|« 
























ORIGINAL PAGE IS 
Of POOR QUALITY 


Usually Che emphasis on performance dictates simultaneous measurement of 
all values of $(r). Figure 2-13 shows such a system. The values of are 

measured at the Intermediate frequency because the correlator Is very wideband 
(GHz). The two IF signals are fed to opposite ends of two transmission lines. 
At the center of the lines the delay in both channels accumulated from the 
point where the signals split in the hybrid, is the same. For taps on one 
side > of the center point t is increasingly positive with increase in distance 
from' center. Going the other way from center t is more and more negative. 

The directional couplers tup power from both lines. The circuitry driven 
by each tap pair is similar to that shown in Figure 2-4 starting at points A 
and B. The processed tap-pair outputs are then multiplexed, digitized and 
fed to a computer whore the Fourier transform of <}>(t) is found. 

If the phase-shifter in Figure 2-13 is correctly indexed, the outputs of 
the square-law detectors can be exactly in phase with the I reference. The 
Q-chnnncl. electronics Is then not needed. 

The digital values of $(t) are fed to a computer where a Fourier trans- 
formation is carried out to yield the spectrum estimate. 

If the bandwidth of the system is narrow enough, the signal can be con- 
verted to base-band, digitized and delayed in shift registers. The register 
outputs are then multiplied by the original undelayed digital signal in digi- 
tal multipliers and Integrated In digital, accumulators to form the set of 
autocorrelation function values. Welnreb contrasts the analog and'digital 
‘approaches in his thesis (Reference 22). His Figure 4, shown here as Fig- 
ure 2-14, gives functional block diagrams for the two correlator types. 

Weinreb ' s correlator used one-bit quantization in order to make the 
digital circuitry simpler. It should be remembered that, in 1963, LSI had 
not arrived yet. A decade later Bowers ct al. report on muitibit correlator 
(Reference 23). 

A tutorial paper on autocorrelation spectrometers is given by Cooper 
(Reference 24). This paper is completely oriented toward digital correlators. 
In Section 3.5.10, Cooper reviews the radio-telescope correlators operating 
at that time. This data is summarized in Table 2-1. 


2.4 IMT£KF£KOH£T£RS 

The simplest type of interferometer is the adding interferometer. The 
output of the antennas at the end of the interferometer baseline ure added in 
a hybrid and passed to the receiver as shown in Figure 2-11. 

The interferometer pattern consists of s ssriss of lobss. Thsy arc all 
positive as indies tad on the pattern diagram in Figure 2-15. 





















FIGURE 204. ANALOG AND DIGITAL AUTOCORRELATION 
METHODS OF SPECTRAL MEASUREMENT 
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TABLE 2-1. OPERATIONAL CORRELATION SPECTROMETERS (1976) 



Clock Rate 





Location 

(MHz) 

BW 

Channels 

Bits 

Comments 

Nit A<), Croon Bank 

0.039 to 20 


. .. 

613 


. 

Mil, Lincoln Laboratory 

50 




Interleaved 

Sampling 

U.S. National Radio 
Observatory 

50 





Jodrell Bank 
Observatory 

* 

10 

256 



C.S. I.li.O. Radio 
Physics Lab 

20 


1024 



Uiiiv. oC British 
Columbia 


16 

128 

1 


Univ. of British 


P 

128 

3x5 

Dual Mode 

Co liunb I ;i 

25 


1024 


i 

Haystack Observatory 

50 

■ 

512 


Trimode 


100 


256 





FIGURE 2-15. SIMPLE ADDING INTERFEROMETER 
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The D.C. component of thiB pattern gives no information about the scene. 
Its presence complicates the use of the data collected. A much better type 
oL' interferometer is the multiplying interferometer. This type can be mecha- 
nized by use of a correlator to carry out the multiplication. The phase- 
switch^ interferometer is also a multiplying interferometer. These possi- 
bilities are shown in Figure 2-16. 

Note that the lobes are alternately positive and negative. There is no 
DC component. Figure 2-17, taken from Krause (Reference 8), shows the output 
of tho adding and of the multiplying interferometer as a source traverses 
their fields of view. 



•I CORRELATION FUNCTION b) PHASE-SWITCHED 


FIGURE 2-16. MULTIPLYING INTERFEROMETER 



•I SIMPLE (UNSWITCHED) INTERFEROMETER OUTPUT 



b) PHASE-SWITCHED INTERFEROMETER OUTPUT 


FIGURE 2-1 7. RESPONSES OF AN ADDING 
INTERFEROMETER 

i 
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A slgnlf leant rot! U fiction on thu use of un interferometer in caused by 
decorrelation of the signals from the two antennas due to their differential 
time of arrival at these antennas. Roughly, this differential time must be 
less than the bandwidth of the signal. 

There are several ways of casing this restriction One popular method 
is used when the instantaneous required FOV is small, but moves through 
large angles over long periods of time. For instance, in the case of radio 
telescopes the FOV is often determined by fairly large parabolic antennas. 

The parabolas sec only a small part of the sky. However, over periods of 
hours, this small FOV moves to another part of the sky. The solution in this 
case is to switch in pieces of transmission line in the path back to the cor- 
relator. by tming tin* right line length the differential timu oC travel to 
the correlator can be made close to zero for e source at the center of the 
FOV. 


in many casus the spectrum of -the signal emerging from the interferometer 
antennas is desired. If the spectrum is found by passage of the signal through 
a bank of channelizing filters or by transform spectroscopy and cross correla- 
tion is carried out after this process, decorrelation is usually not n problem, 
liuch spectrum chaiuiui cun bu made nurrow enough to ensure correlation in the 
presence of the differential time delay between antennas signals existing as 
a result of other design considerations. A sketch of this scheme follows. 

Obviously the urrangumuut In Figuru 2-18 is needed for atmospheric 

sounding. Sounding measurements involve measurement of the spectrums of atmo- 
spheric absorption lines. A radiometric version of frequency scanning will 
bo discussed Inter in tills report. It also uses a multichannel system con- 
figuration. The types of correlation functions which can be generated is 
discussed in E. Tensor Synthesis. 



FIGURE 2-18. MULTICHANNi i. MULTIPLYING INTERFEROMETER 
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For the sake of completeness Intensity interferometry should be 
mentioned. This technique detects the signals in the two channels first. 

Then the envelopes of the two video signals are correlated. It can be shown 
that the correlation coefficient of the squares of the moduli of the antenna 
voltages (detected video) is equal to the square of the correlation of the 
Incident radiation (Reference 25) . As one would suppose, this eliminates the 
affects due to the RF structure of the radiation. Therefore, correlation can 
be obtained over much greater baseline lengths. Also, no phase information 
need be transmitted to the signal receiver modules at the ends of the baseline. 
The detected signals, therefore, can be recorded on tape and correlated at 
soma site remote from the antenna sites. 

This technique is very-long baseline interferometry (VLBI) . 

However, removal of the signal RF structure from the correlation process 
exacts a heavy penalty. The sensitivity is considerably worse than RF cor- 
relators (Reference 26). Even after reuding some review papers, it is not 
clear to the writer that VLBI can produce images (Reference 27), rather than 
measurements of point sources. Therefore, VLBI will not be considered further 
in this study. 


2.5 POLARIMETERS 

I, ti muny cases, the polarization state of received rudiution is desired 
as well as the quantities already discussed. Polarization can be described 
by several methods. One description is the Stokes vector which is favored in 
optics. 

If a^ and a 2 are the amplitudes of the two linearly polarized orthogonal 
components of radiation, then the Stokes vector is 



<•?> + 


I 


<*?> •;> 


M 

s ■ 

^2a 1 a 2 cos y ^ 


C 


<2*1*2 sin 


S 


where y is the phase angle between a^ and S 2 » This assumes that the radiation 
is quasi-monochromatic. Then y can be defined and measured even for noise- 
like signals. 
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Consider three pairs of antennas: let two antennas be linearly polarized 

at 0* and 90°; linearly polarized at 45° and 135°; finally, two antennas right 
and left circularly polarized. 

The powers available from these antennas are denoted by P^.., P^q, , 
Pj 35 , Pr, and P^, respectively. In optics, these same quantities are self 
measured by dichroic filters, which extract various linear polarization com- 
ponents from a beam of light. They are preceded by a retardation plate.. 

Thun 1(0, c) denotes the intensity of light vibrations in the direction mak- 
ing an angle with coordinate axis OX, when the y-component is subjected to a 
retardation e with respect to the x-component. 

Then the components of the Stokes vector are (Reference 28) 


I - 1(0*, 0) + 1(90°, 0) - P 0 + P 9Q 

M - 1(0*, 0) - 1(90*. 0) - P Q - ? 90 

C - 1(45®, 0) - 1(135®, 0) - P 45 - P 135 

S - 1(45°, ir/2) - 1(135®, ir/2) - P R - P L 


( 1 - 2 - 2 ) 


The Stokes parameters have been satisfactory in optics. However, they 
produce complications when used in microwave interferometry. This is due to 
the fact that nl.l conceivable feed arrangements for a two-element interferom- 
eter response to a maximum of two out of four Stokes parameters (Reference 29) . 
(Ko's letter says "minimum" instead of "maximum." This must be a typo. Equa- 
tions (1-2-2) show that four antennas per baseline end are needed.) 

The coherency matrix is better adapted for use with microwave radiometers. 
It needs only two antennas per station for measurement of two orthogonal 
polarizations. If these polarizations are linear, the components of the 
matrix p can be measured by an arrangement such as is shown in Figure 2-19 
below (Reference 30). The power output from dipole A Is p^. The output 
from dipole B is p„y. The real part of p xy la obtained by cross correlating 
the voltages from A and B. The imaginary part of is obtained by delaying 
the signal from A by A/4 and then cross correlating. 


f 
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FIGURE 2-19. CORRELATION POLARIMETER 
WHICH SIMULTANEOUSLY MEASURES THE 
INDEPENDENT PARAMETERS OF THE COHERENT 
MATRIX (THE FEED CONSISTS OF A PAIR OF 
CROSSED DIPOLES) 


Wo hoc Liion Clint Lilt* coherency matrix Iuih the. Conn 


J J 
xx xy 


if ± 

E„E„ E E* 

2 2 x y 



J J 
yx yy 

m 

EE* E E* 

y x y y 

m 

< 4 )_ 


where 


i[«j» (t)-2irvt] 

12 X U) - a x (t) e 


and 


i[*_(t)-2irvt] 
E(t) - « 2 (t) e 2 
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The trace of J is + Jyy » [E x e£] + [E^E*] , The nondiagonal elements 
are conjugates of each other, that is Jxy * Jy*. Therefore, J is Hermitian. 


The interrelation between the Stokes parameters and coherency matrix ele- 
ments are given below. 


xx 


+ J 


yy 


fj^ - 1/2(1 + M) 


M - 


J - J 
xx yy 


C 



J 

yx 


J - 1/2(1 - M) 

yy 

J « 1/2 (c + is) 

xy 


(1-2-4) 


s 




L J yx * 1/2 < c “ is) 


Often we may wish to measure the polarization state of each image pixel. 
If the image is generated by aperture synthesis, the basic interferometer con- 
figuration due to Ko (Reference 30) follows on the next page. 


Both Figure 2-19 and Figure 2-20 show homodyne correlating radiometers. 

In a practical system it is better to use a fully-coherent radiometer. One 
advantage will now become apparent. It was mentioned in the spectrometer 
section that the signal whose spectrum is measured could be contained solely 
in the in-phase signal channel by proper adjustment. This means that a second 
signal can be passed through the quadrature channel of the same instrument. 

Tin* iiaLurnl itigual to iu.ji'ct. ItiLo the two channels arc the pair of orthogonal 
polarizations. In most remote sensing applications the radiant brightness of 
both polarizations are required. So dual use of the same radiometer con- 
siderably simplifies the sensor. 


Passage of the dual signals through mixers is the most subtle operation 
required in order to mechanize the dual radiometer. Figure 2-21 shows one 
embodiment of such a mixer. * 


The LO signals in both channels (separated in frequency by the off-set 
frequency) are injected into the mixer input ports by means of a coupled- 
resonator filter as shown. The normal operation of the mixer then delivers 
E^ to one IF port and Eg to the other port. 
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FIGURE 2-20. INTEFEROMETRIC POLAR IMETER WHICH 
SIMULTANEOUSLY MEASURES THE FOUR POLARIZATION 
COHERENCE FUNCTIONS FOR THE FOURIER SYNTHESIS 
OF THE POLARIZATION ANO BRIGHTNESS DISTRIBUTION 
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Figure 2-22 shows the total dual-polarization radiometer. Note that 
is split into two inputs to the correlator channels by means of an in-phase 
hybrid. E„ is split into its channel input signals by means of a quadrature 
hybrid. Therefore, we have Eg and Ey in the left channel the right channel 
has E x and E y . The sum and difference of these signals are formed in both 
channels by the hybrids which drive the mixer diodes. In the left channel 
E x + j E y goes to one diode and Ex - jEy goes to the other diode. The output 
hybrid then forms Ex and jEy at its sum and difference porta after the signals 
are reduced to the intermediate frequency. Likewise E x and Ey are formed at 
separate IF mixer ports in the right channel. 

The correlators are mechanized as shown in Figure 2-23. 

From Figure 2-22 we see that E X E£, E y Ey, (EyE*>H mag , and (E y E$) real are 
formed. These are the elements of the conerence^matrlx. 

The servos indicated in the block diagram are Hach servo loops which 
have been discussed before. A functional block diagram is given in Figure 2-24. 

The square wave noise bursts are injected into correlator arms with cor- 
rect phase to act like E x , Ey and jEy. The VCO seeks a frequency where the 
duty cycles of the noise diode outputs balance the correlator outputs to zero. 
These are Clapp calibrators. They will be discussed in more detail in a later 
section. 


2.6 TENSOR SYNTHESIS 

Ko (References 29, 30) calls the simultaneous measurement of the many 
coherence functions associated with the incoming radiation simultaneously 
"tensor synthesis." We will not treat the subject as formally as he does in 
Ri'ft’rciict! (30). 

It should be noted that the most general space-time coherence function, 
rCxi, X2» x) , involves cross correlating signals at locations xi and X2 with 
dl ffuruuLlal (.Line delay i . Then l'(xi, x^, t) is u temporal (longitudinal) 
coherence- function. Its Fourier transform is the radiation power spectrum. 

We also have a lateral coherence function r(x]_, X2» 0). Its transform is the 
brightness distribution of the radiation source according to the van Cittert- 
Zurnike theorem. 

Returning to the general coherence function we see that a generalization 
of these transforms implies that r(xi, X 2 , x) has associated with it three 
other functions. Namely, its „ "spatial" Fourier transform, T(mi, V2» x) ; its 
"temporal" Fourier transform r(*i, X2, v); and its total Fourier transform 
r(wi» M2* 
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FIGURE 2-22. CORRAO CONFIGURATION FOR COHERENT MATRIX MEASUREMENT 




COUPLED-CAVITY FILTER 
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FIGURE 2-23. CORRELATOR CONFIGURATION 
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The question arises then ns to how we can generate the coherent functions 
involved and their transforms. In order to do this assume that antennas are 
located at two positions. Let us label the various variables involved as' 
follows: . 


i. i 


k, 

m, n 
p» q 
y» * 


Polarization components 
Time 

Position (m - y position change, n ■ z position change) 
Frequency 

Spatial frequency components in y and z directions 
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Let the correlators at positions m and n consist of pairs of transmission 
lines, to provide delay, together with output circuitry as shown in Figure 2-25. 
Let these correlators have polarization inputs E^ and Ej . The cross correla- 
tion outputs are Tij (xk) and (xi) . The T* subscripts indicate which ele- 
ment of the polarization coherency matrix is represented by the I\ The super- 
scripts indicate the locations of the interferometer baseline. The t subscripts 
indicate the set of differential delays available. The temporal transforms of 

the two correlator outputs are f“ (v ) and f?, (vj. 

ij P ij 9 

If movement of the interferometer baseline generates a series of 
y-position changes, Xq, and a series of z-position changes normal to y dis- 
placements, we have a set of correlation measurements at x m , x n . If correlator 
outputs with a delay of x * ure correlated, we get )’(x n , n, t) . Its 

temporal transform is ?ij(x m , x n ) . Its spatial * transform is „ 

Finally we can transform r“x(xk) and ^(t*) to form fij(Vp) and r"ji(v^). 

If a set of displacements x^ and x^ have been generatedand F D and r measured 
for these displacements, n we can set v„ - ■ v. Then the 0 two-dimensional 

spatial transform gives rijCxn, x n , v). Its transform is r^jCpy, u z , v) . 

A block diagram of these manipulations is shown in Figure 2-25. 



FIGURE 2-25. GENERATION OF COHERENCE FUNCTIONS ASSOCIATED WITH (x, . x 2 , r) 
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The correlation formed at the bottom of the diagram is very important. 

The bandwidth of any frequency channel is AB ■ v p +i - v p - - Vq, This 

bandwidth can be made email enough to ensure correlation no matter what how 
long the interferometer baseline is. This is important if ve wish to build 
an aperture-synthesis system capable of atmospheric sounding. The required 
system bandwidth may have to include the entire O 2 absorption spectrum around 
60 GHz. This is a* bandwidth of 10 GHz or more. 

A method of frequency scanning radiometers will be introduced in a later 
section. It cun provide resolution and mapping normal to an aperture- 
synthesis axis. This method requires 10% to 201 bandwidths. As in the 
sounding case, this large bandwidth would destroy correlation in the aperture 

synthes I m process, if spectrum analysis were not curried on first. 

Measurement of the coherency matrix has an important function even when 
the polarization data is not needed for remote sensing purposes. Often the 
verticul und horizontal linear polarization components of the incoming radia- 
tion are desired. Frequently conical scan is used in ordar to ensure that 
the beam makes a constant angle with the terrain for all scan angles. Such a 
scan can be produced by n rotating reflector as shown in Figure 2-26. 

If the reflector assembly rotates but the feed horn does not, the polariza- 
tion components received by the feed horn change with scan angle. For instance 
In Figure 2-26 above, the electric field vector is shown In the plane of the 
paper for the reflector position shown. This is "vertical" polarization. 
However, when the reflector rotates 90* and the horn doas not, the E vactor 
will ba normal to the ray paths shown on the drawings. It then Is horizontally 
polarized. 
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If Che system is relatively narrow band, a rotary joint can be used 
between the horn and the rest of the microwave plumbing. Then it can rotate 
while the rest of the radiometer is stationary, however, in many systems the 
feed horn is a scalar horn, operating in many frequency bands, spread over a 
very great frequency range. There are no rotary joints which can pans such a 
set of operating bands. 


However, if the radiometers monitoring these bands are designed to measure 
the* coherency matrix, conversion of polarization components to correct coordi- 
nate system can be carried out by the data processing computer in the process 
of image restoration. This is done by multiplying the original Stokes vector 
[IMCS] by a rotator matrix T(C> to produce the required polarization compo- 
nents. This relation is . , 


“l“ 


"l 

0 

0 

0 


1 

M’ 

m 

0 

cos 26 

sin 26 

0 


M 

C’ 


0 

-sin 26 

cos 26 

0 


c 

Jl 


_0 

0 

0 

1 


s 


(1-2-5) 


Thun wt’ art* ah 1 1 * La have l he it intp l. Lclty offerod b/ only rotating the 
reflector by computing the required components in the radiometer data 
processing section. 
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3. CORRELATING RADIOMETER TEMPERATURE SENSITIVITT 


3.1 RADIOMETER FIGURE OF MERIT 


The subject of radiometer temperature sensitivity will be touched on 
very lightly. For some reason sensitivity has attracted more than its share 
of the attention in microwave radiometer studies (References 26, 31, 32, 33, 34). 
In reality, statistical bias errors are more numerous, harder to control and 
usually larger than the temperature resolution error. Absolute temperature 
calibration errors arc also liable to be greater than the resolution error. 

The temperature resolution of a radiometer is given in the two forms below: 


41 " ^ K t sys Ji ■ 


where Tgys ■ system noise temperature, B is bandwidth before detection, b is 
the final bandwidth after detection, t is the integration or smoothing time 
after detection, and K is the radiometric flgurc-of-morl t. The system noise 
Leuipuruturc is usually cast in tiie lonu 


T SYS “ T REC + t ant 


(1-3-2) 


TrEC is the effective noise temperature at the end of the feedline from the 
receiver which connects to the antenna. is the antenna background tem- 

perature. This temperature is the scene radiant brightness, weighted by the 
antenna cross section. Integrated over a sphere surrounding the antenna and 
divided by 4ir. The result of the averaging operation is then corrected to 
account for antenna ohmic losses. So 


T ANT " (L ANT " iy T o + “$T / o<0 * +> B(0 » ♦> dn (1-3-3) 

‘'sphere 
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where is Che ohmic loss in the antenna, T 0 is ambient temperature* B(6* <fr) 

is scane radiant brightness and a (6, 4>) is antenna cross-section in meter 
squared per stcradian. 

The receiver temperature is 


REC 


* » T o + Vj 



(1-3-4) 


where Lp ■ feed-line loss, T^p ■ noise temperature of preamplifier, Grf ■ pre- 
amp gain, Tm - mixer-and-IF noise temperature. With sufficient preamp gain 
the third term becomes negligible. 

When there is no preamp, the mixer and IF noise temperature dominate. So 


T • T 

a rec a m 


(F - 1) T 
n o 


(1-3-5) 


where F r is the noise figure of the mixer/IF-amplifier. 

if iitcro is uo mixer, wo luivc u tunod-radio-f requency (TAF) receiver. 
Then and F n in (1-3-5) refer to the main amplifier. In superhet radiom- 
eters without image rejection both RF sidebands are received and folded on 
top of each other in the RF amplifier. Then the double sideband noise fig- 
ure, FQSBt should be used instead of the single sideband noise figure, FgSB* 
It is often stated that Fqsb “ F S e«/2. However, the correct relaction is 
(Reference 35) 



(1-3-6) 


When Fggg » 1, the approximation is usable. However, when modern low-noise 
receivers are used, the approximation penalizes single-sideband receivers. 

The low-frequency bandwidth, b, in (1-3-1) is the noise bandwidth of the 
filter /integrator. A rectangular filter which gives the same noise output 
power as the actual filter, and has a gain equal to the peak gain of the fil- 
ter, G , has the noise bandwidth, B . That is 

nu n 


G(f ) df 

B - 

n g max 



Area Under Filter Gain Curve 
Maximum Power Gain of Filter 


(1-3-7) 
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I Table 3-1 relates the 3 dB band widths to noise bandwidths for output low pass 

filters. 

Often one sees the noise bandwidth given as the bandwidth to use for the 
high-frequency receiver bandwidth. This is not correct. The correct band- 
width to use is called by Evans and McLclsh the radiometer reception bandwidth 
(Reference 31). This bandwidth is defined as 



(1-3-8) 


TABLE 3-1. PROPERTIES OF LOW-PASS FILTERS 


Filter Type 

B 3 dU 

d noise 

B N0I8E 
B 3 dB 

T 

Ideal Low-Pass 

B 

B 


l 

2B 

RC Filters with Time 
Constant of 

* 


. 


Single RC 

i 

2ttRC 

1 

4RC 

1.57 

2t ac 

Dual RC 

0.643 

2irRC 

1 

8RC 

1.22 

4t rc 

Triple RC 

0.510 

1 

1.16 

6 >C ; 

2irRC 

32R6 

Gaussian (® RC) 

j 



1.07 

i 

2.12 B 3 

Box Car Filter 
(Ideal Integrator) 

0.443 

T 

1 

2ir 

1.13 

! 

* 

5-Pole Butterworth 



1.01 


5-Pole Bessel 



1.04 
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Benda t and Piersol (Reference 36) call the same quantity the equivalent 
statistical b andwidth. The most sophisticated treatment is given by 
Bracevull (Reference 37), who uses the autocorrelation bandwidth. Table 3-2 
lists reception bandwidths defined by Evans and McLeish. 

The most detailed treatment of the sensitivities of correlating radiom- 
eters is given by Faris (Reference 38). For the homodyne type of correlator, 
where the correlation output appears at DC, we find K to be /2. This is .given 
in his Eq (Reference 22). This result is reasonable. K for a single-channel, 
total-power radiometer is unity. The correlator has two channels. Therefore, 
the noise power in the output is doubled. The measured temperature is propor- 
tional to the square root of the power after the signal passes through a square- 
law detector or multiplier in a nutocorrclutor . Therefore, the correlator K 
will be /2 . 

It is unfortunate that Faris did not go on to find the value of K for a 
fully coherent correlating radiometer. The reader will recall that this 
radiometer type first heterodynes the incoming signal down to an intermediate 
frequency. The IF signals are then summed in a power combiner and detected 
by a square-law detector. The DC terms out of the detector arc rejected by a 
bandpass filter. Due to u frequency offset introduced between the IF channels, 
the signal appears on a carrier and is passed by the bandpass filter. The fil- 
ter output is then synchronously demodulated in two phase detectors. The 
phase references arc in-phase and in quadrature with the IF frequency offset. 


TABLE 3-2. RECEPTION BANDWIDTHS OF BANDPASS FILTERS 


Filter Type 

B RECEPTI0N /B 3 dB 

^OISE^ dB 

Cascaded SYN tuned CRTS 

3.14 

1.57 

Number of Stages 



1 

3.14 

1.57 

2 

1.96 

1.22 

3 

1.76 

1.15 

4 

1.62 

1.13 

(«• Gaussian) 

1.50 

1.07 

Rectangular 

1.00 

1.00 
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It should be obvious that synchronous detection rejects the system noise 
in quadrature with the reference signal. Furthermore, when the noise origi- 
nates in two receivers it is evident that the noise power in phase and in 
quadrature with the reference will be equal. Therefore, the in-phase and 
quadrature baseband signals resulting from phase detection will have half the 
noise power which emerges from the bandpass filter. The corresponding value 
of K. is unity. 

However, those who want a more rigorous proof of these facts will find 
them given by Lindsey (Reference 39). 

This reduction of the noise in each of the two output channels has long 
boon known by radio astronomers (References 40, 41) . it is a mystery to the 
writer why this case continues to be ignored in the general literature. 

The value K ■ 1, just given for the fully-coherent radiometer case, 
assumes that there is one input signal that is divided in a hybrid in order 
to drive the two signal channels. However, often we can use two antennas. 

For instance, in the case of aperture synthesis there are two antennas at the 
ends of an interferometer baseline. In general, when probe antennas are used 
to measure the field in a Fourier plane, these antennas receive power from the 
entire field of view. In such a case two antennas mounted very close together 
can both receive power from the entire field -of -view. The two antennas are 
similar to a monopulse pair such as is used in amplitude monopulse. Use of 
two antennas doubles signal power. Therefore, K becomes 1//2. The treatment 
of the temperature sensitivity calculation is entirely different, if there is 
not amplification ahead of the square-law detector. This is the usual state 
of affairs in infra-rod radiometers. This arrangement is also appealing in 
the microwave region in application where a very simple radiometer is desired, 
even at the expense of loss of performance. It also can become attractive 
when a great number of channels must be monitored simultaneously. 


In the usual microwave case, where there is amplification ahead of the 
square-law-detector /multiplier, the performance of a spectrometer using fil- 
ters, and a Fourier-transform spectrometer, using autocorrelation, is the 
nnrno. However, when there is no UF umpiif icution, the transform spectrometer 
temperature resolution is 


at transform 



at filter 


(1-3-9) 


where N is the number of frequency bins. The factor /2/N is called Fellgett’s 
Advantage. 

The sensitivity of transform radiometers without RF amplification is 
treated in Appendix A. 
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We will now summarize Che values of K for various types of correlating 
radiometers in Table 3-3. For comparison purposes the values of K for the 
common types of power measuring radiometers will be given in a purallel 
column. Correlating radiometers with separate antenna Inputs to the two sig- 
nal channels will be identified by an asterisk. Transform spectrometers can 
reduce K by a factor of 1//2 by measuring 4>(t) over a range of delays which 
extends from -* to +T^. The values of $(t^) are added in the in-phase chan- 
nel outputs , whilst the value of <$> (-t^) is multiplied by minus one and added 
to the value of (x^) in the quadrature channel outputs. In both cases the 
signal add amplitude-wise and noise in the two outputs adds power-wise. The 
rusuit is the Lf Jl improvement in K already referred to. 

In Table 3-3, FTS will stand for Fourier- transform spectrometer. 


TABLE 3-3. VALUES OF THE RADIOMETER FIGURE-OF-MERIT, K 


Correlating Types 

K 

Power-Measuring Types 

K 

Phase-Switched , Square-Wave 
Modulation and Demodulation 

2 

Dicke , Squarewave 
Modulation and 
Demodulation 

2 

Phase-Switched, Square-Wave 
Modulation and Demodulation* 

</l 

Dicke, Squarewave Mod Sine 
Wave Demod . 

u//2 

Homodyne (Faris ' case) , 
Total Power 

/2 

Dicke, Sine Wave Mod and 
Demod 

2^2 

Thermal-Grad ien t 


Dicke, Push-Pull (Graham's) 
Squarewave Modulation and 
Demodulation 

/I 

FTS (Positive Delay Only) 

/I 



FTS (Positive Delay Only)* 

1 

Thermal Gradient 

1 

FTS (Symmetric Delay) 

• 

1 

Total Power 

1 

Fully-Coherent 

1 



Fully-Coherent* 

1//2 



FTS (Symmetric Delay)* 





*Star indicates each correlator channel is fed by a separate antenna. 
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3.2 RECEIVER SENSITIVITY 

In order Co calculate the performance of point designs some data on 
contemporary receiver noise temperatures are needed. Two graphs are given. 
One by Weinreb (Reference 42) emphasizes FET amplifiers (Figure 3-1) . The 
other , due to Cardiasmenos (Reference 43) , shows the performance of super- 
heterodyne receivers (Figure 3-2) . This second paragraph has been redrawn, 
since the original did not reproduce well. 


3.3 GAIN FLUCTUATION NOISE 

Power ouL [in l variations duo to receiver gain fluctuation are fully as 
important as receiver noise in determining the system temperature resolution. 
To include fluctuation noise in the temperature resolution calculation, 
observe that for a total power radiometer 


P m T G 
A SYS U 


(1-3-10) 


where P is the RF power input to the correlator or square-law detector, as the 
case may be. The total differential of (1-3-10) is 


AP - + T sys SG 


(1-3-11) 


G may he considered to be a system calibration factor. So if we divide 
(1-3-11) by G, we get the correct answer TsYS " AT when G is zero. However, 
if G is not zero, division by G gives 


AP 

G 


AT + T 


SYS 


(£) 


(1-3-12) 


The second term is error due to gain fluctuation. . Since receiver noise and 
gain fluctuation are both noise-like quantities, their variances add to pro- 
duce the total error. So 



(1-3-13) 


1-3-7 
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notes: 

NOISEFIGURE 10 LOG F, AND NOISE , 

TEMPERATURE T n • 290° (P-1), VERSUS 
FREQUENCY FOR VARIOUS 1SS0 STATE -OF* 
THE-ART LOW-NOISE DEVICES. THE 300 K k 
BIPOLAR TRANSISTOR, FET. AND FARAMF 
VALUES ARE TAKEN FROM MANUFACTURERS 
DATA SHEETS (1) - (3), THE 20 K FET CURVE [ 
IS FROM THE DATA OF THIS FAFBR FL.US DATA 
OF OTHERS (4) - (61 AT 0.6, 1.4, AND 12 GHz, | 
RESPECTIVELY. THE 20 K PARAMP,4 K MASER ' 
(INCLUDING PARAMETRIC UP-CONVERTER AT 

LOWER FREQUENCIES), AND 20 K MIXER 

RESULTS (WHICH ARE SSB AND INCLUDE IF 
NOISE) ARE PROM SYSTEMS IN USE AT NATIONAL 
RAOIO ASTRONOMY OOSERVATORY (NRAO). 

THE NATURAL NOISE LIMITATIONS OUE TO 
GALACTIC NOISE, THE COSMIC OACKGROUNO 
RADIATION, AND ATMOSPHERIC NOISE ARE 
FOR OPTIMUM CONDITIONS ANO ARE TAKEN 
FROM (9) PLUS POINTS AT 22 OHZ ANO 100 GHz 
MEASUREO AT NRAO. 


FIGURE 3-1. PERFORMANCE OF LOW NOISE DEVICES 



FIGURE 3-2. PERFORMANCE OF TRG SUPERHETS 
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In a Dicke radiometer the unknown and reference temperatures ere both 
subject to the same gain fluctuation. Therefore the gain fluctuation error 
is |T S ig ” TrefI* ^ the difference between Tsyg end Tppp is servoed to zero , 
the gain fluctuation error also goes to zero. 

In the correlating radiometer case receiver noise in the two receiver 
channels does not correlate when the correlation function is generated. The 
autocorrelation functions of the receiver noises appear at DC and are 
rejected by the bandpass filter between the square law detector and the I&Q 
phase detectors. 

Of course, if the radiometer is operated in a total-power mode, the 
noise from both channels will not be rejected and will appear in the radiom- 
eter output as a pedestal upon which the signal rides. Then gain fluctuation 
will produce just as large errors as in the noncoherent, total-power radiom- 
eter case. 

In the frequency offset mode the correlating radiometer will still 
experience gain fluctuation errors due to variable amplification of the sig- 
nal Itself. However the signal is usually smaller than the receiver noise 
(e.g., 200*K vs. 1000*K) . So the gain fluctuation noise is reduced 
accordingly. 

However, it will be found when calibration is discussed, that a refer- 
ence signal can he injected into the two correlator cluinnels in u manner 
that causes the correlator output to be proportional to (TsiG “ Tapp) . 
Therefore if this difference is driven to zero, gain fluctuation noise will 
go to zero. Thin behavior is very similar to the case of u Dickc with a 
servoed reference temperature. 


1-3-9 
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4. CALIBRATION 


4.1 CLAPP CALIBRATOR 

Almost all of the calibration methods used by noncoherent radiators can 
bo adapted for use in n correlating radiometer. Uowever , a calibration method 
proposed by Clapp (Reference 44) is clearly superior to other ways of calibrat 
ing coherent radiometers. This method can only be used this radiometer type. 

A functional block diagram of this calibration method, taken from Clapp and 
Maxwell's paper is shown in Figure 4-1. 

Following Clapp, the inputs to amplifiers and G 2 are 



A 1 + (B P r + B nr + B pi + B ni> 


A 2 - (B pr “ B nr ’ * B pi + J B ni> 


(1-4-1) 



FIGURE 4*1. SYSTEM CONCEPT, COMPLEX 
* CORRELATION RADIOMETER 
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The antenna signals are 
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The 


j A l " A lr + JA li 
\ 2 " A 2r + JA 2i 

cross-correlation function F can be defined by 
P ■ V2 " F r + 


Its real and imaginary parts are given by 
* A lr A 2r + A 1 I A 2 i ^ 

<A li A 2r “ A lr A 2p 



(1-4-2) 


(1-4-3) 


(1-4-4) 


A complex error signal will be developed by the pair of product detectors 
shown in Figure 4-1. This error signal has the form 


E - (G 1 S 1 ).(G 2 S 2 )* - (G 1 G*)(F r - H r ) + j (G^*) (1-4-5) 


H is the contribution of the complex noise source whose real and imaginary 
parts are 


/H - B B* - B B* 

1 r pr pr nr nr 

(h - B B* - B B* 

\ i pl pl ni ni 

Cross product terms which average to zero have been dropped. 


(1-4-6) 


1-4-2 
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The real part of the error signal actuates a servo loop which controls 
noise sources B pr and B nr . Likewise the imaginary part of E controls B p i and 
B n i> When these servo loops have been brought to equilibrium the noise sources 
will be held at levels which can be seen to be 




. (1-4-7) 


In a correlator which measures only a single temperature only the noise 

sources B _ and B _ are needed, 
pr nr 

4.2 HACH SERVO LOOP 

Hubbard and Erickson (Reference 45) used thermionic noise diodes as 
noise sources when they mechanized the radiometer analyzed by Clapp and 
Maxwell. However, the Hach servo loop discussed in Section 2, which uses 
semiconductor noise diodes, is better for orbital use. The phase tracking 
between the two channels may be quite poor without preventing the servo loops 
from rcuching equilibrium. Let and $2 be the cumulative phase 'shifts 
through the two channels from the input antennas or input power-dividing 
hybrid. Then ($] - $ 2 ) is the relative phase shift between channels. It is 
a measure of phase tracking. 

Clapp and Maxwell computed the convergence of the servo loops and 
plotted their results on the complex plane. Their plots, showing servo con- 
vergence for relative phase shifts of +45° and -84*, are reproduced here as 
Figure 4-2. In actuality the relative phase shift will not exceed 10* to 20°, / 

for a well designed radiometer, under the worst of conditions. 

Figure 4-3 shows experimental results of Hubbard and Erickson's (Refer- 
ence 45) loop stability tests. Part (a) shows the radiometer response with 
various relative phase shifts. At +84° one can see an increase in servo 
noise. At +90* the servo breaks into oscillation as one would expect. 

Part (b) shows the result of varying servo gain ±5 db around a reference gain 
of 0*. Again the system input is 1200*K. Only slight variations in system 
response can be seen. Figure 4-4 shows receiver response vs correlated sig- 
nal input. 


1-4-3 
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NOTE! CONVEHQENCE OF COMF1. IX NOISE SOURCE, 

H, FROM INITIAL. VALUE, H 0 , TOWARD FINAL VALUE, F 


, a) THE SOLID CURVE IS THE SPIRAL FOLLOWED 
WHEN THE RELATIVE PHASE SHIFT OF THE 
TWO AMPLIFICATION CHANNELS IS «4S 9 . THE 
DASHED LINE SHOWS THE CONVERGENCE 
WHEN THE RELATIVE PHASE SHIFT IS 0° 


b) THE SOLID CURVE IS THE SPIRAL FOLLOWED 
WHEN THE RELATIVE PHASE SHIFT OF THE 
TWO AMPLIFICATION CHANNELS IS -E4°. THE 
DASHED LINE SHOWS THE CONVERGENCE 
WHEN THE RELATIVE PHASE SHIFT IS 0° 


FIGURE 4-2. SERVO CONVERGENCE 


QUADRATURE PHASE CHANNEL 



1200°K 1200°K 1200°K 1200°K 

SO 0 84° 46° 0° PHASE SHIFT 


OUAORATURE PHASE CHANNEL 


1 

a 


1200°K 1200°K 1200°K 

0 OS +S OS -8 OS 



IN PHASE CHANNEL 



1 T \ \ 

HPRRM : i immA WWFS mm 


1200°K 1200°K 1200°K 1200°K 

SO 9 S4° 46° 0° PHASE SHIFT 


IN PHASE CHANNEL 



0 DB +8 OB -8 OB 


a) RECSIVER RESPONSE VERSUS RELATIVE b> RECEIVER RESPONSE VERSUS 

AMPLIFIER PHASE SHIFT AMPLIFIER OAIN CHANGE 


FIGURE 4-3. RECEIVER RESPONSE 
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FIGURE 4-4. RECEIVER RESPONSE VERSUS 
CORRELATED SIGNAL INPUT 


4.3 REFERENCE NOISE-SOURCE STANDARDS 

In the past noise diodes and gas-discharge tubes have been popular 
standard noise sources (Figure 4-5). The thermionic noise diode will not be 
considered Further. A complete treatment of temperature-limited vacuum tube 
slot noise is in the Rad Lab Series (Reference 46). 

Gas-discharge sources iuive been widely used. They operate at frequencies 
up to tans of gigahertz. The effective temperature is very close to the kinetic 
temperature of the plasma electrons in the discharge. This is of the order 
of 10,000 # K. The most coonum gas used is argon, although neon has been used 
at times. Again these tubes belong mainly in the past. 




9.17 a 10 ## K 





* , * 

-- ■ 1 

9.17 x 10* °K 


: [ 6.79*10 ; 


1 1 
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FIGURE 4-5. SYSTEM NOISE SOURCES - TYPICAL NOISE OUTPUT 
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Ac fche present cine Che most popular standard noise source is tha 
semiconductor noise diode. Microwave SemlconducCor Corp. Is a supplier of 
these devices. «Some pertinent data follows abstracted froa their literature. 

General operating characteristics: 


Operating temperature 
Temperature sensitivity 
Voltage sensitivity 
Current sensitivity 
Switching speed 


-55* to 100*C 
0.01 db/*C 
0.1 db/% 

0.03 db/% 

1.0 ysecond 


Tlie highest frequency diode is quite flat in response. It appears chat it 
can be used out to 100 GHz or so. 


With proper temperature control and good stabilization of applied power 
the semiconductor noise diode can hold calibration for long periods of tine. 
However, calibration against an absolute standard is easier to provide than 
for a noncoherent radiometer. The reason lies in the fact that a correlating 
radiometer, using a hybrid power splitter at its input, measures the differ- 
ence between the sum and difference ports of the hybrid shown in Figure 4-6. 

When calibration is required the switch-pair are thrown to "calibrate" 
position. The servo loop then produces a frequency which sets the diode 
noise output waveform to a value such that the 200*C temperature difference 
between the 100°C and 300°C loads is balanced out by the diode noise injected 
into the mixers in the calibration system by directional couplers 3 and 4. 

When the switch-pair returns to "measure" position, the servo loop 
balances out the unknown temperature difference between the input and the 50°C 
load on the hybrid difference channel. It may be seen that the ratio of the 
unknown difference, present in the measure mode, und the known difference, 
present in the calibrate mode, is the same as the ratio between noise wave- 
form duty factors present in two codes. So the unknown temperature differ- 
ence is this ratio times the temperature difference between the calibration 
loads. Adding the temperature of input hybrid delta to the calculated differ- 
ence gives the unknown temperature. 



S'SSSSK 



FIGURE 4*6. CALIBRATION OF NOISE DIODE AGAINST AN ABSOLUTE STANDARD 

i ' 


4.4 DISCUSSION OF THE CALIBRATION CHARACTERISTICS OF FREQ-DOMAIN IMAGERS 

If the radiometer measures the autocorrelation of the signal, the use 
of a servo loop of type just described remedies a fault of autocorrelators, 
the large dynamic range they require. This results because the 2 ero differ- 
ential delay bin of the autocorrelator measures the total signal power. Uhen 
the signal has a wide bandwidth, as most radiometric signals do, the autocor- 
relation function may drop 10 or 20 db in the first couplo of delay bins 
adjacent to the x * 0 bin. As a result a large portion of the system 
dynamic range is used up just to accommodate the x ■ 0 bin. 

The servo reduces the output of the x * 0 bin to zero. This greatly 
reduces the required dynamic range. The source spectrum autocorrelation 
function is subtracted from the signal autocorrelation function not only in 
the x • 0 bin but also in all the other bins holding the autocorrelator 
output. 
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If the noise-source spectrum is flat, the autocorrelator output will have 
the form shown in (b) of Figure 4-7 . To left of (b) diagram the original 
autocorrelation is shown before subtraction of the noise-source autocorrela- 
tion. The spectrum of the unknown signal is shown above the two correlation 
functions. 

■After (b) is Fourier- transformed the known noise-source spectrum is 
added in to recover the signal spectrum. 

Very often the average shape of the unknown spectrum is known. The 
fluctuations of actual spectrum around the average may be smal compared to 
the average spectrum power. This is true in the case of temperature sounding. 
Tike spectrum calculated from the standard atmospheric temperature profile is 
the average spectrum. Its power density expressed as a temperature, is in 
the range of several hundred degrees Kelvin. The rms fluctuation of actual 
spectrums about the average spectrum never exceeds 5*K. Assume that the sys- 
tem can accommodate 3o variations in the spectrum. We find that the ratio of 
average spectrum to spectrum fluctuations is about 210*K/30* "7°, 17 ab. 

In cases of this sort we can achieve considerable improvement in system 
performance by shaping Cite noise source spectrum before its output is injected 
into the radiometer. If ill is shaped to resemble the average spectrum and 
subtracted from the signal spectrum, the residual will be much smaller than 
Che total spectrum. Thin in shown .in Figure 4-8 where the total and residual 
spectrum*. are sketched together with their correlation spectrums. , 


POWER DENSITY SPECTRUM OP SIGNAL 


| SUSTACT NOISE SOURCE 

| POWER OtNSITV EPECTRUM OF SSRVQSO NOISE SOURcIM 


POURIER TRANSFORM 


FOURIER TRANSFORM 


a) AUTOCORRELATION OR THE 
SIGNAL SPECTRUM 




n il ■ **■ 


b) AUTOCORRELATION OP SIGNAL 
SPECTRUM MINUS REFERENCE 
SPECTRUM 


FIGURE 4-7. MODIFICATION OF OUTPUT AUTOCORRELATION BY REFERENCE NOISE SPECTRUM 
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FIGURE 4*8. TOTAL AND RESIDUAL SPECTRUMS AND THEIR AUTOCORRELATION FUNCTIONS 


It may bu .seen that .since the residual spectrum has much less power in it 
thun the total spectrum, its correlation function will reach the noise level 
of the radiometer for a smaller differential delay than the value where $(x) 

Of the total spectrum reaches zero. It is obvious that correlation measure- 
ments should bo truncated for a value of t where $(t) subsides into the noise. 
Beyond that point more noise would be introduced into the transform computa- 
tion than signal data. 

Elimination of the delay bins in the to T2 regions makes the auto- 
correlator smaller and simpler. 

Elimination of the bins in the xi to X£ region also reduces the noise 
introduced into the output spectrum. This reduction can be an order of mag- 
nitude or so, since the reduction in spectrum power is an order of magnitude 
or more. The net result is that the noise variance of the system tends to be 
proportional to the real unknown spectrum, namely the residual spectrum. 
Nothing is lost. We can get the total spectrum by adding the average spec- 
trum back into the residual spectrum. We havo used our a prior knowledge 
about the spectrum to reduce noise variance. 
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Even in mapping, where absolute temperature values are not important, 
removal of the- average temperature can improve performance. A block diagram - 
of a servo loop to accomplish this is shown in Figure 4-9. 

The loop is biased so that increase in error signal reduces the injected 
reference. The composite signal formed by injecting the noise source into 
the signal channel with negative polarity will adjust itself so that a few 
negative going excursions will hold the loop in equilibrium as shown in • 
Figure 4-10. 

This arrangement permits display of the signal with maximum dynamic 
range. Dynamic range is always in short supply in display systems. Signal 
control which produces better image dynamic range always causes dramatic 
improvements in image quality. 

When precision measurements are required the elimination of the sky-horn 
or cryogenic load as the low temperature reference is important. Elimination 
of the sky-hom makes it possible to design better radiometer installations 
aboard spacecraft. When power-measuring radiometers are used, the sky-horn 


AVERAGE 
IMAGE FLUX 



RMS SIGNAL AMPLITUDE TO 
AUTOMATIC OAIN CONTROL CIRCUIT 


FIGURE 4-9. REFERENCE NOISE INJECTION SERVO LOOP FOR TERRAIN MAPPING 


▲ AMPLITUDE 



AVERAGE 
IMAGE PLUX 


' OCCASIONAL NEGATIVE 
PEAKS PROVIDE SERVO 
ERROR SIGNAL 


FIGURE 4-10. STABILIZED MAPPING VIDEO 
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oust see the sky well above the horizon In order to produce a predictable 
temperature. This often drives the designer to put the radiometer antenna 
on top of the spacecraft. Now the scanning beam must look down outboard from 
the craft — peering over the edge as it were. 

This causes restrictions on field-of-view, since the beam cannot look 
through the spacecraft. 

What is worse, the top mount causes undesired interactions between 
antenna and spacecraft structure. If one visualizes the beam as simply emerg- 
ing from the antenna and continuing to the ground, it is hard to see why the 
interaction occurs. 

However, the beam forms out of a haze of partially reactive electromag- 
netic field extending out to about d^/X. It is this haze which couples to 
the antenna surroundings. Rapid changes in the antenna pattern can make it 
very hard to correct for antenna pattern bias by deconvolution. 

It is also possible to get beam dropouts in directions where there is 
tight nntcnna/structuro coupling. As a result, most aircraft antenna installa- 
tions are checked out on model ranges in order to eliminate unpleasant 
surprises. 

The Clapp calibrator makes the sky-horn unneeded. Referring to Figure 4-1 
we see that all that is required to get a zero or negative temperature is the 
reduction of Bpr and increase in B nc until sufficient negative temperature has 
been injected into the system to produce the desired condition.' Since there 
is no sky-horn, the scanning antenna can be mounted on the bottom of the space- 
craft where field-of-view problems are reduced. Coupling of the antenna field 
to the spacecraft is also much reduced. 

There is one oilier characteristic of the correlating radiometer which is 
important enough to warrant discussion. In correlating radiometer designs 
where correlation functions are formed, this function will be Fourier trans- 
formed to obtain a spectrum. If the function is a cross correlation formed 
by aperture synthesis, the transform will be an image, not a spectrum. 

In either case, the transform calculation produces noise in the spectrum 
bins or in the Image pixels which originates in all the delay bins of the cor- 
relator. The transform process consists of multiplications of the values by 
phase shifting exponentials of the form, exp(-j2Tr<j>) . The phase shifted terms, 
when summed, form the sequence of value output by the transform algorithm. 

The noise in the input bins comes along and also is multiplied by the exponen- 
tials. However, being noise, this has no effect. The summing steps in the 
transform then adds all the noise samples power-wise. 


WHQjNAl PAGE IS 

OF POOR QUALITY 


As a result the transform algorithms acts like a food blender. All output 
resolvable elements have the same noise mix from the correlator delay bins. 

Each signal element output by the transform algorithm is a mixture of all $(x) 
values. This process ties the whole instrument together giving it one response. 
Variations in response or noise level in the autocorrelator bins do not appear 
in the spectrum or pixel bins . One calibration operation calibrates the total 
output . This ensures smooth spectrums and images devoid of streaking and 
uneven response. 

This property of the transform is so valuable it is purposely built into 
some image transmission systems. The image is transformed to the frequency 
domain before transmission and converted back on reception. If interference 
or channel drop-outs produce glitches in the spectrum during transmission, 
the result after inverse transformation on reception is a slight drop in the 
image signal-to-noise ratio distributed throughout the image. This is much 
more desirable than the defects in the image itself, if it were transmitted 
without transformation. 
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5. SPATIAL FREQUENCY CONCEPT 


5.1 INTRODUCTION 

Almost from the beginning of electronic technology amplifiers, transmission 
lines and other devices have been analyzed by calculating their behavior in 
thu frequency domain. This involves finding the spectrum of the applied sig- 
nal. Modification of the spectrum by passage through the device is then cal- 
culated. The properties of the device regarded as a filter are the character- 
istic of interest in doing this analysis. 

When thinking about antennas, the concept of the antenna as a filter is 
equally useful. In the case of the electrical filter, mentioned first, the 
signal is a function of time. It is expanded into a spectrum of temporal 
frequencies. In the case of antennas and lenses the signal is a two-dimensional 
function whose independent variables are two angles. This function can be trans- 
formed into a two-dimensional spectrum whose independent variables are two 
spatial frequencies. 


5.2 HEURISTIC DISCUSSION OF SPATIAL FREQUENCY 

The spatial frequency spectrum describes the distribution of some property 
over the field-of-view of the antcnna/lens. For instance, let the property be 
radiant brightness. Figure 5-l(a) shows a one-dimensional version of brightness 
distribution. Let us treat this distribution just as if it were a function of 
time, not angle. Let its Fourier transform be found. The new independent vari- 
able is a spatial frequency in cycles per radian or cycles per degree and the 
new function is a spatial frequency spectrum. Figure 5-1 (b) shows a sinusoidal 
component of the brightness distribution present in 5-1 (a). It corresponds to 
a small portion of the total spatial frequency spectrum which falls within a 
small frequency increment. At zero angle it may be seen that the sinusoidal wave 
representing this component is neither at zero nor at a maximum. This means 
the waveform is neither a pure sine wave nor a pure cosine wave. However, it 
can be represented as the sum of a sine wave and a cosine wave, each of proper 
amplitude. The cosine component is a part of the real spectrum. The sine 
component is part of the imaginary spectrum. If the spectrum is observed with 
a fully coherent radiometer, the real spectrum is the portion of signal which 
peases through the in-phase synchronous demodulator (phase-detector). The 
imaginary component is recovered from the total signal by the phase detector 
in quadrature with the reference. 
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FIGURE 5-1 . RELATION OF THE SPATIAL FREQUENCY 
SPECTRUM TO SCENE BRIGHTNESS DISTRIBUTION 


Note that (c) is an amplitude distribution. The radiant brightness in °K 
is proportional to rf power density (watts/Hz/steradian) . However, a cor- 
relator, like a square law detector, outputs a voltage proportional to rf 
power level. So the video amplitude is proportional radiant brightness and the 
video variance has the dimensions of degrees Kelvin squared (°K^) . This point 
is important because the amplitude spectrum is just another form of the image 
signal. An inverse transform will recover the brightness distribution. The 
power density spectrun, S(F)S*(F), has lost the phase information. So the 
scene brightness distribution cannot be recovered. In short, the power density 
spectrum is a piece of statistical information about many scenes, while the 
amplitude spectrum of a given scene is that scene in a coded form resulting 
from Fourier transformation. 


5.3 DETERMINATION OF THE POWER-DENSITY SPECTRUM OF THE SCENE 

The statistical data stored in the video power density spectrum is of 
first importance in designing radiometers and predicting their performance. 

This will become apparent when we consider frequency-domain imagers in detail. 

At present, data about radiometric -scene statistics is very sparse. In 
order to have a reasonable spectrum for use in this study, estimates of the 
video power density spectrum were made by three different methods. The details 
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of the process of making these estimates are given in Appendix C, "Microwave 
Radiometer Image Statistics." 

The first spectrum estimate is based on the work of Mandelbrot (Refer- 
ence 48). His spec Crums have the form P(F) - Const. F~2H-1. These spectrums 
are associated with scenes whose statistics are independent of scale. This is 
a property of many natural scenes. Let these scenes be observed over a range 
of magnifications. The resulting scenes are said to be self-similar with 
respect to a set of ratios r-j when the statistics of all the scenes are the 
same . 


A spectacular example of this behavior was the series of images of Jupiter 
sent back by Voyager as it closed to a shorter and shorter range. As range 
decreased Jovian weather patterns will become too big to be seen as separate 
entities in the image. However, little pattern details which were not resolved 
before now could easily be seen. The new patterns had the same overall appear- 
ance as the old patterns. 

Mandelbrot's final choice for H was arrived at by comparing the properties 
of statistically-generated landscapes with actual landscapes. These computer 
generated landscapes are shown on pp 210-215 of reference 48. His final choice 
for the value of H is 3/4. 

Another means oF determining II is based on the Fact that the spectrum 
shape is a measure of redundancy in the scene. If the spatial frequency spec- 
trum were flat, each new pixel value would be completely unpredictable from 
knowledge of past pixel values. Let the high frequencies in the spectrum be 
ruducud. As this reduction increases und the spectrum drops off more and more 
steeply redundancy increases. 

There are ways of coding the image which decrease the bandwidth required 
to traimiuit images at a given rate by removing the image redundancy. The 
amount of image compression obtainable depends on the redundancy of the original 
image. This in turn depends on the value of H in the general spectrum given 
for P(F) . 

In Appendix C.2 it was estimated that satisfactory image compression 
codes can reduce images with a dynamic range of 1,000 to formats with an 
information density of one bit per pixel. The value of H which agrees with 
these facts was then found to be 0.70. 

A third approach used experimental microwave radiometer spectrum data 
supplied by Mr. Charles Hawthorne of the Navy Ordnance Test Station, Inyokern. 

A regression algorithm was used to fit the experimental data to the power-law 
spectrum form. The value of H obtained in this manner was 0.60. Figure 5-2 
shows the raw spectrum data. 

Extension to the actual two-dimensional scene spectrum is made by 
expressing the spectrum in polar coordinates. The one-dimensional spectrum 
is then used as the spectrum function along coordinate radials. No variation 
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FIGURE 5-2. POWER SPECTRUM VERSUS FREQUENCY 


In spectrum value occurs in angle. A value of H - 3/4 was adopted for this 
study. 


5.4 ALTERNATE STATISTICAL DESCRIPTIONS 


Other power-density spec tr urns, as well as alternate statistical formats, 
have been used to describe scene statistics. They are reviewed in this section. 

Viterbi and Omara (Reference 50) treat the case where the scene is iso- 
tropic and has a space correlation function which depends only on the total 
distance between two points in the two-dimensional space. That is 


V(r ,r ) 
T x y 



(1-5-1) 


Changing to polar coordinates, let 

V" rcos6 r * r y " rsin6 r 

and 

u ■ wcosO, . u> - usinfl 
x m y a) 
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Then Che polar spec Crum is 
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Vicerbi end Omara also report a spectrum used by Sakrison and Algazi for 
television image analysis. It is 


2(u 


+ (<i>) - 


2 2 
uT + ID 

o 


(1-5-3) 


The corresponding autocorrelation is 
♦(r) - exp(-[r|d) 

where d is the coherence length of the field. 
Berger (Reference 51) uses the spectrum 


(1-5-4) 
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where A is chosen so that (2v) 
n * 
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Pratt (Reference 52) assumes a Markovian source (memoryless source) whose 
elements have equal variance, o^ and which possess adjaceut pixel correlation 
factors of and p c along rows and columns. 

5.5 SPATIAL FREQUENCY VERSUS TEMPORAL FREQUENCY FILTERS 

* • It should be obvious at this point that an antenna acts as a filter of 
spatial frequency spectrums in the same manner as electrical filters operate on 
temporal frequency spectrums. 

A detailed development of this parallel between temporal and spatial 
filters is given by O'Neill (Reference 53). Their parallel properties are 
ttliown in O'No ill's Table 2-1, which in reproduced here as Table 5-1. The 
quantities in Table 5-1 are defined in Table 5-2. 


TABLE 5-1. SPATIAL FILTERS VERSUS TEMPORAL FILTERS 



Time filter* 

Space filteni 
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TABLE 5-2. DEFINITIONS OF FILTER FUNCTIONS 
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5.6 TREATMENT OF RANDOM WAVEFORMS 

For a random signal existing over an infinite interval, the Fourier 
transform does not exist. Therefore, a statistical description is required. 

In general, the analytic attack consists of defining the operation of a finite 
interval. The statistical functions is then the limiting value of the output 
of .the operation as the interval becomes infinitely long. The expressions for 
Random Input in Table 5-1 stem from this rationale. For example, we have 


♦ «(p) ■ <o(p')o(p+p')> » lim -~ 

00 T-h. 4L 2 


This is mathematically correct, however, there Is a problem. Since we cannot 
in uctuuiity averuge over a record of infinite length, the equation (1-5-6) 
really instructs us to calculate the mean values over a finite interval T or L. 
Then we obtain an estimate of the spectrum by taking the Fourier transform of 
the correlation and letting the interval become indefinitely long, or at least 
long enough to ensure that the measured S(u>) will meet accuracy requirements. 
Thus 


+L +L 


a 


o(p)o(p.+p’)dp' 


*«( T ) " 

XX 


<X(t)X(t+r)> - 


lim 

T-*» 


+T 

•h J 1 

-T 


X(t)X(t+r)dt 


S(uO t" s t ( w) M T > = <X(t)X(t+r)>_, (1-5-7) 

X ” H * T-h» T**» 1 

Unfortunately, the variance of S (w), in the limit, is very large (Refer- 
ence 54). Swartz cites Jenkins and Watts (Reference 55) for proof that 


lim Var 8 (m) - S 2 («) (1-5-8) 

T-MO i 

Using the ratio of mean to standard deviation as a measure of signal- to- . 
interference ratio, we find that 

E[S _(«)] 

> tpsr ^ 1 (1 - 5 - 9) 

This says that the true spectrum is only as big as the noise (uncertainty in 
S t (m)). 
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The way out of this difficulty is to chop the signal record up into a 
set of segments. The Fourier transform of each segment is calculated individually 
to get a set of spectrum estimates, S n (u>). The N estimates, S n (u>) are then 
averaged together to get the final estimate of S(u). 

* 

Now in practice the total observation time, T, which is available is 
not Infinite. The question then arises, "How many segments of duration T/N 
should we use?" As N gets larger and larger, the variance in the estimate of 
S (w) decreases due to better averaging. However, the resolution of our esti- 
mator becomes poorer, because the length of each observation segment is getting 
shorter. This increases the mean square difference between the true spectrum 
and estimated spectrum due to loss of resolution. This error is called bias. 

It should be evident that there is an optimum N that produces the m-tn-iimim 
total error. This total error is noise variance plus bias squared. We will 
drop this subject for the time being. However, we will see as we progress 
chat preserving a balance between variance and bias is a central theme in 
analysis and a design driver in mechanization. 

The important point at this time is the observation that each signal 
record is a transient. Therefore, the relations for transient inputs in 
Table 5-1 can be used. The Random Input relations in the table are no longer 
mathematical identities, but approximations. Their continued use is important 
because they eliminate all phase information concerning the successive random 
signal records. These functions are therefore needed to analyze apd design 
system configurations .which must operate in an optimum manner when driven 
by various members of a class random signals. In contrast to this state of 
affairs, functions such as x(t) and X(u) or y(t) and Y(w) are the signals 
themselves. 


5.7 INTERRELATIONS BETWEEN FUNCTIONS IN TABLE 5-1 


It may be seen that there are three classes of functions involved for 
both time and space filters. One class concerns the input signals. A second 
function cIumh concerns tiic output signals. The third class involves the signal 
transformation process produced by the filters. Stripping away the details of 
limits, normalization and the like. Tables 5-3, 5-4 and 5-5 show these rela- 
tions. A star stAnds for convolution. Dual stars indicate a two-dimensional 
convolution. Likewise, a dot means multiply. Two dots mean carry out a 
two-dimensional multiplication. The symbols and Indicate one and two 
dimensional Fourier transforms, respectively. Likewise, single and dual 
pentagons, ★ and ** signify single or dual autocorrelations or cross corre- 
lations as the case may be. 

The relation, v™* <I>(u) Is an expression of the Weiner-Xhinchin 
theorem. This theorem states that the Fourier transform of the signal auto- 
correlation function is the power-density spectrum of the signals. Convolutions 
in the space /time domain are always paralleled by products in the frequency 
domain as a result of the convolution theorem. 
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TABLE 5-3.* FILTER INPUT RELATIONS 
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TABLE 5-4. FILTER OUTPUT RELATIONS 
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Space Filters 
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TABLE 5-5. FILTER INPUT/0UTH7T RELATIONS 
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5.8 EYE RESPONSE TO THE SPATIAL FREQUENCY SPECTRUM 

» 

Just as in the case of temporal filters, one of the great advantages of 
analysis in the frequency domain is the fact that a series of modulation 
transfer functions can be multiplied together to obtain the MTF of the total 
system. A good introduction to this technique is afforded by examination of 
response of the human eye to the spatial frequency spectrum. This will be 
followed by a description of the calculation of eye response to a complete 
photographic system (Reference 56). 

Unlike almost all other optical elements, the eye MTF is that of a band- 
pass filter.. Figure 5-3, taken from Price’s paper referenced above, shows this 
response. The low spatial- frequency spectrum components are removed by the 
retinal neural network which processes the output of the light sensitive elements 
in the retina (rods and cones). Apparently the reception of the low-frequency 
components of the scene viewed by the eye had little survival value (or even 
an udverse effect). These components, therefore, were eliminated by evolu- 
tionary processes. 

It has been found that optimum photographic systems have a total MTF, 
including the eye MTFj which has the greatest area under the composite response 
curve (Reference 56). This condition is subject to various constraints such as 
constant available light flux and rules for selection of viewing distance. 

Figure 5-4 is taken from Price’s paper. T.t shows how the total MTF can be 
found. • 

•The use of this technique in order to generate radiometer images for 
photointerpretation will be taken up later in this study. 


5.9 CONCLUSIONS 

The high degree of parallelism between temporal and spatial filters 
would indicate that the concept of spatial frequencies and spatial frequency 
spectrums should be very important when designing lens systems or antenna systems. 
This indeed is true in optics. The modulation transfer functions of lenses 
and other optical components has become the standard method of specifying their 
performance. Design of lenses by use of frequency-domain concepts have become 
the dominant design methods. Recursive algorithms for computer lens design use 
very sophisticated descriptions of the response of the human eye to determine ..... 
what changes in the lens design should be made on each recursion (Reference 56). 

Radio- astronomers have made extensive use of spatial-frequency concepts 
in order to design radio telescopes. Unfortunately, most antenna engineers 
have not followed suit. Often this has led to designs which appear primitive 
compared to those of optical engineers and radio astronomers. This study will 
concern Itself almost entirely with these frequency-domain design techniques 
as a means of designing radiometric mappers . This concluding section seems to 
be a good place to mention one other interesting correspondence between spatial 
and temporal frequencies. In addition, some differences between space and time 
filters will be mentioned which are not apparent from the tables. 
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In optical/antenna systems light distributions fall on both sides of the 
coordinate axes. Indeed, these distributions often are symmetric around these 
axes. This means that the causality condition, which exists when t is the 
independent variable, does not exist in the spatial frequency domain. There- 
fore, the real and imaginary parts of the transfer function, t(u.,w ), are not 
a Hilbert transform pair as in the temporal case. * n 

Also, in general, in the temporal case frequency is a complex variable. 
That is S ■ o + iw. Then h(t) and H(u) are related by a generalization of 
the Fourier transform called the Laplace transform. Spatial frequencies are 
always real, so the Laplace transform is never required. 

Neither of these differences harm the usefulness of the correspondence 
between temporal and spatial frequency functions in aiding the analysis and 
design of optical and antenna systems. 

There is an interesting correspondence between the uncertainty relations 
in both domains. We all know that for temporal frequencies the uncertainty 
relation is 

~ AwAr >_ 1 (1-5-10) 

The corresponding relation for spatial frequencies is (Reference 57) 

~ ApAx 1 1 ' (1-5-11) 

where Ax is the size of a pixel resolved by the optical spatial channel and 
Ap is the spatial bandwidth in radians per unit distance. 
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6. ANTENNAS AS SPATIAL FREQUENCY FILTERS 


6.1 RELATION OF ANTENNA MTF TO APERTURE EXCITATION 

In Section 5 the concept of the modulation transfer function (MTF) of a 
lens or antenna was introduced. It was shown in Table 5-5 "Input-Output 
Relations" as the Fourier transform of the filter memory function in the 
temporal case. It also was shown as the Fourier transform of a spread func- 
tion in the space-filter case. If the device is an antenna, the spread function 
is called the antenna power pattern, P(4>). 

It has been shown by Booker and Clemmov (Reference 59) that the antenna 
power pattern is proportional to the complex autocorrelation function of the 
aperture distribution function. Thus 

« 

t(x^) « E(x A )*E*(x A ) (1-6-1) 

where t(x^ Q ) is Fourier transform of the antenna power pattern, E(x^) equals 
the aperture distribution, x^ ■ x/A * distance across the aperture face in 
wavelengths. 

One would have expected Eq. (1-6-1) to be a convolution, rather than auto- 
correlation. Kraus (Reference 59) shows that when Eq. (1-6-1) is used, the 
antenna response to a brightness distribution function is 

S(* c ) - B(*)*P(4») (1-6-2) 

f\j 

where P(4>) is a mirror image of the power pattern P($). This is equivalent to 
the space filter relation in Table 5-5. However, an angular response is found, 
instead of a response over an image plane. This mirror image pattern is used 
so that the characteristic folding and inversion operation required by con- 
volution will not be needed in Eq. (1-6-1). The autocorrelation replacing 
the convolution is more difficult. Table 6-1 shows two tableau. The first 
shows where Eq. (1-6-1) comes from. The second justifies the corresponding 
convolutional form of Eq. (1-6-1). ^he 
function formed by inversion; e.g., P(4> 
in Eq. (1-6-1) is a spatial frequency. 


hat symbol ^ indicates a mirror- image 
- ♦) - P(4> - to) • The variable 
°It will be replaced by the symbol, F. 
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TABLE 6-1. CALCULATION OF THE MODULATION TRANSFER FUNCTION 
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The correlation theorem is not as well known" as the convolution theorem. 

It is discussed by Brigham (Reference 60) , If E is symmetric, the distribution 
between autocorrelation and convolution disappears. 


6.2 ANTENNA MTF FOR SOME SPECIAL EXCITATION FUNCTIONS 

Figure 6-1 shows a graphic method of finding the MTF when the excitation 
is given. All autocorrelations slope off to zero in manner shown. The equi- 
valent of the general tableau shown in Table 6-1 is given in Figure 6^2 for 
uniform weighting'. & 
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1/ax - ANCLE SUBTENDED BY 1 CYCLE OF 
CUTOFF SPATIAL FREQ 

FIGURE 6-1. GRAPHIC CALCULATION OF 
AMTF 



FIGURE 6-2. FUNCTIONS ASSOCIATED WITH UNIFORM 
APERTURE EXCITATION 


Note that TRI(F) reaches zero at +a. This is called the cutoff spatial 
frequency. The MTF is zero outside this interval. As noted, the null- to -null 
beamwidth of the power pattern, sine 2$ , is 2/ax radians. Since the period 
of the cutoff frequency is 1/a^, the cutoff wavelength executes two cycles 
per null-to-null beamwidth. The Nyquist criterion then says that if the 
antenna beam is scanned, at least four samples of the antenna output must be 
taken during the time the beam is scanned through an angle equal to the null- 
to-null beamwidth of the beam. 
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It can be seen from Figure 6-1 that the cutoff spatial frequency will 
always be equal to the aperture size in wavelengths, a^. As an example of a 
more complicated excitation function, let ECx*) ■ TRI(xx). He have found in 
Figure 6-2 that; 


, TRI(y) -w sinc 2 (z). 

# 

So the tableau for this case must be as shown in Figure 6-3. Tables of 
modulation transfer functions for more complicated types of aperture excitation 
weighting have been tabulated and are available (Reference 62) . 

All filled apertures act like bandpass filters with response which goes 
to zero at +ax* It is obvious that in order to pass the spatial frequency 
spectrum with no distortion from zero out to the spatial frequency a^, the 
antenna MTF should be rectangular. If the total spectrum of the scene fell in 
this 0 to a^ range, there would be no signal distortion. In the more likely 
case where the spectrum extends beyond a^ the rectangular MTF produces the 
least distortion possible without increasing system resolution. Bracewell 
calls the resulting image signal the principal solution. 

At first blush it would appear that regardless of what aperture weighting 
is chosen, the MTF of an aperture with continuous excitation must slope down- 
ward until a value of zero is reached at a^. Probably the closest approach 
to a rectangular response that can be obtained is produced by use pf Taylor 
weighting with fairly high sidelobes. For instance, the MTF of Taylor weighted 
antenna with -5 dB sidelobtes shown in Figure 6-4 comes as close to a rectangular 
respousu us is posslbiu (Keiereuce 62) . 



FIGURE 6-3. FUNCTIONS ASSOCIATED WITH TRIANGULAR 
APERTURE EXCITATION 
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FIGURE 6-4. MODULATION TRANSFER FUNCTIONS OF TAYLOR WEIGHTED APERTURES 


6.3 MTF OF THE COINCIDENT ARRAY OF UTUKURI AND McPHIE 

The trouble lies in the power pattern of an antenna with a rectangular 
MTF. This principal power pattern is the Fourier transform of the rectangular 
MTF. This power pattern is a sine $ response. It has both positive and 
negative sidelobes. However, a radiometer using a square-law detector cannot 
be used to obtain the princiapl power pattern because the detector output 
cannot be negative. 

However Utukuri and McPhie get around this difficulty by the use of two 
interleaved array structures. Instead of using a square- law detector, the 
two array outputs are multiplied together. This arrangement is shown in 
Figure 6-5 (Reference 63) . 

It can be seen that Utukuri and McPhie have replaced the power-measuring 
radiometer associated with a single antenna with a correlating radiometer. 

The two radiometer channels are fed by the pair of interleaved arrays. 

Figure 6-6 shows the calculated response of a five-element array. The 
sidelobes of alternating polarity may be seen. 
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FIGURE 6-5. COINCIDENT ARRAY OF UTUKURI AND MACPHIE 



FIGURE 6*6. FACTOR PATTERNS AND PRINCIPAL POWER PATTERN 
OF A FIVE ELEMENT COINCIDENT ARRAY 



6.4 MTF OF A MULTIPLYING INTERFEROMETER 

We are now ready to go on to the graphic calculation of the modulation 
tranafer function of a multiplying interferometer. Figure 6-7 showing this 
graphic solution is due to Kraus (Reference 59 ) , and is his Figure 6-22. 


1 — 6—6 
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FIGURE 6-7. GRAPHIC DETERMINATION 
OF THE MODULATION TRANSFER OF A 
MULTIPLYING INTERFEROMETER 


Thus we see that Che MTF of a multiplying interferometer has the spatial 
frequency response of a bandpass filter. 

Kraus in his Eq. (6-100), pg 184 gives the interferometer power pattern 
function for the euso where the antenna elements at the ends of the interferom- 
eter ere 'dissimilar. It is ’ 

P(4») - E^($)E 2 (iJ>)e + "^ + E^($)E 2 ($)e”*^ (1-6-3) 

This refill 1 1 In proved in Appendix li "Cunttruiixcd l.ntorforouieCcr Mapping Uchuvior." 
Eq. (1-6-3) of this section is given in Eqs. (B-15) and B-17) of Appendix B in 
terms of the gain functions of the two antenna elements with their phase 
centers moved to the coordinate origin at the midpoint of the baseline. The 
functions E^(^) and E2(4) Are the field patterns of the two Interferometer 
elements and tji - 2irS^sin$. 

If the elements are identical, that la E^($) ^ £ 3 ( 4 ) ■ E 0 (6), then the 
power pattern becomes 


P(4») - |E o (4») | 2 cos(27fS x sin«|») 


(1-6-4) 


Finally, if the element patterns are symmetrical, so that E^($) » Ej($) and 
E 2 Cd >5 ■ £ 3 ( 4 *) « the relative power pattern becomes 


P($) • E 1 (4)E 2 (^)cosi/i 


(1-6-5) 


The power pattern is seen to be proportional to the product of the field 
patterns of the Individual elements. 
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6.5 METHODS OF VARYING THE EFFECTIVE BASELINE LENGTH 

These mapping radiometers must vary the fringe pattern given by Eq. (1-6-4) 
in order to carry out their mapping functions. This can be done in a variety 
of ways. One method of changing the fringe pattern is by changing the baseline 
length, S^, by moving the antenna elements. We can also change S^ for a fixed 
antenna separation by sweeping the frequency. Both of these methods change 
the spacing between lobe peaka, as well as lobe position. 

If a phase shifter is Inserted in series with one element, ^ becomes 
(2irSxsin$ + $ s ) . If is continually increased by the phase shifter, a source 
at fixed position will produce a sinusoidal voltage at the interferometer output. 
Alturuutcly we cun put u variable length transmission line in series with one 
of the elements. Then 


i|< ■ 2n~ (Ssiwj) + Afc/p) (1-6-6) 

c 

Where Li ■ differential length of the transmission line feeds from the 
radiometer, p » transmission line phase velocity ■ v/c, v • velocity of a signal 
along the line. 

The last two methods are equivalent if Li is small. However, they behave 

t*»tl roly d I ITorciitly wIhmi M In large. 

l 

If the Interferometer elements were points, then the bandpass filter 
response would become an impulse. However, if the elements are not points, 
then the MTF passband has a finite width as shown in Figure 6-7. This will 
produce a bias error whan measuring the magnitude of the spatial frequency 
spectrum at frequency S^. Correction of this bias error will be considered 
in a la bur suction. 
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7. FREQUENCY DOMAIN IMAGERS 


7.1 FOURIER OPTICAL SYSTEMS 

Fourier optics has become an important subject in recent years. 
(References 64, 65) These techniques Involve conversion of optical Images 
to spatial frequency spcctrums, recording such spoctrums in the form of photo- 
graphic records called holograms, and filtering these spectrum records. The 
results of these operations are used to reconstruct modified images by trans- 
formation back to the space domain. 

Usually the objects dealt with in optics are recorded as images on photo- 
graphic film. Manipulation of the image starts with the photo record. In our 
case, the input signal is in the form of microwaves radiated by a scene at a 
great distance from the radiometer. We want, to record data about the radia- 
tion over a plane corresponding to one surface of an imaging lens as shown in 
Figure 7-1. 

In the optics case the image on the photo transparency is converted into 
a spectrum by use of a lens. Use is made of the fact that lenses in the 
correct optical configuration can act as two-dimensional Fourier transform 
devices. The transformed image appears as a spectrum in a plane which is 
called a frequency plane or Fourier plane. 



FIGURE 7-1. DIRECT IMAGING SYSTEM USING A LENS CONTRASTED WITH FREQUENCY-DOMAIN IMAGER 
USING TWO RAOIATION-FIELO PROBES 
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' In our case Che distant scene produces a frequency- domain field at the 
surface of the imaging lens. This surface, shown in Figure 7-1, is a frequency 
plane. The lens acts as a Fourier transformer which converts the data to the 
space domain. This converted signal is Che image which appears in the image 
plane. 

1 ’ We intend to record the frequency-plane data. Then the image can be 
generated by carrying out the Fourier transform by data processing in a com- 
puter. ‘I'lict question then a rises ns to exactly what format should he used Co 
record the frequency plane data. ' SAR radars give some guidance on this point. 


f 


HYNTIIF.T IC-AI'HU'i'URIi KAUAKS AS ACTIVE FREQUENCY-DOMAIN IMAGERS 


SAR radars are active frequency-domain imagers. The coherent video, 
which is recorded as the synthetic aperture is generated, gives us a clue to 
required format in che passive case. The SAR antenna receives the radiation 
reflected by the scene. The phase change from pulse to pulse of the signal 
return is the most important part of the data about the radar return. In order 
to tie the data from successive returns into a single unbroken record of the 
phase changes there must be a master timer signal. This signal is generated 
by very stable crystal oscillators. It is called the reference signal. Its 
smooth, predictable phase change is compared with the unknown phase change of 
Che return by limans of a phase detector (synchronous demodulator). This is 
shown in Figure 7-2. 
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< The coherent video record is Fourier transformed to get the image or 
Fresnel transformed when the scene is in the near-field of the synthetic 
aperture (the Fresnel transform is a generalization of the Fourier transform). 
The radar signal is linear in amplitude , while the radiometer signal is linear 
In power.. However, this is not important in making this comparison. The 
mathematics is the same. 


7.3 TIIIJ VAN CITTERT-ZERNIKE THEOREM 

The radiometer imager in Figure 7-1 has no internal reference like the 
SAR. However, nl.1 we need is n 'reference with a stable phase dcrivod from 
Che radiation from the scene. 

This can be provided by a second probe with a fixed location in the 
frequency plane. The fact that the probe receives noise,, rather than a sine 
wave, is not important. As long as the noise bandwidth is small compared to 
the center frequency of the noise spectrum a phase comparison can be made. 

The generalization of the phase detector, when noise-like signals are used, 
is the correlator. There the outputs of the reference and traveling probes 
are multiplied together and averaged. 

We readily will recognize the output of the correlator as a lateral 
coherence function, (tjr, r^) . Its Importance stems from the fact that 

it is the Fourier transform of the object intensity distribution shown in 
Table 5-3. This relation 1 b the Van Cittert-Zernike Theorem. This theorem 
is the basis for the analysis of the propagation of partially coherent radia- 
tion. A proof of this theorem, abstracted from Horn and Wolf (Reference 66), 
is given os Appendix E. 

This theorem says that if the linear dimensions of the source and the 
distance between the two frequency-plane probes are small compared to the dis- 
tance between the frequency plane and the source, the degree of coherency between 
these probes, r., (?», x n ), is equal to the absolute value of the Fourier 
transform of the intensity function of the source. 

In Table 5-1 relations between functions in the object plane, in the 
image plane, and relations between these function sets were given. We now 
need to prepare a similar tableau for object-plane/frequency-plane relations. 

This can be done with the aid of the Van Cittert-Zernike theorem. This is 
done in Figure 7-3. The symbol «■ a + indicates that the two functions differ 
only by a constant factor. 
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FIGURE 7*3. SOURCE/FREQUENCY PLANE RELATIONS 


7. A DUAL INTERPRETATION OF FREQUENCY-PLANE INTERFEROMETER MEASUREMENTS 

Since, o (£, n) is the transform of Fl2 and d (w*?, m n ) is a transform of 
n (C, il), d (“>r. Hi,,) muHt he proportional to l'j.2* *W»e spatiul frequency 
spectrum S (u, ‘v) is not the transform of because it is not a ‘spectrum 
, description of the variations of as and are varied. As shown in 
Section V, S (o, v) is the angular spectrum of the incoming radiation. In 
Figure 7-4 « geometric construction shows the d (u>£, w n ) and S (o, v) are 
proportional. Therefore. r i2 t^) and S (u, v) are proportional. 
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Figure 7-4 shows the two probe antennas separated by a baseline of 
length Z. This baseline length produces a set of lobes with angular frequency 
of vl (cycles per radian). There will be a component of the scene brightness 
which fluctuates at a linear frequency of u p 1 cycles per unit length. For an 
interferometer distance from the scene of OT; the interferometer lobes are 
synchronous with the brightness component of frequency A large output 

of magnitude S (u% o) or d (uA o) appears at the interferometer output port. 
All other components of the object brightness distribution average to zero at 
the output port. This occurs because the interferometer lobes are not 
synchronous with any of these other spectral components. 

We see than that the correlation function measured by the frequency plane 
probes lias u .dual interpretation. It is cither the lateral coherence function 
of the frequency plane or the spatial frequency spectrum of the scene. In 
either case Z must be varied through a range of values In order to measure the 
total function. These arguments prove the validity of the bottom and right- 
hand proportionalities in Figure 7-3. 

It is now obvious that by moving a probe over a region in the frequency 
plane, and measuring as we move, we are in a position to find the object 
intensity distribution. This is done by taking the two-dimensional Fourier 
transform of the measured lateral coherence function, F^ (Try r n ). This 
process is called aperture synthesis. The probe/correlator device, plus the 
Fourier transformer, constitutes a frequency-domain imager, 

7.5 RECIPROCAL RELATION BETWEEN. SPACE-AND FREQUENCY-DOMAIN IMAGERS 

Transform pairs always exhibit a reciprocity in width. As one member 
of the pair becomes smaller the other member becomes larger and vice versa. 

The reciprocal width theorem states thr.s property more exactly. It says that 
the equivalent width of a function is equal to the reciprocal of the equiva- 
lent width of Its transform (Reference 67) . . 

If the frequency-domain imager is contrasted with a pencil-beam scanner 
(space-domain imager), there are two transform pairs whose comparison is of 
interest. One pair is probe travel required to produce the image. The other 
pair is antenna size. We will make the comparison using a microwave lens in 
the pencil-beam, scanner case. This is not a restriction on bhe relations 
involved. The lens merely eliminates the folded ray paths produced by 
reflector antennas. This makes for a neater diagram. Figure 7-5 shows the 
two cases. 

It is Ko (Reference 68) who gives the relation D ■ A0” 1 , where D is the 
diameter of the disk sampled in the Fourier plane in wavelengths and A6 is 
angular resolution* He also states that in order to meet the Nyquist sample 
crit erion, so there is no aliasing, the sample spacing in wavelengths must be 
However, in order for the probe anten na to see the entire scene, it 
must have an aperture of no greeter then F6V“^ wavelengths. Therefore the 
probe moves by its aperture width, between successive samples. 
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FIGURE 7-5. COMPARISON OF SPACE AND FREQUENCY-DOMAIN IMAGERS 


It may be seen Chat Che principle of transform reciprocity favors the 
frequency-domain imager. This device requires a small antenna and large probe 
motion. A mapper where only large motions are needed to obtain high resolution 
is ideal for space applications. There is a lot of real estate in orbit. 
Motion over the frequency plane is limited only by the ingenuity of antenna 
engineer and the vehicle designer. 
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7.6 COMPARISON OF THE SENSITIVITY OF SPACE-AND FREQUENCY-DOMAIN IMAGERS 

in Che Introduction heuristic arguments were given to prove that the 
sensitivities of these transform twins were the same. The time has come to 
p^ove this on a more formal basis. Rather than dealing with brightness, as 
we have up to now. let us consider the total brightness of the scene when 
viewed at the location of the radiometer antenna. Brightness is a spectral 
density. Total brightness. B (6. 0 is the integral of the brightness over 
the rail I omo ter operating band. It haw the dimensions of watts iif 2 rud” . 

Since we are not dealing with a particular scene, we will replace B(8, 0 
with an average total brightness, B, over the field-of-vlew of the antennas. 

The antenna has u capture cross section A (0, 0 . The incremental power 
received from direction (8, <0 *•» & (0» 0 A (6, 0 dfl, where dfl is a.differ- 
entifij. solid angle. This cross section can be related to the antenna power 
pattern by the equation 

A (8, 0 - A e P N (8, 0 (1-7-1) 


This normalized power pattern is dimensionless. Its peak value is unity. A fi 
is the antenna effective aperture for radiation angle-of-arrlval in the 
direction of peak gain. 

t 

A very important property of antennas is the fact that the integral of 
A (8, 0 over a sphere surrounding any antenna is always X 2 (Reference 69). 
This is true regardless of the antenna size. 

Let us assume that the beam efficiencies of both the pencil-beam scanning 
antenna and the frequency-plane probe are 100 percent. Then in the first case 
the integral ol* A (0, 0 over the punei I-beam will be A 2 , In the second case 
Che integral of A (0, 0 over the field-of-vlew seen by the frequency-plane 
probe will be X 2 . The solid angles corresponding to these cases will be 
designated flj, and flj. 

The brightness of any emitting surface is given by the relation 


B (8, 0 - — (0 ; +> e (6, 0 

X 2 


(1-7-2) 
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where T is the surface temperature and e is the emissivity. The average 
brightness is given by 

¥ - f f T (G, <J>) «dA - (1-7-3) 

a * , 

scene 


where T is the average radiant temperature of the scene. 

Thun Tot l><> tli the puncil-buaiu antcnnu und fruqueney-plune probe cases 
the thermal power received is 


1 or 2 


' \ ; f l ' 

n, or R 0 


(0, 4>) A (6, ♦) 


jo a 

dn ■ 7 


// 

fl, or n. 


A (0, <j>) dfl 



/ 2kTAv \ , 2 

\~irr ' 


kTAv 


(1-7-4) 


The factor 1/2 appears because the antennas can receive only one of two 
orthogonal polarization components of the radiation. 


It may be seen that the large pencil-beam, scanning antenna and the small 
frequency-plane probe receive the same power. In actuality the frequency- 
plane imager receives twice the power ■ found in liq. (1-7-4) . This occurs 
because there are two frequency-plane probes. However, the number of samples 
in the Fourier plane is double the number of pixels in the image scanned by 
the pencil beam. Call this number N. Then the 2N spectrum samples can be N 
complex numbers. The total power distributed between the real and imaginary 
parts of the complex number is ?2 Eq. (1-7-4). However both real and 
imaginary bins contain the same noise power as that in a pixel bin. This 
results because there Is the same time available fur measurement of the' N 
samples or N pixels. Therefore the sample noise power is twice the pixel 
noise power. But the two Fourier-plane probes receive twice the signal power 
compared to the power received from a pixel. 

The final conclusion is that the space-domain and frequency-domain 
imagers have equal sensitivities when all system parameters are the same, 
in both cases. 
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1 

7 ‘,7 SAMPLING THE FREQUENCY PLANE 

Up co this point Che sampling procedure followed by Che craveling 
frequency-plane probe has been vague. It has only been specified chac Che 
sample grid be dense enough Co meec Che Nyquist cricerion. Usually Che sam- • 
Pling sequence is divided inCo cwo seeps. The firsC seep consiscs of moving 
Che craveling probe out along a baseline, scarcing at the reference probe and 
ter m i n ating at some maximum baseline length. A new baseline orientation Is 
Chen selected by stepping the baseline drive mechanism through a small. angular 
increment. Probe motion along the baseline is then repeated and another step 
in baseline orientation is taken. These motions continue until the baseline 
has rotated 180° . The final result is the sample grid in the Fourier plane 
shown in Figure 7-6. The half-wave spacing indicated ensures that the Nyquist 
criterion i - . 'aet. 

The r -r- ',ess of sampling along a baseline is often called "aperture- 
synthesis.. The terms "rotation synthesis" or "supersynthesis" are reserved 
for the e impling procedure using the rotating baseline. In the literature 
these o^t. rations are treated merely as successive portions of the Fourier- 
plane imaging process. 



FIGURE 7-6. STRUCTURE OF FREQUENCY PLANE IMAGER SAMPLE GRID 
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8. BASIC APERTURE SYNTHESIS SYSTEMS 


8.1 RESPONSE OP AN INTERFEROMETER TO AN EXTENDED, NONCOHERENT SOURCE 

This section is not intended to be a mathematical analysis in depth of 
aperture-synthesis, s Luce this has been done many times before (Reference 70, 
71). The Intent is to make a more heuristic examination backed by appropriate 
references. Comment will also be made about the applicability of the technique 
to the problem of orbital radiometric mapping. 

The response of the aperture-synthesis interferometer to an extended 
source is analyzed in Appendix B. There it is shown that only an incoherent 
source produces a linear interferometer response. Eq. (B-16) of Appendix B 
shows this response, P(X), to the power distribution over the scene, P Q (X), is 

P(X) - ?P o (X)*K(X) (1-8-1) 

where K(X) is the antenna power pattern of the interferometer. Then Eq. (B-20) 
shows the relation between the space-domain and frequency domain functions 
first introduced in Section 6, Antennas as Spatial-Frequency Filters. Namely, 


P(X) - P (X) * K(X) 

I Jt J1 «- 8 - 2) 

S(F) - 0(F) • MTF(F) 

where S(F) and 0(F) are the complex spatial frequency spectrums of the image 
and object. 

Let G;l(X) and G 2 (X) be the antenna field patterns of the two interferom- 
eter antennas. Let Ai(n) and A 2 (n) be the excitation functions of the probe 
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ape rtures . Then Eq. (B-21) of Appendix B shove the relation of K(X) and 
MTF(F) to these four functions as 

K(X) - G 1 (X) • G* (X) 

. . It U U (I “ 8 “ 3) 

MTF(F) - A.(n) * A*(n) 

t 

* 

As explained earlier , Kraus uses the correlation theorem to replace the 
convolution by a correlation to 'obtain 

K<X) » G-^X) • G 2 *(X) 

\t It . (1-8-4) 

MEF(F) - A 1 (n) ★ a 2 * (n) 

/ , 

Graphic construction of the MXF per Eq. (1-8-4) was shown in Figure 6-7, .< 

and is repeated here for easy reference. 

Kow let us vary in Figure 8-1 from zero to a maximum value D. From 
Eq. (1-8-2) it may be seen that the interferometer acts as spectrum analyzer as 
it scans over the spatial frequency spectrum of the scene 0(F). The output is 
the image spectrum, S (F) . 

Often the interferometer is treated as an inferior antenna which is 
ambiguous due to response of its many lobes to the object scene. Methods are 
then devised to resolve the ambiguities. The effort to resolve the ambiguities 
is made because the interferometer is much lighter and less costly than a 
filled aperture whose main beam equals the interferometer lobe width. Fig- 
ure 8-2 shows the variable baseline interferometer as a spatial-frequency 
spectrometer. 
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• INTERPORMETEfi IS SPECTRUM DISSECTOR OR 
SPATIAL-FREQUENCY SPECTRUM ANALYSER 

• IN CONTRAST, PILLED APERTURE IS IMAGE 
DISSECTOR: ANALYZES SCENE 1 PIXEL AT A TIME 
BY ANOLS SCANNING 

FIGURE 8-2. THE VARIABLE BASELINE INTERFEROMETER AS A 
SPATIAL-FREQUENCY SPECTROMETER 


Tt is obvious from previous discussion that, when operated In the * 

frequency domain, Che in to rf urometer is not ambiguous at all. In fact it is 
:| the pencil-beam antenna which xp ambiguous. 

T.ot the Pi 1 l<*d-n per Cure width he 0, ho chat Che pencil-beam and inter- 
ferometer lobe have the same width. Then the low-pass MTF of the filled- 
aperture has a passband which extends over the entire spectrum out to frequency 
D. All the components of the spectrum are passed resulting in complete 
ambiguity. 

The conclusion is that both the interferometer and the pencil-beam 
scanning antenna arc both useful, unambiguous mappers when operated in their 
proper domains. 


8.2 DISCUSSION Of* LATERAL DECORRELATION EFFECTS 

The above analysis of interferometer operation as a spatial-frequency 
spectrum analyzer may be confusing to those who have read radio astronomy 
textbooks and papers where the response of the interferometer to small uniform 
sources is analyzed. It is stated that when the source gets too big there is 
little response. A source subtending an angle of 1/S^ is too big. 

Kraus in his Figure 6-15 on page 175 shows why this is so by plotting 
the response of an adding interferometer to sources with three different 
angular extents. When the source becomes equal to or larger than the lobe 
spacing, the average of tKe source times the portion of the lobe pattern 
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falling on the source is nearly constant* There Is then no usable output 
signal. Tho bohnvior for n multiplying radiometer is similar, except that DC 
term Is not present. Kraus' Figure 6-14 is also given, since it shows how the 
point source response formed by multiplying the interferometer angular 
response by the antenna pattern of the interferometer antennas. 

• * From Figure 8-3 it may be seen that as the baseline length increases and 
the lobes become thinner and more closely packed, the response of the inter- 
ferometer to a finite source will decrease. For a uniform circular source of 
radius p located on the perpendicular bisector of the baseline at distance R, 
the normalized interferometer response is 


P 12 ‘V 


Hi « 

v 



(1-8-5) 


where i|» is the phase difference between radiation at the ends of the baseline 
nnd 


V “ 2ttS xr “ ” S x Aa (1-8-6) 

The angle subtended R by the source is 4a. • 



FIGURE 8*3.~RESPONSE OF AN INTERFEROMETER TO SOURCES OF VARIOUS SIZES 
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In Figure 8-4 the plot of due t0 Bom end Wolf (Section 10.4*3, 

pg 512, Figure 10.7) is shown. 

The abscissa is not V, but stems from an illustrative optical experiment. 
The first null occurs for v ■ 3.83. Let us assume that the interferometer 
output is decorrelated when ( P 12 1 Co 3 db below its peak value. Since 

|Ml2j is the amplitude of the correlator output, -3 db occurs when | V*12 1 “ 

1 /v2. At this point v ■ 1.6. , 

Substituting in tiq. (1-8-6) and solving for Aa, we get 


Au 


1 AO 

2SX " 2 


(1-8-7) 


Eq. (1-8-7) says that the disk cannot have an angular extent greater 
than half a, lobe width without suffering significant decorrelation. 


It may appear hard to reconcile these results with the earlier explana- 
tion of the interferometer performance in terms of a narrow bandpass filter 
which scans the spatial frequency spectrum of the scene. However these calcu- 
lations concerning the decorrelation of a disk source are in truth a help in 
understanding the earlier analysis. 


I*i*1 



FIGURE 84. NORMALIZED LATERAL 
COHERENCE FUNCTION OF A 
NONCOHERENT CIRCULAR DISK 
SOURCE 
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' The key to this puzzle is the fact that the terrain is not a featureless 
surface like the disk* It Is composed of structures of all sizes ranging 
from continental dimensions down to facets only a wavelength or so in size. 

The facets can be vegetation, wave facets, ice patterns, soil and rock patterns 
and many other things, including emissivity variations over the surface of 
man-made objects. Let us model this terrain structure as a surface over which 
a multitude of disks are strewn. The disks have random emissivities clustered 
around the average terrain radiant brightness. 

Let the disk diameters have a uniform probability distribution ranging 
from the operating wavelength up to the diameter of the mapped area. Disks 
larger or smaller than these limits cannot be seen. Call these limit 
diameters and d^x. Also let the interferometer be located at 

distance R from the surface. The angles subtended by the disks range from 
the angular field-of-view downward to X/R. That is 


t “ 

max 

FOV, 

Aa * 

max 

1 

rad 

i m 

min 

X/R, 

Aa . < 

min 

1 

rad 


Aperture synthesis starts with the travelling probe very near the reference 
probe. In order for the largest disks to contribute to the output of the 
interferometer correlator, Eq. (1-8-7) requires that 



2A« 


max. 


(T— 8-9) 


when S A ■ S A _ , all the disks contribute. As S A increases the contributions 
of the larger Sisks drop out. 


Let X ■ s x/S Amin ; then 1 <X<». Then the disk size corresponding to X is 

•»> • t>t)M • * 
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Assume that all disks smaller than ot(X) correlate completely ant I all disks 
larger than ct(X) do not correlate at all. The contribution to the correla- 
tion function of a disk will be proportional to <»2(X). Therefore 


aa oo m 2 

( d(X) - j a 2 (y)dy - b 2 J - -b 2 - | - ~ , 1SXS» (1-8-11) 

x x y 

We will set 0 (X) equal to zero in the region it cannot be measured, namely 
zero to one. Then from Vol. I, .Sec. 1.2, Eq. (3) of the Bateman Tables we 
find that the cosine Fourier transform of (8-11) is [Ci(y)]. Therefore 


f(x) 

«* 0 , 

0<X<1 





► 

f(x) 

1 

“ X" 

1<X<» 





-Ci(y) 


( 1 - 0 - 12 ) 


Figure 8-5 shows a plot of Cl(y) taken from the Bureau of Standards 
Tables (pg 232, Figure 5.6). This is the spectrum associated with the corre- 
lation function (lateral coherence ‘function) that was just derived. 

The spectrum does not extend below y ■ 1 because #(X) was not* measured 
for X less than ±1. This portion of the correlation function carried the 
information about the low frequency portion of the spectrum. • This gap in 
d(X) also produces the oscillatory behavior of the spectrum plot. The oscilla- 
tions would not be present in a more realistic model. Such a model would be 
more random than the one used here. 



FIGURE 8-5. PLOT OF THE COSINE 
INTEGRAL FUNCTION 
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» 

‘v However Che absolute values of the maximums and minimums of Ci(y) act 
as envelope which shows the behavior of the more random model. From pg 244 
of the Bureau of Standards Tables we find that 

! Ci(y oeak ) l " |CiC(n + l/2)ir]«-- 1 3 (1-8-13) 

. • peaK y peak y peak 


Remembering that fl(X) is represented as a voltage at the output of the 
correlator. It may be seen that the spatial frequency spectrum is of the form 


S(F) * F 2 (tf(X) } * | Cl (X peak ) | = i-j *• F" 2 (1-8-14) 

y peak 


This result is not too far different from the spectrum shape selected in 
Appendix TIT on the basis of other considerations. This spectrum was 
S(H) - K" 2 * 5 . 


8.3 DETERMINATION OF MAXIMUM BASETJNE T.ENOTlT 

| 

In most of the literature the maximum baseline length is treated as a 
specified parameter. A perfect mapping system is considered to be one having’ 
a MTF which is unity out to the critical frequency of the system. Bracewell 
defines the signal emerging from such a system as the principal solution. 

Since there is lots of room in orbit,' a less fatalistic attitude will be taken. 
It will be assumed that the baseline can be as long as we care to make it. 

As the baseline lengthens the spatial frequency spectrum of the scene will be 
measured with higher and higher resolution. This is illustrated by Figure 8-6 
due to Steel. (Reference 72) 

As spectrum resolution increases the rms difference between the signal 
representing the true scene and the signal produced by the mapping radiometer 
decreases. This rms error is one which does not decrease as the integration 
time available increases. In statistics an error of this type is called a 
bias error. The bias squared, b 2 , is the integral of this video power density 
of the scene from D to infinity. 

However as the maximum baseline length increases, the number of spectrum 
samples Increases. Since this maximum length is the radius of the disk over 
which samples must be taken, the required number of samples increases as the 
square of the maximum baseline length, D. The total observation time avail- 
able is fixed by the mapping geometry and the satellite velocity. Therefore 
the time available increases as D 2 as the maximum baseline length Increases. 
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FIGURE 8-6. PROGRESSIVE INCREASE IN 
SPECTRUM RESOLUTION WITH LENGTHENING 
OF THE CORRELATION RECORD 


2 2 

The total mean square error is the sum of AT and b added power-wise. 

It is obvious that this sum goes through a minimum for some value of D. Call 
this optimum baseline length D 0 t - This is a common type of statistical 
calculation. Figure 8-7 shows such- an error for a temporal filter due 

to Papoulis. (Reference 73) 


Oppenheim and Shafer (Reference 74) show plots of a spectrum where the 
total signal has been chopped up into records with lengths of 14, 51, 137 and 
452 samples. The total number of samples is 3,164. So there are about 226, 
62, 23 and 7 records, respectively, having the above lengths. A spectrum 



FIGURE 8-7. PLOT OF NOISE VARIANCE, AT 2 , 
BIAS SQUARED, b 2 , AND TOTAL MEAN SQUARE 
ERROR, c 2 FOR A TEMPORAL FILTER 
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estimate is made for each record; then all the records of a given length are 
averaged together. The results arc shown in Figure 8-8 (Oppenhcim and 
Shafer, Figure 11.3). 

Here the noise is the result of the fact that the variance of a spectrum 
estimate obtained by Fourier transforming a correlation function approaches 
a nbn-zero constant as N -* •>« The limiting variance of a spectrum estimate 
such as this was given in Eq. (1-5-8) of this report. 

It may be seen that N ■ 51 is near to being the optimum record length. 


In Appendix G, "Performance of Aperture-Synthesis Imagers," the optimum 
disk diameter as well as error expressions are derived. In Appendix XI, 
"Recapitulation of Radiometric Mapping Performance," plots of the rms noise 
error (square root of the variance) and the bias error are given for two 
aperture-synthesis systems. One system is in geosynchronous orbit. The 
other system is in low earth orbit. Plausible parameters are used in both 
cases. We have for GEO: AT - 30°K - rms terrain fluctuation, D 0 pt ■ 1.8 km> 
f ■ 1.2 GHz, ST ■ 3°K ■ temp res, b * 6.25°K « bias, e « 7°K “total error. 

We sec that for LEO: D opt - 1.5M, f - 20 GHz, <VT - 3°K, b - 8 K, c * 9°K. 
These temperature resolutions are poor compared to mapping radiometers being 
designed at the present time. This occurs because many designs ignore bias 
error completely, both bias due to imperfect spectrum reproduction of spatial 
frequencies less chan D and the finite-resolution error which is the present 
subject of discussion. • 
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FIGURE 8-8. PERIODOGRAMS FOR A FIXED TOTAL NUMBER OF SAMPLES AND 
RECORDS OF LENGTH N 
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' Probably the designer knows Intuitively that he cannot build an antenna 
large enough to reach optimum size, although this does not excuse failure to 
consider bias errors in the frequency range of 0 to D. This portion of the 
spectrum is the signal. Distortion of this portion of the spectrum is distor- 
tion oT the signal. 

* 'In many cases the system is used in applications which do not require 
resolution of detail smaller than some critical size. Then bias is found 
by integration from frequency D out to a frequency corresponding to this size, 
rather than out to infinity. This will decrease D 0 _ t and also the correspond- 
ing bias and noise errors. 

A second comment about this situation is that once bias error is 
acknowledged the system can be designed to reduce the error. Aperture 
synthesis does not turn in too good a performance record because bandwidth is 
limited by decorrelation effects. Figures 9 and 10 of Appendix show that 
truG-timo-dclay antennas do bettor. But the optimum 12M antenna in LEO is 
certainly at the limits of the state of the art; the optimum 3KM antenna in 
GEO is certainly impossible. However SPXNRAD (Figures 5 and 6) and RADSAR 
(Figures 7 and 8) can be built with optimum size. They also turn in fair 
radiometric performance (iST ~ 1.5 « 2°tC and bw 2.5 - 3°K) . Contemporary 
design trends in radiometry such as Hadamard transform techniques (Refer- 
ence 75) will probably reduce the errors values just quoted considerably in 
the future. 


8. A BIAS DUE TO FINITE RESOLUTION OF THE SPATIAL FREQUENCY SPECTRUM 

Up to now we have only been concerned with bias due to finite angular 
resolution. There are other bias errors which are present and should be 
mentioned. One bias is due to the finite bandwidth of the bandpass MTF of 
the interferometer. In Figure 8-2 it may be seen that as the MTF scans over 
S(F) the result is the convolution of the two functions. The multiplicative 
relation, S(F) ■ d(F) * MTF(F) holds for a fixed baseline ’length. Let S(F) 
be the estimate of S(F) generated by the scanning MTF of the interferometer as 
the baseline length is changed. Then 


3(F) - d(F) * MTF(S A ) (1-8-15) 


This result is completely analagous to the situation encountered when a 
temporal spectrum is measured by a scanning temporal filter. The bias for 
this case is found by Bendat and Piersol (Reference 76). It is 


B 

b [s(F)3 - S"(F) (1-8-16) 

i 

where B is the filter bandwidth, 
e 
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The Courier transform of the MTF passband shape is the element pattern 
of the two probe antennas alone. The interferometer pattern appears only 
When the pattern origin is made the center of the baseline and the probe 
patterns have phases referred to that origin. The power pattern for the 
interferometer is then generated as shown in Figure 8-4. The element pattern 
is given by 


E(X) « G(X) • G*(X) 


(1-8-17) 


The relations (8-2) for a fixed baseline are therefore replaced by the 
relations in (8-18). §(X) is the estimate of the intensity distribution 

function for the scene. % 



P Q (X) • E(X) 

ir jr 

tf(F)*MTF(S.) 


(1-8-18) 


Note that $(X) is quite n good reproduction off P () (X), if the probe 
.'incciiitmi Iwivu m» aper ture size wliicli is small eompured to the muxi|iium buse- 
line length, D. Then MTF(S^) is nearly an impulse. Also E(X) is nearly 
uniform across the FOV subtended by the scene, except near the edges. 


It would appear that P 0 (X) could be recovered by dividing P(X) by E(X). 
However this inverse filtering operation will amplify noise greatly near the 
edges of the scene where E(X) is small. The solution to this difficulty is 
use of a Wiener filter, rather than an inverse filter. The Wiener filter 
weights the degree of inversion carried out by the signal- to-noise ratio. 
This will be discussed in more detail in a. later section. 


ft might be mentioned that tiie response of the aperture synthesis system 
over the spatial frequency band from zero to D is almost rectangular. This is 
much closer to perfect rectangular response than the response of a filled 
aperture of size D (shown in Figure 6-1) . As a 'result the aperture synthesis 
system has almost twice the resolution of a filled aperture having the same 
size. 


8.5 ROTATION SYNTHESIS 

* 

As mentioned in Section VII, "Frequency-Domain Imagers," the lateral 
coherence function is measured over a disk in the Fourier plane by rotating 
the interferometer baseline while repeating the aperture synthesis procedure. 
These repetitions are made often enough to meet the Nyquist criterion. The 
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Van Cittert-Zernike theory says that the two-dimensional Fourier transform of 
the lateral coherence function is the source brightness distribution 
(Appendix E). 

There is a second relation between 1^2 and B which is less obvious. Let 
r 12 (t) be the lateral coherence function measured over displacement t. Then 
the one-dimensional transform of (t) is the brightness distribution of -i 
the scene measured alon^ a cut through S(u, v) which has the same orientation 
as the displacement x in the frequency plane. This is the projection theorem. 

, This theorem shows that B(x, y) can be built up out of a series of 1-d 
transforms of Tx 2 measured along a set of diameters of the Fourier disk. 

No matter how B is computed from the process of measuring r^ 2 along 

a set of disk diameters in the frequency plane followed by computation of 
B(x, y) is called rotation synthesis. 


8.6 ALTERNATE WAYS OF SCANNING THE INTERFEROMETER MTF IN SPATIAL FREQUENCY 

Up to now wo have examined aperture synthesis where the length of the 
interferometer baseline has been varied. There are two other means of varying 
the effective baseline length. 

The first alternate method is variation of the operating frequency. This 
changes the variable S\ ■ S/\ by varying X instead of S. Sheridan (Refer- 
ence 78) built a system of this sort for mapping the sun (Reference 78). It 
is apparent that the source must have the same emissivity over a very wide 
range of frequencies when this technique is used. Also', the radiometer must 
be located where it will not receive radiation from man-made source. This is 
a very hard condition to meet. 

The MFASR avoids this problem by using a discrete set of frequencies '<U 
which fall at places in the radio spectrum where there is little interference. 
The question then arises as to whether there are enough samples available to 
permit a reasonable reconstruction of the image. This problem is treated in 
later sections. 

I 

The other possible method of varying S,\ is rotation of. the inter- 
ferometer end-over-end. This is illustrated in Figure 8-9. 

The axis of rotation is at right angles to the axis oi' rotation for 
rotation synthesis. If the field of view is centered on the zenith-nadir 
line ZN, the effective interferometer baseline length equals D cos 8. This- 
varies from a value of D when the baseline is horizontal to zero when it falls 
on ZN. 

The field-of-view must be held fixed by rotating the radiators on the ' 
baseline ends in the opposite direction from the base rotation. This is 
indicated in the drawing by the circular arrows marked -0‘. 
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FIGURE 8-9. ROTATING INTERFEROMETER BASELINE 
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9. IMAGE FORMATTING AND OPTIMIZATION 


9.1 INTRODUCTION 

Frequency domain Imagers rend to invite a much larger amount of data 
processing than conventional imaging systems. In the case of conventional 
systems, the image in formed by hardware. Image-data processing is then an 
additional step. If the system designer elects not to take this step, the 
system still delivers images to the user. 

In the case of frequency-domain images, data processing is mandatory. 

Even if nothing else Is done, the spectrum data must be Fourier transformed 
in order to produce the image. Once data is digitized arid fed to a computer 
iu order to effect the frequency-to-spacc-doniain transition there is strong 
motivation for additional data processing in order to improve image quality. 

Image manipulation by computer tends to quantify image quality to a 
degree which never occurs when a lens or pencil-beam scanner forms the image. 
This quantification of image properties leads to definition of the operations 
on the image data which are needed to produce the best possible images. This 
in turn leads to specification of the performance of the rest of the imaging 
system required to get the best results out of the image-data processing. 

Although there are grey areas, image processing can be broken down into 
three areas. One area I will call formatting. This involves arrangement of 
the image data for best display. Examples are image annotation, dynamic range 
adjustment to match the capabilities of the display, contour line generation, 
and map-projection generation. 

A second area is image manipulation in order to improve quality according 
to some criterion. 

A third area is measurement of image parameters and comparison between 
images in order to quantify image quality. 

Fortunately, there is a software section VICAR which meets most routine 
image processing needs. The rest of this chapter and Volume II are concerned 
with image processing associated with specific Fourier-plcne imager systems. 
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9.2 THE SOFTWARE SYSTEM - VICAR 

This software is described by Castleman (Reference 117) in his Chapter 4 
and Appendix K, VICAR PROGRAM INDEX. Appendix M includes a library of FORTRAN 
language imago-processing programs. 

* ‘The VICAR system, and its documentation, may be obtained from the 
following. • 

COSMIC 

Computer Center 

112 Barrow Hall 

University of Georgia 

Athens , GA 30601 

(404) 542-3265 

The VICAR system is filed under number NPO-13415. A set of image- 
processing programs for the PDP-11 computer are also filed there under NP0- 
14892 and the name Mini- VI CAR. 

Appendix K. is a reproduction of Castleman* s Appendix II where a list of 
VICAR programs are given together with brief descriptions. Asterisks indicate 
that the program was available from COSMIC at the time Castleman* a book was 
published. 


9.3 IMAGE QUALITY CRITERIONS 

The most obvious criterion for image quality is the root-mean-square 
difference between the final image and the original scene. The total rms dif- 
(Vrititcc coutiincti of two components. One component is system noise. Its mag- 
nitude is given by AT, the system temperature resolution. The second component 
is called bias. 

In statistics bias 1c defined as error which does not go to zero as the 
averaging time increases without limit. One source of bias, the finite resolu- 
tion of the system, has already been discussed in Section 8.4. There a method 
of minimizing the sum of the system noise variance and the square of the finite 
resolution bias was presented. As the linear array length or maximum inter- 
ferometer baseline length increases system noise increases. This occurs 
because the number of spectrum samples required increases. However, the resolu' 
tion bias decreases as the system resolution Improves. It was found that there 
was an optimum array length where the sum of the noise variance and bias 
squared went through a minimum value.* 

Radiometer mapping performance for this optimum configuration was com- 
puted for three types of mapping systems in appendices to this report. These 
three mapping types are aperture synthesis, pencil-beam scanners using series 
fed antennas and true-time-delay antennas. 
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“ . There are additional bias errors due to causes other than finite resolution. 
They will be discussed' later in this chapter, A common type Of bias error is 
nonuniform response of the system MTF between zero and the cutoff spatial fre- 
quency. That is the system does not produce the principal solution. In 
Bracevell's terminology the principal solution is obtained when the mapping 
system exhibits a rectangular MTF from zero out to the cutoff frequency. 

» » 

Bias error due to motion of the scene is an important error type. It will 
be found that such motion produces a nonrectangular MTF which must be multi- 
plied into the chain of transfer functions. This idea will be developed in 
more detail in a little while. Corrections for motion effects is therefore 
carried out as part of the process where deviation of the system MTF from 
rectangular is corrected. 

The optimum filter for carrying out this correction is the Wiener filter. 

It will be discussed in the next section. If the image is to be used for visual 
inspection, rather than measurement, it will be found that the Wiener filter 
emphasizes low frequencies too much. The Wiener restored image appears dull 
and flat, rather than crisp. Some of the finer detail appears blurred. 

In Section 5.8, the ideu was introduced that the optimum image for visual 
inspection maximizes the area under the system MTF. The MTF of the eye is 
included in the chain of MTF's whose product is the system MTF. A filter which 
produces this result is the parametric Wiener filter (geometric mean filter). 

Tli I a I’ M Lor . /rtf wu|, J. mi other filters for improved image uppeurnnee, are dis- 
cussed by Hunt (Reference 118) and by' Stockbara (Reference 119). Both papers 
give examples, of the improvement of image quality compared to Wiener-filtered 
images. 

, A final Criterion of image quality is the entropy of the image. Informa- 
tion theory shows that the entropy of a signal record is a measure of the infor- 
mation it contains. Therefore, if the conversion of the image signal from its 
spectrum .form to the final space-domain image can be carried out so as to 
maximize image entropy, it can be argued that the best possible image results. 

Presumably, this "best" image can now be processed so as to meet any other 
quality criterion since it contains more data about the original scene than any 
other possible image. 

t i 

The Importance of these observations lies in the fact that algorithms 
exist which csm. maximize image entropy in the process of reconstructing the 
image from spectral data. These algorithms go under the appellation, "maxi- 
mum entropy method," or MEM. They are very important when the image must be 
reconstructed from under sampled spectrum data. 

Since this is the problem we face in mechanizing the MFASMR system, these 
algorithms are treated in detail in Volume II of this report. Image reconstruc- 
tions using MEM are presented there. Only a brief introduction to these 
methods will be given in this chapter. 
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The. last type of optimum filter to be considered Is the matched filter. 
This filter exhibits the maximum possible response when driven by the Input 
waveform it was designed to match. The filter response does not resemble the 
input waveform In any way. Instead, the output is a pulse which approaches a 
delta function as the time-bandwidth product of the waveform it matches 
approaches infinity. This filter maximizes the probability of detecting the 
matching waveform in the presence of Interfering noise. 


9.4 WIENER FILTERS 

If the modulation transfer function of the radiometer system is not 
rectangular, it would nppeor that nil that had to be done to make it rectan- 
gular would be to add a filter to the MTF chain whose MTF is the reciprocal 
of the MTF of the rest of the system. This filter is called an Inverse filter. 

Unfortunately, in regions where the system MTF is zero the Inverse filter 
transfer function becomes infinite. Figure 9-1 shows that even when the system 
MTF has no internal zero regions, noise becomes dominant at the high frequencies 
where the system is tending toward zero. The result is a noisy image. 
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FIGURE 9-1. NOISE PROBLEM IN DECONVOLUTION 
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Like the Wiener filter, the PSE filter is phaseless (real and even). Both 
filters reduce to straight deconvolution in the absence of noise, and both 
cutoff completely in the absence of signal. 

However, the PSE filter (sometimes called a homomorphic filter) does not 
cutoff at internal zeroes of the point-spread function F(u,v). 

The other filter related to the Wiener filter is the parametric Wiener 
filter. It is of the form 


°PW (U > V) 


F*(u,v)l ? 

r F*(u,v) I 

Jf(u,v)| z _ 

_ |F(u,v) | 2 Pg(u,v) + P N (u,v) J 


(1-9-3) 


This is a generalization of the filters already discussed. When a * 1 
Eq. (1-9-3) becomes the inverse filter. If a * 0 and y ■ 1 we have the Wiener 
filter. When Y is not unity the result is a parametric Wiener filter. The 
constant y may be selected for any desired amount of Wiener-type smoothing. 
When y ■ 0 we are back to the straight deconvolution filter. If o = 1/2 and » 
y ■ 1, Eq. (1-9-3) defines a filter which is the geometric mean between ordi- 
nary deconvolution and Wiener deconcolution. The parametric Wiener filter 
with y < 1 and the geometric mean filter with the same constraint define the 
range -for y where the best visual images are obtained. 

* 

In the construction of error budgets needed to judge the suitability of 
a system for a given application we need to be able to compute the residual 
rms error after Wiener filtering. Both signal and noise power density spec- 
trums enter into the determination of this error. The error expression is 


2 

e 



“p s (u,v)P n (u,v) 
Pg(u,v) + P N (u,v) 


dudv 


(1-9-4) 


This is a generalization of two dimensions of an expression given by Bendat 
(Reference 120) for th £ rms error of a Wiener filter. 


Now Eq. (3) , Appendix G gives the signal spectrum in polar coordinates 


as 
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where <ST is the rms radiance fluctuation of the object scene, FOV » 
field-of-viaw < 1 rad. 


Tor a uniformly illuminated antenna the normalized transfer function is 
(1-p), 0<n<l. Lot the noino spectrum ba flat with typical rms amplitude of 
0.5 to JokT Since <$T varies from 20°K for rural areas to A0°K or so for urban 
areas, ice/water or land/water images, (AT/5T) 2 is much less than one. So in 
polar coordinates we have 
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f V 2lT p s (p,6)(i-p)At 2 p 
IJq Fg(p ,0) (1-p) + A*/ 


1-Ap 
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4rrvT5vp 5/2 /ci-P> [(AT/6T) 2 ] 


(1-9-6) 


As long aft p is n l ittle loss chan 1 the second term of the denominator 
is small. Therefore, 


a" 2 

- 2irAT 4 T* 


1-Ap 


trAT 


(1-9-7) 


Therefore, the total rms error is approximately AT^. 

It is important to note that Wiener filtering cannot be used directly for 
image reconstruction when the spectrum is undersampled. This difficulty arises 
because missing data samples are represented by zeroes. The Wiener filter will 
cut off in any spectral region with zero spectral power density. An estimate 
of the original spectrum with undersampling eliminated must be generated by 
interpolating in the regions where samples are missing. Some sort of two- 
dimensional spline function is needed. 

If polynominal splines are used, the polynominal coefficients roust be 
determined so that the spline gives known values Of the spectrum where samples 
exist. As an alternate wo can use an expansion into a set of orthonormal func- 
tions. As in tho polynominal spline cose, the series coefficients are picked 
to make the expansion give the correct values at points where samples exist. 

If Wiener filtering is carried out after the spectrum has been filled in 
by interpolation, it is not obvious that we still have a system which minimizes 
rms error. The performance of Wiener filtering when they input data is under- 
sampled is treated in more detail in Volume lit where computer results are 
reported. 


1-9-7 

\ 


ORIGINAL PAGE IS 
OF POOR QUALITY 



i 



9.5 MAXIMUM ENTROPY RECONSTRUCTION METHOD 


It will be found in the next chapter that the amount of distortion 
calculated for the MFASMR is probably too great to make the system very useful. 
These calculations are based on a discipline called "rate distortion theory." 
This theory assumes that the image reconstruction process is linear. However , 
there are more general non-linear methods which are more powerful than linear 
image construction. 

As mentioned in the introduction to this chapter, one of the most promi- 
nent processes of this type is the maximum entropy method, MEM. This process 
is recursive. Rather than requiring a priori statistical data about signal 
and noise, MEM effectively collects the required statistical* data during the 
recursion cycles. It goes without saying that MEM is nonlinear. 

In order to determine whether MEM or similar techniques can reduce dis- 
tortion in MFASFMR to a tolerable level, computer experimentation was used. 
This work is reported in detail in Volume II. Only a few comments will be 
made about MEM at this point in order to paint a complete picture of the image 
reconstruction problems - albeit a sketchy one. 


An excellent place to start a survey of MEM techniques is the IEEE 
Reprint Volume edited by Childers (Reference 121) . The promise of MEM is well 
illustrated !>y a compar (noil of spectrum estimates derived from the. spectrum 
autocorrelation function made by Abies (Reference 122). Abies' plots in Fig- 
ure 1 to 4 are reproduced here as Figure 9-3. It may be seen that the MEM 
spectrum estimate is considerably better than either the weighted or the 
unweighted Fourier transform of the autocorrelation function of the signal. 





FREQUENCY 


FIGURE 9-3. COMPARISON OF FOURIER — TRANSFORM 
SPECTRUM ESTIMATE WITH MAXIMUM ENTROPY 
SPECTRUM ESTIMATE 


1-9-8 



ORIGINAL PAGE IS 
OF POOR QUALITY 


In spite of fine results obtained with HEM, the method has generated 
considerable controversy. Some MEM algorithms have been unstable. As the 
number of recursions increase some MEM algorithms produce artifacts. Typically, 
these artifacts are sharp spikes which grow on the edges of spectrum lines as 
recursion continues. An example of this behavior, taken from a report by 
King (Reference 123), is shown in Figure 9-4. The fourth recursion is the best 
estimate of the original spectrum. The spikes, which are well developed by 
the seventh recursion, obviously increase the entropy of the signal. Apparently, 
the MEM estimate increases entropy beyond that of the true signal by producing 
artifacts. 

However, Wernecke (References 124, 125, 126, 127) has developed a stable 
MEM algorithm for image reconstruction. Stability is ensured by imposing a 
constraint on the calculation of the estimated image during the recursion pro- - 
cess. This constraint prevents any point in the image estimate, being formed 
during a recursion cycle, from taking on a negative value. This condition is 
of course consistent with the nature of image signals, since such signals are 
always positive everywhere. 

This condition may not appear to he a very severe constraint. However, 
it is sufficient to ensure stability of the MEM reconstruction process. Almost 




ANGLE, DEGREES 

FIGURE 94.' GROWTH OF~ARTIFACTS DURING RECURSION BY 
AN UNSTABLE MEM ALGORITHM 
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*11 possible unstable processes involve oscillatory behavior which grows with 
time. This oscillatory disturbance drives part of the image estimate negative. 

I cons training the image estimate to be always positive damps out this oscilla- 
tory behavior. Then recursion cycles of all orders are stable. As recursion 
continues the entropy of the images estimates asymptotically approaches a maxi- 
mum value. 

Volume II will shows results for real radiometer image-data. However, a 
"sneak preview" will be given here where the original object scene is two cubes 
stacked on top of each other. This is similar to a test object used by Wernecke 
in Reference 124. Wernecke 's "stacked blocks" object is shown in his Figure 3.7. 
In Figure 9-5, which follows, the stacked blocks phantom used by Dr. Ming Chang 
Jn shown. I)r. Cluing carried out the* computer study of image reconstruction for 
MFASMK reported in Volume II. The following data comes from this study. 

Following Figure 9-5, Figure 9-6 shows the amplitude spectrum of the 

I stuckud blocks object. Figure 9-7 shows the residual spectrum after samples 
are taken over the spectrum plane as shown in the upper right hand corner of 
the figure. We now have an undersampled spectrum. 

Figure 9-8 shows the reconstructed image obtained by taking the two- 
dimensional inverse Fourier transform. As one can see, the resulting image is 

I quite poor. 

Figures 9-9, 9-10, and 9-11 show the MEM image reconstruction process 
taking place. The results of 2, 4, 6, 10, 18, 22, and 25 iterations are shown. 
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Although the Image is far from perfect, after six to ten Iterations, it is 
considerably better than the direct inverse Fourier transform. The thin ridges 
extending out from the cube corners are similar to the impulses which appear 
at function discontinuities when the Fourier transform is taken without under- 
sampling. This behavior is called Gibbs' phenomenon. Since this Gibbs arti- 
fact occurs in both reconstruction methods, it does not appear reasonable to 
fault the MEM algorithm for this generating this imperfection. 

* f 

Continued iteration up to 25 recursions shows that the reconstruction is ’ 
stable. It does not have spikes which grow indefinitely as was the case for 
King's MEM algorithm (Figure 9-4). 


9.6 IMPACT OF MEM RECONSTRUCTION METHODS ON SYSTEM DESIGN 

Even when there is no undersampling MEM reconstruction is an attractive 
method of forming the image. The MEM method is adaptive in that we do not 
have to know a priori statistical data about signal and noise. The signal 
statistics directly controls the MEM recursion as the image is built up. 

As we have seen, Wiener filters require quite a bit of signal and noise 
statistics. This is a severe drawback, since the mapped terrain hardly quali- 
fies as a stationary source. The statistics vary as we look at farmland, for- 
rest, the sea, ice fields, etc. There are extensions of Wiener filter which 
can adapt the filter to a locally stationary signal (Reference 128 ) ( . However, 
these systems are complicated. MEM reconstruction appears to be tiie better 
answer. 

Now it stands to reason that if' we are going to maximize the entropy of 
the final image by use of MEM reconstruction, we should design the rest of the 
system to produce a signal with maximum entropy rate at the input to the MEM 
process. Appendix G,. "Optimizing Antenna Illumination and Bandwidth for Maxi- 
mum Radiometric Information Transfer," by Dr. Jack Gustinsic shows how an entire 
radiometric mapping system can be designed to maximize signal entropy. 

It was hoped that mapping systems, other than the series-fed array system 
selected by Gustincic, could be treated in the same manner. However, there is 
not time enough to do this at present. Also the size of this report, which is 
probably too big already, would increase even more. However, it is possible to. 
comment on some general design principles. 

First Gustincic found that maximum information rate occurred for approxi- 
mately uniform antenna weighting. This was for a system using a series fed 
array. In the other mapping systems considered, which use true-time delay 
antennas or transform arrays, it turns out that exact uniform weighting is 
optimum. Similarly, the spatial frequency spectrum measured by aperture syn- 
thesis should be uniformly weighted in order to obtain maximum information rate. 
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The overall ‘transfer function of the radiometer receiver, treated as a 
filter, is subject to a similar constraint. Goldman (Reference 129) shows that 
the entropy loss through a filter is zero only when the filter has a rectangular 
paseband. The virtues of rectangular passbands for radiometers has been recog-" 
nized for quite a while, even though improvement has been given as a decrease in 
temperature resolution rather than increase in Information rate (Reference 130). 


9.7 MOTION COMPENSATION 

* 

If the image data Is not under sampled, there is one type of error which 
should be removed by Wiener filtering before MEM image reconstruction* This 
error is due to image motion. Removal of blur due to motion is very important 
when the mapper is in low-earth orbit. 

Let us examine the signal at an image plane, such as a focal plane of a 
lens. This image can be decomposed into an assembly of points. Let these 
points move at a constant velocity v due to motion smear. Let the time 
required to sense the image be t. This can be the exposure time of a camera, 
the time required to measure the frequency spectrum of the scene or the time 
needed to complete a raster scan with a pencil beam. The displacement of the 
point during the observation time is AX ■ mvt. The smeared point of bright- 
ness I is described by the expression 


f (x) 


A, 


_£2L < x < .AX 

2 - X -2. 


1 ^ 0 , otherwise 


( 1 - 9 - 8 ) 


where A is proportional to I. Thus each point is converted into a rectangle. 
The corresponding Fourier transform is 


C MOTION 


sin (FAX/2) 
*** (FAX/2) 


(1-9-9) 


We can conclude from Eq. (1-9-9) that image motion introduces a sine transfer 
function into the system chain of MTF's. It may be seen that, in conformance 
with earlier comments, the motion error cannot be removed by introduction of 
the inverse filter A/G(F) . Its poles at F ■ 2mr/Ax render it unstable. The 
"correct” image will be dominated by noise. 

The removal of image blur by digital and electronic filtering has been 
treated by Sawchuk (Reference 131, 132, 133). Optical filters for removal of 
motion blur by optical processing is treated by YU (Reference 134) and by 
Tsujiuchi (Reference 135). In Figure 9-12 the rectangular response to the 
moving point is shown in a plot from Tsujiuchi' s paper. The restoration of 
the point image is accomplished by use of a dual filter. One filter section 
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FIGURE 9-12. CALCULATED IRRADIANCE OF A 
RESTORED POINT IMAGE BLURRED BY LINEAR 
* MOTION 


modifies only amplitude; the other section modifies only phase. The plot 
shows the effect of each section separately and their combined effect. The 
smeared point is labeled original. The system impulse response to a point 
source is labeled "Combination Amplitude and Phase Filter." 

In the following discussion it is important to note that the error was 
not produced in the sensor. It was the image data itself which is degraded. 
Let us examine the response of pencil-beam scanning imager to this degraded 
data. 


Jacobs (Reference 136) -computes the degradation of the spatial-frequency 
response of a scanning pencil-beam imager due* to image motion. He uses the 
product of the pencil-beam scan rate . and the radiometer time constant as 
a parameter which characterizes the rate of image motion. The scan rate is 
tied to image motion by .the requirement that the scan process must be rapid 

enough to keep up with the image motion. 

* 

The degraded MTF due to scanning is found to be 


MTF'(iu) 


•~=M=== 

>/l + (tvw)^ 


(1-9-10) 
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where MTF(w) Is the transfer function of the antenna and MIF' (<o) is the 
transfer function degraded by scanning. The degradation factor [l + 
was plotted by Jacobs and shown below as Figure 9-13. 

We find that pencil-beam scanners also suffer degradation due to Image 
motion. For some reason, degradation due to image motion is usually not cal- 
culated when the system performance of such mappers is estimated. Also pencil- 
beam scanner systems almost never use restoration filters to compensate for 
motion. However, these systems are using the same Motion-degraded images that 
frequency-plane imagers use. Degradation due to the motion is just as' severe. 
Motion correction will be just as beneficial. 

Consideration of these effects when evaluating both imager system types 
will lead to better comparisons of their performances. 


9.8 BIAS DUE TO SPECTRUM CURVATURE 

In Figure 8-2 the bandpass MTF of an interferometer is shown measuring 
the spatial frequency spectrum of the scene. In the next chapter Figure 9-4 
shows how measurement of the spectrum of the signal received by a frequency- 
scanning array produces an image of the terrain. 

In both eases these, spec trams lire measured by f'! 3 tors with Finite pnss- 
batidu. A Filter producing no error would give the power density of tike spectrum 
at a frequency equal to the center frequency of the filter passband. 



FIGURE 9*13. IMAGE DEGRADATION FACTOR ASA 
FUNCTION OF NORMALIZED SCAN RATE 
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Due to spectrum curvature and the finite bandwidth of the spectrum 
measuring filters they do not find the midband spectrum density. The measure- 
ment is perturbed by a bias error. According to Bendat and Pier sol (Ref- 
erence 137) this bias error is given by 


b 2 t£(f)J 


If 0"<£) 


(Ir9-ll) 


where G M (f ) is the second derivative of the spectrum density function and B is 
the bandwidth of the measurement filter. There is no point in computing these 
biases for the raw spectrums. Their shape and curvature will be altered by 
filters used to correct for image errors in the image restoration process. 

However, this bios should be computed along with the other system biases 
in the process of drawing up a bias budget for the purpose of estimating sys- 
tem performance. 
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10. UNDERSAMPLING AND IMAGE RECONSTRUCTION 
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10.1 INTRODUCTION 

At this point, enough background material has been introduced to proceed 
with an analytic evaluation o£ the Multifrequency Aperture Synthesizer Microwave 
Radiometer (MFASMR) . This frequency-imaging system varies the baseline length 
of its interferometers by stepping the operating frequency through a set of 
values scattered over most of the microwave spectrum. The frequencies selected 
are: 2.375, 4.75, 9.5, 19.0, and 38.0 GHz. 

A typical MFASMR configuration is two interferometer baselines at right 
angles to each other. Baseline lengths of 10M with 1M dishes is also repre- 
sentative. Such a system is shown in Figure 10-1. 

In Section 7.5, a paper by Ko, (Reference 68) "Coherence Theory in Radio- 
Astronomical Measurements," was referenced. In this paper, Ko shows that the 
maximum baseline length in wavelengths is the reciprocal of the system resolu- < 
tion in radians. Ko also states that, in order to prevent aliasing, the sample 
spacing in wavelengths is the reciprocal of the FOV in radians. Section 7.5 
point's out that this implies the moving interferometer antenna is stepped 
along the baseline in steps equal to the antenna width. When this observation 



FIGURE 10*1*. TYPICAL MFASMR 
CONFIGURATION 
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is applied to the configuration in Figure 10-1 , we see that it generates a 
10-line image (10 x 10 pixels) , since the number oE pixels in the final image 
equals the number of complex spectrum samples. 

Since five complex samples along the two baselines can be taken, using 
five operating frequencies, we see that 10 samples out of a grid of 100 samples 
can be taken. Only 10 percent of the raw image data is retained. The problems 
associated with undersampling the spectrum of the scene dominate the processing 
of the MFASMR data. As a result, it was decided to incorporate the analysis of 
its performance in this chapter on undersampling and image reconstruction when 
undersampling is present. Volume II describes a computer study of the same 
problem. This study is more detailed than the analytic attack used here. 


10.2 COMMENTS ON THE MFASMR ANTENNA CONFIGURATION 

Before proceeding to MFASMR performance analysis, it is necessary to take 
a look at the MFASMR antenna arrangement. This system is, of course, a fre- 
quency-plane imager. The antenna arrangement, plus the designation of the 
operating frequencies, determines the sot of samples of the scene spectrum 
available for image reconstruction. 

The first comment concerns the use of two baselines at right angles to 
each other. As mentioned in the Hughes proposal, thi6 arrangement conveys 
a strong feeling that the spatial frequency spectrum being sampled .is a 
separable function. A separable function is the product of two functions, each 
of which is a function of a single independent variable. Let F(U,V) be the 
spectrum function, if it is separable, it is of the form R(U)T(V), where R(U) 
and T(V) are known over two boundaries of the region in which F(U,V) exists. 

Wa know from theoretical physics that in order for separation to occur, the 
set of eigenfunctions of R and T must have nodes over their corresponding 
boundary surfaces. This rarely happens and then only for a restricted set of 
coordinate systems. In such cases, specification of R and T over their 
corresponding boundaries forces these functions to assume certain values 
throughout the region occupied by the field, namely, values which are the sum 
of the eigenfunction set with appropriate weighting. 

It is easy to see that F(U,V) fails this test. Specification of F(U,V) 
along the interferometer baselines does not constrain F(U,V) over the rest of 
the. frequency plane in any way. 

A more formal proof that F(U,V) is not separable stems from the projec- 
tion theorem. This theorem showed that there is a one-to-one correspondence 
between profiles generated by cuts through the space-domain function, say 
g(8,R), and corresponding frequency-domain profiles obtained by a cut through 
F(U,V), which has the same orientation as the cut through the object function, 
f(x,y). These two profiles are a one-dimensional Fourier-transform pair. The 
above relation between F(U,V) end f(x,y) implies that we must measure F(U,V) 
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profiles over a set of orientations distributed over 180° with Nyquist sampling 
density in order to be able to calculate f(x,y) exactly. These measurements 
of the spatial frequency over a set of orientations is the process of rotation 
synthesis, which by now has become familiar to us. 

The minimum number of orientations which can be used with any hope of 
image reconstruction are cuts in the U and V directions and along orientations 
at ±45° with respect to the principal axes. This is the sort of sample set 
used by Werneke (Reference 124) shown in Figure 3.8, page 96 of his doctoral 
dissertation. The same sampling scheme has been used by Dr. Chang when 
checking out Werneke 's MEM program. It is shown in the upper right-hand cor- 
ner of Figure 9-7 in this report. 

In order to obtain such a sample set, another antenna must be added to 
the original MFASMR configuration as shown in Figure 10-2. 

Thu diagonal baselines Aii and AD contribute absolutely essential additional 
samples. Mow the samples, although sparse, are indeed distributed over the 
entire disk in the frequency plane corresponding to the scene. The potential 
baselines DC and DB are not used. Since translation of a baseline without 
orientation change does not affect the value of the mutual coherence function, 
these baselines merely repeat the measurements made using baselines AC and AB. 

We have added ten more samples to the ten available with the original con- 
figuration. So we now have 20 samples out of a possible 100. The undersampled 
spectrum is 20 percent of the total spectrum. This is also a statement that 
the information rate of the undersampled data stream is 20 percent of the infor- 
mation rate of the sample stream associated with the Nyquist sampled spectrum. 



FIGURE 10-2. MODIFIED MFASMR CONFIGURATION 
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10.3 RATE DISTORTION THEORY 

Once the information rate of the coded image signal is known, we can ask 
what error is produced in the reconstructed image when the information rate is 
less than the information rate for Nyquist sampling. Theoretical techniques 
have been devised which find an upper limit on this error. These techniques 
go under the name of "rate distortion theory." (Reference 140) 

The mean-square error introduced into the signal as a function of the rate 
of source coding and information transmission through a communications' channel 
is best expressed in parameteric form. For this mean-square-error criterion 
and for a Gaussian source this parametric representation is analytically 
tractable. (References 141, 142) 

The analysis assumes optimum coding. If the coding is optimum, then the 
distortion is an upper bound on the error for all sources with the same spec- 
trum, but with other than Gaussian statistics. 

The parametric form for the rate of source coding and for information 
transmission used the parameter 0. In Figure 10-3 it can be seen that 0 is 
a decision level in a -sorting process . 

As indicated by the hatcher code, the portion of the power-density spec- 
trum above the level 0 is retained and transmitted. The crosshatchcd area, 
which can be interpreted us the frequency distribution of the error process, 
x(t) -y(t), is often called the error spectral density. In ideal mean-square- 
error (MSE) encoding, no attempt is made to preserve any of the information con- 
tained in those regions of the frequency axis in which 4»(aj) <_ 0, where 
#(o>) > 0, a fraction [*(w) - 0]/4(m)o*f the source power generated in an infin- 
itesimally narrow band centered at u is retained in the ideal reproduction. 

In other words , one should try to reproduce the portion of the source output 
generated at frequencies such that *(w):> 0 wi v ,h a MSE of 0 per unit bandwidth. 
For 0 somewhat higher than illustrated in the figure (hence with distortion 


PRESERVED SPECTRAL DENSITY “ 
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FIGURE 10-3. SPECTRAL DENSITIES OF SOURCE, ERROR, 
AND REPRODUCTION FOR OPTIMUM CODING 
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somewhat larger) , the portion of $(«) centered at to * 0 is the only one that 
an ideal system attempts to reproduce, for such 0 one may put the process 
through an ideal low-pass filter prior to encoding. In general, a filter with 
several passbands and stopbands is required. 

Then in terms of this parameter 6, from Berger, Eqs (4.5.51) and (4.5.52) 
the mean-square-error and information rate are 


d(0) - 


R(e) - 
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The variance of the process, which is used to normalize d(0) , is 
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If the spectrum is one-sided, the integral limits of ±® are replaced by 
the limits 0 and -H». 

In the two dimensional case the differencial frequency element in the 
Fourier plane is of the form (p d<J>dp). Starting with the Cartesian form of 
the two-dimensional integrals from Reference 142, Eqs (8.4.45) and (8.4.46), 
we see that in polar form all the integrals are of the form 
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Before looking at our case let us examine some typical results from the 
literature. As an example, Berger uses a Gaussian source with spectral density 
of the form 


(4.5.76) «(u) 
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The equation number to the left is Berger's equation number. 
Berger chooses A r so that 
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Letting n ■+■ « in Eq (1-10-5) yields an ideal Gaussian source with constant 
spectral density 
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Substitution of this 4>(u) into Eqs (1-10-1) and (1-10-2) and elimination of 6 
yields the result 



( 1 - 10 - 8 ) 


where R(D) is the information rate as a function of the MSE distortion, D. 
Hence, the MSE rate distortion function of a stationary white Gaussian source 
of average power a 2 limited to frequency band |u) | /2 it <_ B is given by Shannon's 
classic formula “ 


(4.5.78) R(D) 
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, Berger observes Chat the same R(D) curve applies to any stationary 
Gaussian process of average power o^ whose spectral density is flat over a set 
of positive frequencies of total measure B and vanishes elsewhere. In parti- 
cular, it is important to note that the one-dimensional results also apply to 
a two-dimensional raster. Since there is no correlation from pixel to pixel 
for a white bandlimited Gaussian source, there is no correlation between raster 
lines. Therefore the raster lines can be unraveled to form a one-dimensional 
process of the type we have just been examining. 

Berger plots the results of evaluating distortion as a function of rate 
and plots the results in his Figure 4.5.4 for the spectrum defined in 
Eq (1-10-5). This plot is reproduced below as Figure 10-4. 

Note chat as the normalized mean-square distortion approaches unity the 
information rate approaches zero. The white-noise source drops the most 
rapidly since it is the least correlated signal and therefore the most damaged 
by bandwidth compression. 

Pratt (Reference 143) treats a common source model for image signals, the 
separable Markovian source whose elements possess equal variance and possess 
adjacent pixel correlation factors and p c along rows and columns. He finds 
the rate distortion function for the one dimensional case to be 


R 1 (D) - \ log 2 j ^ .. I PR.>j ' (1-10-10) 
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The corresponding function for the two-dimensional case is 


RjCD) 
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( 1 - 10 - 11 ) , 


Figure 10-5 reproduces Pratt's Figure 7.7.2. There Pratt plots distortion 
functions for particular values of Pr and p c . 



FIGURE 10-5. RATE DISTORTION FUNCTIONS FOR ONE AND 
TWO DIMENSIONAL CODING OF MARKOV IMAGE SOURCE 


10,4 APPLICATION OF RATE-DISTORTION THEORY TO CALCULATION 
OF MFASMR PERFORMANCE 

In Appendix I the MSE distortion is calculated for the oT 3 ^ 2 spectrum 
which we have selected as the spectrum of microwave radiometer scenes. Again 
both one- and two-dimensional cases are treated. 

The cases treated in Section 10.2 are for signals which are continuous 
in amplitude. If the MSE distortion D becomes zero, the information rate 
becomes infinite. However, MFASMR will use digitized signals quantitized in 
amplitude. This limits information rate even if there is no distortion. 

Figure 10-6 shows that good quantizers turn in a performance not much 
worse than the classical Shannon rate-distortion curve. 

Not only is the MFASMR image signal not continuous in amplitude, it is also 
not spatially continuous. The image is characterized by the number of raster 
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FIGURE 10*6. MAXIMUM INFORMATION RATE VERSUS 
NUMBER OF QUANTIZATION BITS 


Ilium |> rune ii t . Th.Ltt uImo limits Lu L'oriuu ticm race, in order to accurate rate 
limitations due to the finite number of image lines and limitations due to 
undersampling , the spectrum used in Appendix I was normalized. The number of 
image lines was selected as an image parameter. Then we can write 


.0) 

max 


Mu 


min 


( 1 - 10 - 12 ) 


The maximum information rate versus number of image .lines was then calculated. 
The result is shown here as Figure 10-7. 

The information rate was then normalized as a percent of the maximum 
possible rate for that number of raster lines. The distortion was normalized 
as a percent of the variance. The resulting rate distortion curves are plotted 
with the number of raster lines as a parameter. The result is shown as 
Figure 10-8. 

The general region in which MFASMR systems will operate is indicated on 
Figure 10-8 as an oval centered on an information rate of 20 percent of the 
maximum rate obtained with Nyquist sampling. It may be seen that the mean- 
square distortion will be about one-third of the total signal power. 
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10.5 COMMENT ON MFASMR RATE DISTORTION CALCULATIONS 

The reaction to the calculations in Section 10.3 is one of, "What went 
wrong?" In the original Hughes proposal in Figure 2-6 a comparison image 
distortion versus image compression ratio. This ratio is the reciprocal of 
the fractional undersampling. A compression ratio of 3:1 was estimated at that 
time and a distortion of 0.1 percent was read from the graph. This figure is 
reproduced here as Figure 10-9. We now know that the compression ratio* is 
about 5 or 6 corresponding to the information rate of 20 percent of maximum 
rate. The lower result is due to the fact that we did not -count samples from 
baselines which are translated, but not rotated. In this report. We now know 
these samples are redundant. 

Figure 10-9 shows that for the less esoteric coding schemes and about 5:1 
compression we should expect about 1 percent distortion. Now. Figure 10-9 is 
data obtained with optical images. Such images usually have 1,000 lines or so. 
Figure 10-8 shows that allowable data compression is a strong function of the 
number of image lines. In Figure 10-8 we see that 20 percent undersampling 
gives about 2.5 percent distortion when the image has a thousand lines. This 
is in rough agreement with Figure 10-9. If we could increase the number of 
lines in the MFASMR image to 1,000 or so, image distortion would be reasonable. 


I 



FIGURE~10<9. DISTORTION VERSUS COMPRESSION RATIO FOR VARIOUS 
COOING ALGORITHMS 
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Unfortunately there are not enough frequency bands available for use by 
radiometers to allow us to shrink the antennas and increase the number of 
lines while maintaining an information rate which is a fixed percentage of 
maximum rate. 

The reluctant conclusion is that other frequency-domain Imager designs 
will be more satisfactory. 

Hopefully in the process of understanding MFASMR, we have arrived at a 
point where our understanding of frequency-domain imagers is complete enough to 
produce more successful designs. 

We see that the probe antennas must be very small compared to the base- 
line length. This makes the signal record long enough to measure the spatial 
frequency spectrum with good accuracy. This in turn leads to a high resolution 
imago after Fourier transformation of the record. 

This is completely analogous to the SAR radar case. The SAR radar is an 
active frequency-plane imager. Just as in the passive case, the resolution of 
the SAR improves as the antenna h±zq decreases. 

As we have seen in earlier chapters satisfactory frequency-imaging sys- 
tems for orbital use can be designed which do not use undsrsampling. Perhaps 
then wo should start out by developing such systems. This Is consistent with 
the history of radio-telescopes. The early instruments did not undersample. 
Then we can decide whether undersampling is beneficial in orbital mappers after 
exhausting the possibilities of systems meeting the Nyquist criterion. 
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1. INTRODUCTION 


,1.1 KUOIJM’JMH OK IMACJli RESTORATION 

As presented in Part I, digital image restoration is a mandatory element 
in designing a frequency domain imaging system. This Is especially true if 
the system considered corresponds to sampling below Nyquist rate. 

In the conventional imaging system, the spatial./ temporal Images are 
formed directly by hardware. Depending on the application of tliu linages, the 
tcitoval of the modulation transfer function (MTF) of the antenna or the point 
spread function (PSF) of the sensor by means of image restoration may or may 
not be necessary. In the case of frequency domain imager, the spectral data 
must at least he Fourier transformed to produce the spatial/ temporal images. 


For the time being, assume that the imaging system is linear. Following 
the notations in Reference 1, let f(x,y) be the Imaging scene, and let 
h(x,y,u,v) the system transfer function, then the sensor output g is obtained 
by 


«<x,y) 


CO CO 



h(x,y,u,v)f (u,v) 


dudv 


or, 

T{f} g 


(II-l-l) 


The problem is to find T”^ such that T~ {g} -> f. For a fully sampled 
frequency imager, h(x,y,u,v) * exp [- j 2ir (xu+y v) ] . In this case, T corresponds 
to a two-dimensional Fourier transform, and is then the inverse transform. 

In an undcrsampled system, a nonlinear sampling function s(*) is applied 
to g to produce measurements m^, i ■ 1,2,3, ... M 

» i - s(g(x,y)) + n ± 


( 11 - 1 - 2 ) 


.jmmal PAGE !?* 


* 

when n-[ is the measurement noise/error. To reproduce f(x,y) from is no 
longer a trivial task. Reference 1 pointed out that the problem of image 
restoration is an ill-conditioned problem at best and singular at worst. 
may not exist as in the case of under sampled system (singular) nor does it 
have to be unique. Even if IT 1 exists and is unique, (as in the case of fully 
sampled frequency domain Imager), it may be so ill-conditioned that for a small 
perturbation e in the measurement field g, 

T” 1 {gfc) -)■ f + d (II-1-3) 


where d is not arbitrarily small and is not negligible. Chapter 6 of Refer- 
ence 1 gives an excellent discussion of above problems. 


1.2 RATE DISTORTION THEORY FOR DATA COMPRESSION 

In cases where perfect reconstruction of the image is not required, i.e., 
under certain distortion constraints, rate distortion theory is the mathe- 
matical basis for data compression (Reference 2), Tfc should be noted that 
data compression exploits the statistical properties of the information source 
and removes redundancy, subject to a distortion measure. In this way, useful 
information is not arbitrarily discarded. 

Consider a discrete memoryless source X. It is desirable to encode X 
such that Y, an estimate of X, can be reconstructed. See Figure 1-1. 

Let p(x) be the source probability density and d(x,y) > 0 be a single- 
letter distortion measure. If p(x,y) is the joint probability of X and Y, then 


R(D) - min I(X;Y) 


) II-1-4) 


,y) “ p(x) and ^ ^p(x,y)d(x,y) £ D} 

y x y 

where I(X;Y) is the Shannon average mutual information between the two random 
processes X and Y (Reference 2). R(D) represents the minimum information rate 
that an encoder-decoder pair can operate on X such that the average distortion 
between X and Y, an estimate of X, is no longer than D. 


(p(x,y):\p(x 


X 

CAi/^nnc Mr 


DECODER 

T " A 


tNCODtn 




FIGURE 1-1. SOURCE CODING BLOCK DIAGRAM 
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It is easier to understand the mutual information if it is expressed by 
I(X;Y) - H(X) - H(X|Y) (IT-1-5) 

where H(X) is the entropy of X, H(X) - Ip(x)log(p(x)) ; and H(x|Y) is the 
conditional entropy (l.c., the remaining uncertainty of X when Y is completely 
given) ♦ To encode-decode X noiselessly, Y is required to supply all the infor- 
mation about X; therefore H(X|Y) » 0 or given Y there is no more uncertainty 
in X. Equation (II-1-5) yields R(0) - H(X) , which is the Shannon noiseless 
coding theorem; one needs the rate equal to the source entropy to code it 
noiselessly. If Y is such that X and Y are independent, then H(X|Y) ■ H(X). 

In this case, Y supplies no Information about X and thus coding/decoding process 
suffers a maximum distortion. This can be illustrated by a simple example. 

I.ot X be a coin toss process with a fair coin and Y be an independent coin toss 
process. On the average, X and Y are going to have the same outcome (head 
or tail) half the time. If d(x,y) - 1 for x ^ y and d(x,y) » 0 for x “ y 
(Hamming distortion), then the average distortion is 0.5. Therefore no infor- 
mation about the X process has to be encoded. At the decoder end, one only 
has to perform an independent experiment which has the same probability distri- 
bution as X to achieve the maximum distortion of 0.5 (R(0.5) ■ 0.). 


For a continuous amplitude stationary process (analog source), the 
entropy rate analogous to the discrete entropy goes to infinity. The extension 
of the rate distortion function to such sources is not trivial. It can be 
shown that for a zero moan discrete-time stationary Gaussian source with mean- 
square distortion measure (d(x,y) « (x-y)^), R(l>) can be parameterized by 


K(l>) w 


f max 10, lug J dw 

-IT 

■■if 


uiinl.0 ,<l»(w) ]dw 


where 4*(w) is the spectral density of the process; and 
ir 




(w) dw < 


(II-1-6) 


(II-1-7) ' 


is the variance of the process. Figure 1-2 shows few typical examples of R(D) 
curves for such Gaussian sources with various correlation properties (correspond, * 
to different spectral densities). These curves represent the optimum encoder/ 
decoder pair that one can operate and achieve the desired distortion. Equa- 
tion (II-1-6) can be extended to the time continuous case by taking the limit 
of the integrals to infinity. Figure 1-3 gives the graphic^ interpretation 
of Eq. (II-1-6) . In ideal MSE econding of the Gaussian source, it is necessary 
to reproduce only the portion of the source output at frequencies where *(w) > 6. 
If 6 is such that only the eenterlobe of 4>(w) is above 0, then an ideal low 
pass filter is the only necessary component in the source encoder. 


II-1-3 
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NOTE: NEAREST NEIGHBOR CORRELATION (M- 1/2), 

AND MARKOV CORRELATION (r - 0.78). 

FIGURE 1-2. COMPARISON OF MSE RATE DISTORTION 
FUNCTIONS FOR UNCORRELATED DATA 


PRESERVED SPECTRAL DENSITY 
XXXX ERROR SPECTRAL DENSITY 



FIGURE 1-3. SPECTRAL DENSITIES OF SOURCE, ERROR, 
AND REPRODUCTION FOR OPTIMUM CODING 
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Equation (I' 
Gaussian random 

H(D) - ~ 

8ir' 


4ir 


where u and v are the spatial frequencies of the two-dimensional Gaussian 
random field. Equations (II-1-6) and (II-1-8) which define R(D) for Gaussian 
sources, are the upper bounds to all other rate distortion functions of non- 
Cassian sources with the same spect ral den sity 4>. Part 1 uses the above 
formula applying 4>(p) ■ p"- 5 ' , p ■ / u^+v z , as the spectral density of the 
radiometer image. It gives an indication as the achievable distortion if the 
data are undersampled. 

There are few restrictions when applying the theory. The source statistics 
have to be known although work has been done to extend the theory to universal 
source coding to apply to sources with completely or partially unknown statistics. 
More severely, the theory provides only the theoretical limitations of expected 
performance. It does not in any way provide the necessary algorithm to achieve 
this theoretical limitation. As suggested by the minimization problem of 
Eq. (11-1-4) over probability distributions, the mapping from X to Y is ' 
probabilistic and is governed by p(xjy) in achieving the minimum. Often, one 
can apply Eq. (11-1-6) using the estimate of the one dimensional spectral 
density of the random process. This estimate may be obtained by average 
period Lgraiu or maximum entropy techniques. Part £ discussed some of the pros 
and cons of these techniques. However, in order to apply Eq. (II-1-8) to two- 
dimensional images, the two-dimensional power spectral density has to be 
estimated, Roucos and Childers (Reference 4) recently extended the maximum 
ontropy spectral estimation method to the two-dimensional case but does not 
provide a guarantee of the convergence of solutions. Finally, a meaningful 
distortion criterion for subjective image quality has yet to be defined mathe- 
matically so that one can incorporate it in the minimization of Eq. (II-1-4). 

Image data compression techniques have undergone extensive research 
(References 5 through 7) . In every case , the redundancy of the images is care- 
fully studied for its particular application. It is emphasized that the objec- 
tive of data compression is to produce a mapping (or transformation) T, 

(which may not be linear) , such that a suitable T*^ can be found to reconstruct 
the original images. Image restoration, however, usually deals with only the 
second half of the above problem to remove inherent system degradations. 


E-l-6) can be extended to two dimensions in order to apply to 
fields (Reference 3) , 


00 00 

; J' j' mux Jo, log dudv 


ft) tr 

If 


(n-i-8) 


min [6 ,<t>(u,v) Jdudv 
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1.3 FOURIER SYNTHESIS TEST CASES 


In simulating a frequency domain imaging system, a spatial frequency 
plane is obtained from an image f(x,y) via 


F(u,v) 


- j^y*h(x,y,u, 


v)f(x,y)dxdy 


(II-1-9) 


where h(x,y,u,v) ■ exp[-j2n(xu+yv) ) and D is the field of view. Note that 
since f(x,y) is real, F(u,v) is conjugated symmetrically. Measurements, mi, 
have to he made only on half of the frequency plane: 


UK * S{F(u,v)} + e^ i ■* 1,2,3, ... M 


(II-1-10) 


where S is the sampling function and e^ is the measurement error. Three 
sampling strategies arc used in studying the reconstruction techniques: 


Sl.fF(u,v)l 


F(u,v) if u«v, u“-v, u«0, v«0 
0 else 


(tt-l-H) 


This is similar to multi-A system's sampling of frequency axes and diagonals 
in Cartesian coordinate. This is an extremely undersampled system. For a 
64x64 image, Sl{F(u,v)} only retains 249 samples from F(u,v) over the spatial 
frequency domain coverage. 


S2fF(u,v)} - 


F(u,v) if |u| £ (1/2) u^^., 1 v | £ (1/2) v 


max 


(TI-1-12) 


else 


This is similar to ideal lowpass filtering, only the frequency samples near 
zero spatial frequency are preserved. 


S3{f (u,v)} - 


F(u,v) if |u| £ (l/4)u mnv , t |vl £ (1/4) v m ^ 


else 


(II-1-13) 


For a 64x64 image, this corresponds to low pass filtering with about same 
number of measurements as SI (256 vs. 249)., 
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Figures 1-4 and 1-5 illustrate the 64x64 stacked cube test patterns 
along with their two-dimensional Fourier transforms. Case 1 clearly has its 
spatial frequency contents concentrated along axes and diagonals. Case 2 is 
rotated slightly with respect to Case 1, in order to disperse its spatial 
frequency contents. Both cases are essentially the checker board test patterns. 

Two basic reconstruction methods are investigator 1 * The simplest one is 
the direct transform technique which includes inverse filtering and Wiener 
ft 1 terlng. For s Ingular reconstruct f.ons , Inverse and Wiener filtering are not 
expected to provide any improvements. They are included in this report because 
of their classical role in image restoration. Examples are given to demon- 
strate that both methods result in an improvement over direct inverse trans- 
forming when Lite system Ls not undursampied. The second method, due to 
Wernecke (References 8 through 11), is iterative in nature and is called 
maximum entropy method (MEM). Brief descriptions of each method are described 
in their respective chapters. In Section 4, Part II, three radiometer images 
are used to study the reconstruction algorithms and their performances are 
compared. Thoughout the report, squared-errors are computed by taking the 
sum of the squares of the differences between original and reconstructed 
images. The sum Is then divided by the total, number of pixels to produce 
estimated mean-square-error. 
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2. DIRECT FOURIER TRANSFORM TECHNIQUES 


2.1 DIRECT FFT INVERSE 

Given a set of frequency domain measurements ny. * F(ui,Vj[), 

1 * 1,2,3, ... M, the simplest reconstruction would be to perform inverse 
Fourier transform on these measurements. This is to assume that in an under- 
sampled system, the missing data take on values of zero; 

M 

A y ' 

f (x,y) « 2 _, a i (m i explj2n(xu i + yv^]} (II-2-1) 

i-1 

whore a.£ are chosen similar to the ''lag window" in classical power spectral 
density estimation or the "weighting" on each data m^. £(x,y) so defined is 

always real since it is known a priori that F(u,v) is conjugated symmetrically. 
Any set of measurements can be augmented with their complex conjugates at 
proper negative spatial frequency locations. Obviously, f(x,y) obtained 
from Eq. (II-2-1) will have any kind of resolution only if the measurement 
coverage of the (u,v) plane is almost complete. Equation (II-2-1) is aJso 
lacking in taking into consideration that an admissible reconstruction 17(x,y) 
must be nonnegative. 

Figures 2-1 and 2-2 show the reconstructions of test cases 1 and 2 for 
each of the three sampling functions SI, S2, and S3. Sampling functions S2 
and S3 represent spatial "brlckwall" filters with impulse response of sin(x)/x 
nature. The spatial convolution of the original cubes with the sin(x)/x 
response is apparent. Note that although sampling function S3 has approxi- 
mately the same number of measurements as SI, it enables u much better direct 
inverse transform reconstruction both subjectively and in squared-error sense. 
Since both -of the test cases have their frequency contents concentrated 
around zero, S3 perserves much more useful information than SI; and in some 
sense, is a better data compression algorithm. 


2.2 INTERPOLATION OF MEASUREMENTS 

There exist two major problems in direct transform techniques. The obvious 
one is that missing data in an undersampled system are assumed to have zero 
values. Secondarily, the measurements may not fall on rectangular lattice 
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e) 

FIGURE 2-2. DIRECT INVERSE TRANSFORM. TEST 
CASE 2 


TT-2-3 
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to facilitate computational efficient Fast Fourier- transform (FFT) technique. 

For example, in rotation synthesis interferometry, measurements lie on a polar 
rather than rectunyular grid. Thompson and Bracewell (Reference 12) studied 
the method of partitioning the (u,v) plane and performing "cell summing" on 
the measurements in each cell to obtain a set of data on rectangular lattice. 
Stark (Reference 13) proved a sampling theorem in polar coordinate where image 
can be reconstructed from spatial samples. It requires nonunlformly spaced 
samples correspond to scaled zeros of Bessel functions along radial direction. 

The problem of transforming these nonunlformly spaced samples onto rectangular 
grid still exist. 

Very little effort has been placed on searching for an optimum interpolation 
strategy since the inulfci-A system is significantly undersampled. Any recon- 
struction method involving severe interpolation introduces intolerable arti- 
facts. A simple method is tested (Reference 14). 

Assuming that the frequency data along axes and diagonals are sufficiently 
sampled (beyond the Nyquist rate) , therefore the entire axes and diagonals are 
completely specified, via Nyquist sampling theorem. . Measurements can be arranged 
in (6,f) plane as Illustrated by Figure 2-3. Every point on the rectangular 
grid inside the circle with radius f n, aY can be interpolated via a simple 
bicubic spline on (8,f) (Reference 15). Computational efficient FFT algorithms 
then can be applied on interpolated data. Since frequency measurements are 
complex, the real and the imaginary parts arc Interpolated independently. Fig- 
ure 2-4 shows the resulting reconstruction. on Case 1 and sampling function SI. 
Figure 2-5 has a sampling function such that angular 8 is specified every 
2 degrees from -180 to 180 as opposed to every 45 degrees for SI. The. two 
figures cieuriy indicate the relationship between artifacts and severity of 
interpolations. It is interesting to note that in order to completely specify 
every point on the 64x64 rectangular grid without any interpolation, each 
frequency quadrant (i.e. , 0-90 degree) must contain 651 distinct nonunlformly 
spaced angles. It roughly corresponds to an angular resolution of 0.14 degrees. 



f max 



b) <0. F) TRANSFORMATION FOR BICUBIC SPLINE 


FIGURE 2-3. BICUBIC SPLINE INTERPOLATION 
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2.3 INVERSE FILTERING AND WIENER FILTERING 

I live rnu and Winner filtering techniques are two classical restoration 
methods. Their usefulness is mainly to remove linear PSF or MTF distortion. 
Wiener filter is, superior to Inverse filter vhen the data are noisy and it 
reduced down to inverse filter in Lite noiseless case. Extending these classical 
filtering techniques from one dimension to two dimensions Involves applying 
theory of Tceplitz and clrculant matrices (Reference 1) . Results quoted here 
are quite intuitive, however. 

Let F(u,v) be the complex two-dimensional Fourier transform of f(x,y) 
and let H(u,v) be the system transfer function. Let G(u,v) be the measure- 
ments field and 


G(u, v) - H(u,v)F(u,v) 


(II-2-2) 


If H(u,v) is nonsingular everywhere over the (u,v) plane, then F(u,v) can be 
estimated from G(u,v) by 


F(u,v) • G(u,v)/H(u,v) 


and 


?(x.y) - // F(u,v)exp [ j2ir (xu+yv) ]dudv 


(II-2-3) 
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If all measurement (u,v) pairs on G(u,v) are on rectangular lattices, £(x,y) 
can be computed by inverse FFT algorithm by discretizing Eq. (II-2-3) . This 
inverse filtering technique is also called least-square method. 


Apparently, if H(u,v) - SI as defined in Section 1, H(u,v) is singular 
and takes on zero values over most of the (u,v) plane. Although G(u,v) - 0 
also whenever H(u,v) - 0, G(u,v)/H(u,v) gives an indeteruined ratio. To show 
an example where inverse filtering works for nonsingular H(u,v), let SI'. be 
such that 


Sl f {F(u, v) } 


j\F(u,v) u«v, u—v, u«0, v»0 
L i 1 0 else, where c small 


or, 

U(u,v) 


’1 u*v, u»-v, u»0, v*0 
. c else 


(II-2-4) 


(II-2-5) 


Oh v ion, sly liq. (11-2-2) and (11-2-3) yLcid F(u,v) ■«* F(u,v). 'flic inyerse filter 
l/h(u,v) restores f(x,y) perfectly. 

There art* two Inherent iihortcuuiLags La thLs technique. There will be 
artifacts around the edge of the reconstructed image from the circular con- 
volution wrapped around effect in calculating Eq. (II-2-3) . Also if the data 
are noisy, the noise components along with nieaHitrements are amplified by 
J./!l(u,v) when ll(u,v) is small. Figure 2-6(b) shows the noisy reconstruction 
using inverse filter 1/H(u,v) where H(u,v) is defined by Equation (II-2-5) . 

Wiener filter is derived to minimize mean-square-error between f and f 
assuming the knowledge of both the image and the noise- statistics. It can 
be shown that (Reference 1) 


A, 

F(u,v) 


and 


H*(u.v>C(u.v) 

|H(u,v)| 2 + [P n (u,v)/P f (u,v)J 


£<x,y) 


//• 


F (u , v) exp [ j 2ir (xu+yv) ] dudv 


(II-2-6) 


where P n (u,v) and Pf(u,v) are the power spectral densities of the noise and 
the image respectively; and H* denotes complex conjugate. 
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There is no ill-conditioned behavior of singularity problem associated 
with Wiener filter. Note that if H(u,v) 0, F(u,v) •+• 0 also. This implies 
that missing data are treated as zeros as in the case of direct inverse trans- 
form. If ? n + 0, i.e., noiseless, Wiener filtering of Eq. (II-2-6) reduces to 
Eq. (II- 2-3 ) , the Inverse filtering approach. Furthermore, if Pf - 0 (no signal) 
F(u,v) 0. This implies that one can not expect to recover Information at 
(u,v) where noise is dominant. 

* > 

. Figure 2-6 (a) illustrates the improvement on noise problem using Wiener 
filtering over inverse filtering technique. P n is assumed to be white Gaussian 
noise and Pf is estimated from the two-dimensional FFT of the orlgianl test 
pattern. It is well known that Wiener filter is a smoothing filter j recon- 
structed image appears flaf: and fine details arc blurred. 
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3. MAXIMUM ENTROPY RECONSTRUCTION 


3.1 DISCUSSION OF MAXIMUM ENTROPY METHOD 


In view of inadequate reconstruction techniques when the problem is 
singular, Wcrnor.ke proposed and implemented an Iterative algorithm utilizing 
the principle of maximum entropy [References 8 through 11]. When the recon- 
struction problem is singular, one has to find a solution of image that is 
consistent with the known measurements. Of course, ' solution neither has to 
oxinl nor has in he unique. A measure to pick the "heat" solution is needed. 
In addition, the .a priori knowledge that the solution (image) being non- 
negative is a very strong constraint and should be considered in the recon- 
struction process. 

Maximum entropy principle has been applied successfully to estimate the 
power spectral density S(v) of a stationary random process from measurements 
of the correlation function. There J.u no direct mapping between 8(v) and 
the entropy rate of the random process in the Shannon sense. However, the 
random process with spectral density S(v) can be viewed as the output process 
of a linear filter with squared magnitude frequency response of S(v) and 
driven by a white input process. Bartlett showed that the entropy gain of 
the filter (difference of the input and output entropies) can be expressed 
by [References 16 and 17] 


E 



log S(v) dv 


(II-3-1) 


It is not required to know the complete statistics of the input except that 
it is a white process. The entropy .of the output process is maximized if 
Che entropy gain of the filter E is made as large as possible. Naturally it 
is necessary to have the constraint that the estimated spectrum S(v) has to 
be consistent to the measured data. Maximum entropy method can be viewed as 
the "most random" extension to the known data. In the case of one-dimensional 
MEM spectral estimation, the solution is achieved via polynomial spectral 
factorization. Fart I discusses few pros and cons of this method. Two- 
dimensional extension only has limited success [Reference 4] . 
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Wernecke observed chat there are great similarities between reconstructing 
radio brightness distribution in Fourier synthesis and spectral estimations. 
Both reconstructions are subject to nonnegative constraint and are related to 
available measurement data via Fourier- transform relationship. Defining a 
similar entropy measure H, 


'■II. 


log. f (x, y) dx dy 


(II-3-2) 


the problem cun be formulated as 


if 


max I / log f (x, y) dx dy 
D* 


such that 




y) exp[-j2 ir(xu i + yv^)] dx dy 1 - m 1 (u i , Vj) 


and 


f(x, y) >. 0 for all (x, y) in D (II-3-3) 

where D is the field of view.- If the measurements are not noiseless nor error 
free , mi, i - 1, 2, .. H become m^, and additional constraint 

M 

^ ^2 !®i “ ®±| 2 " M ' ' (II-3-4) 

i-1 °i 
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has to be added to Eq. (II-3-3) . o^'s represent the variances of the 
measurements. This constrained optimization problem may not have an explicit 
solution and is hard to solve numerically. Furthermore, given a set of mea- 
surements, there may not be a nonnegative solution Shat optimizes Eq. (XI-3-3). 
The problem can be turned Into an unconstrained optimization of 


//■ 


max I f log £(x, y) dx dy - 


- *ii 2 

Q-a 


(II-3-5) 


where X plays the role of Lagrange multiplier and is chosen so Eq, (II-3-4) 
is satisfied to a sufficient accuracy. Equation (II-3-5) has a nonnegative 
solution f(x, y) for any X > 0. A numerical algorithm can be developed by 
partitioning f(x, y) into pixels (discretize the integrals) so that 
Eq. (II-3-5) can be solved iteratively [Reference 8]. A user's manual and 
program listing supplied by Wernecke is attached in, the Volume II, Appendix. 


3.2 EXAMPLES 

Figures 3-1 through 3-3 show the reconstructions of the two stacked cube 
tost cases with sampling functions Si, S2, and S3. Since uil three sampling 
functions have large number of zero values, the entropy given by Eq. (II-2-3) 
is modified slightly to stabilize the computations: 


II’ 



log (1 + f(x, y)) dx dy 


(II-3-6) 


Approximations of Sluinnon entropy associated with these images and reconstruc- 
tions are given along with the squared-errors. They are computed by using 
relative frequency approximation of the probability densities where 


■ f^/F f^ » ith pixel intensity 


F - y ^ f ^ Is the total intensity 
i 
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•I TEST CASE 1 ORIGINAL H, - 0.7 



b> TEST CASE 2 ORIGINAL H, - 0.7 

FIGURE 3-2. MEM RECONSTRUCTION, SAMPLING FUNCTION S2 
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and 


H 1 " "S P i log(p i ) CII-3-7) 

i 


In all Che cases, MEM method yields stable solutions that are both good 
ir. the subjective and squared-error sense. When the data are not severely 
undersampled (S2) , it converges rather quickly. Table 3-1 summarizes squared- 
errors for MEM and direct- transform reconstructions on the two test cases. It 
is interesting to note that for Case 2 and SI (Figure' 3-1) , the phase rotation 
of the original is lost in the reconstruction. 


TABLE 3-1. SUMMARY OF SQUARED-ERROR FOR 
TEST CASE 1 AND 2 



SI 

S2 

S3 

CASK 1 

MEM 

7.4 

8.4 

21.8 

FFjf 1 

25.5 

12.3 

19.2 

CASK 2 

MEM 

30.3 

10.2 

22.7 

FFT -1 

34.4 

12.9 

19.3 


3.3 ALTERNATIVE SOLUTIONS TO MAXIMUM ENTROPY METHOD 

Wernockc also obtained a set parametric expression using Calculus of 
Variations for various entropy criteria [Reference 8] . Tables 3-2 and 3-3 
summarize his results for both error-free measurements and measurements with 
errors. Reconstruction can also be done based on numerical solution of these 
nonlinear parametric equations. Iterative method gi.veu in Section 3.2 and 
Appendix L is more preferable since the image, f(x, y) , is manipulated 
directly rather than the parameters which has no bearing on the image quality.' 
Any one of these parametric equations can be incorporated into the MEM algo- 
rithm with only minor programming modifications. 
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TABLE 3-3. PARAMETRIC EXPRESSIONS EOR ERROR-TOLERANT MAXIMUM ENTROPY 

IMAGE RECONSTRUCTIONS 


ENTROPY/SMOOTHNESS CRITERION 


// 


in f(x, y) dx dy 


PARAMETERIZED RECONSTRUCTION 


2A 2 2 h i^ x * y) 

i-1 a i 


f f M A 

- J Jim, y)/F] 4n[f(x, y)/F] dx dy exp -H q /F - 2XF^ h i (x, y) 


i-1 w t 


■ if'*- 

jjm.. 


y) P.II r(x, y) dx dy ux|> -I - 2A £ ii^Cx, y) 


i-1 a i 


y) dx dy 


4X 2 "2 h i (x * y * 


i.=l a l 


PARAMETER DEFINITION 


J / h k (x * 

V 

■° ■ * Lh 


y) f(x, y) dx dy - k ■ I, 2, ...» M 


y) in f(x, y) dx dy 




y) dx dy 


M A 2 
r— « A . 


- 


i-1 °1 


* (if the value of A is not chosen a priori 
by the user) 
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4. RADIOMETER IMAGE TEST CASES 


Three radiometer images arc chosen Co test Clio cocoas Cruc Cion algorichms 
and sampling functions. Test Image 1 mainly consists of fields and roads. 

Test Image 2 has a mixture of fields and city area. Test Image 3 lias a more 
detailed city outline. The original images are shown in Figures 4-4 (a), 

4-L>(a), and 4-6 (a) . * Figures 4-1 through 4-3 ure Che spatial frequency dis- 
tributions of the images obtained by two-dimensional FFT. Frequency-domain 
data are simulated by applying sampling function SI, S2, and S3 to these 
Fourier-trn ns Formed Images . 

Figures 4-4 through 4-6 show the reconstruction techniques under investi- 
gation for sampling function SI. the simulated multi-,\ system. Figures 4-7 
through 4-9 show Clio results of sampling function S2 and Figures 4-10 through 
4-12 siiow the effects of sampling functions S3. Table 4-1 summarizes the 
results in terms of squared-error for all the sampling functions and two recon- 
structions. The low error values are somewhat deceiving since the original 
images have relatively low dynamic range. Figures 4-13 and 4-14 rearrange Che 
displays of the results of Test Image 1 to contrast the effects between vari- 
ous sampling strategies. 


As iu the case of s tueked-eube, S3 is clearly a better sampling function 
than Si for approximately same number of data. As a matter of fact, under the 
assumption that the power spectral density of radiometer images is unimodal 
(♦(p) ■ p“5/ 2 ) ^ rate distortion theory discussed in Section 1.2 suggests that 
unit should try Lo preserve as much low frequency region as possible. 

Maximum entropy method did not provide any improvements in all the cases 
tested except for stacked-cube case with SI. Its usefulness should not be dis- 
counted, however. The method is iterative and in all cases, it converged 
rather quickly. Table 4-2 shows the convergent reconstruction of Image 1 
with sampling function SI. In Table 4-3, MEM is given all the spatial fre- 
quency samples. It converges to a solution with insignificant error. Since 
the MEM has built in measurement error/noise and nonnegative reconstruction 
constraints, it may be a more viable restoration technique than direct inverse 
transforms even the system is not severely undersampled . 
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FIGURE 4-1 . FOURIER TRANSFORM OF IMAGE 1 



FIGURE 4-2. FOURIER TRANSFORM OF IMAGE 2 



FIGURE 4-3. FOURIER TRANSFORM OF IMAGE 3 
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FIGURE 4-4. SAMPLING FUNCTION Si APPLIED ON IMAGE 1 
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C) MAXIMUM ENTROPY dl BICUBIC SPLINE 

FIGURE 4*6. SAMPLING FUNCTION Si APPLIED ON IMAGE 2 
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FIGURE 4-6. SAMPLING FUNCTION SI APPLIED ON IMAGE 3 
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FIGURE 4-7. SAMPLING FUNCTION S2 APPLIED ON IMAGE 1 
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FIGURE 4-8. SAMPLING FUNCTION S2 APPLIED ON IMAGE 2 
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c) MAXIMUM ENTROPY 

FIGURE 4-9. SAMPLING FUNCTION S2 APPLIED ON IMAGE 3 


II-4-8 


•10344 lit 




■nnon - *.3 

M ^ • 4.3 

•) MAXIMUM KNTROPY 

FIGURE 4-10. SAMPLING FUNCTION S3 APPLIED ON IMAGE 1 


ORIGINAL PAGE IS 
Of POOR QUALITY 


•I ORIGINAL 


b) INVERSE EFT 


ERROR • S.1 
M • 4.9 


II-4-9 



ORIGINAL PAGE IS 
OF POOR QUALITY 


•) ORIGINAL bl INVERSE RPT 

■ RAOR-«.e 
• 4.3 



c) MAXIMUM ENTROPY 

FIGURE 4-11. SAMPLING FUNCTION S3 APPLIED ON IMAGE 2 
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oJ MAXIMUM ENTRORY 

FIGURE 4-12. SAMPLING FUNCTION S3 APPLIED ON IMAGE 3 
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FIGURE 4-14. MAXIMUM ENTROPY WITH VARIOUS SAMPLING FUNCTIONS 
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TABLE 4-1, SQUARED-ERROR SUMMARY WITH TEST IMAGES 


IMAGE 1 


'MAGE 2 


IMAGE 3 



TABLE 4-2. CONVERGENCE TEST FOR MEM RECONSTRUCTION 
TEST IMAGE 3 AND TEST CASE 1, SAMPLING FUNCTION SI 



Steps where 
meuuiugJL'ui 

images were 
produced 


Test Image 3 

Squared-Error 

Entropy 

12.13 

3.09 

11.30 

3.55 

11.10 

3.65 

11.04 

3.72 

10.91 

3.83 

10.89 

3.89 

. 

10.88 

3.92 

10.87 

4.00 


Test Case 1 


Squared -Err or Entropy 


36.9 

22.5 

15.1 

12.4 

10.7 

9.9 

9.7 

8.5 

8.3 

8.0 
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TABLE 4-3. CONVERGENCE TEST FOR MEM RECONSTRUCTION, 
TEST IMAGE 3, FULL SAMPLE 



Squared-Error 

Entropy 


4.14 

4.27 

Steps whore 

0.09 

4.22 

meaningful 
images were 

0.02 

4.23 

produced 

0.02 

4.23 


0.02 

4.23 
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5. CONCLUSION 


This report is by no means a complete summary of image restoration 
techniques. Several topics in the initial outline of work have not been 
incorporated in the report due to lack of manpower and time. The high recp** 
lutlon radiometer images from Chinn Lake arrived too late to be of use in 
tills computer study. Analysis of antenna HIT and sensor PSF was not carried 
out, although it is noted in Chapter 2 and Reference 1 that such linear 
effects, if nonsingular, can be removed effectively by inverse filtering 
technique for noise-free data, and Wiener filtering Cor noisy data. 

In general, undersampling is quite undesirable. Many algorithms have 
been developed to compress image data; soma of them achieve very high com- 
pression ratio for the particular class of iinuges that they are designed for. 
In ail cases, the statistical properties of the images are studied* carefully 
and algorithms designed to remove redundancy as the rate distortion theory 
suggests. One cannot expect' to recover the original source resolution when 
useful information is arbitrarily discarded. This is clearly illustrated by 
the two sampling strategies where one is simulated multi-X system and the 
other is lowpass filtering system. Both systems have approximately the same 
number of data samples. Given thnt all three test images have their spatial 
frequency contents concentrated around zero frequency, the latter preserve 
much more useful information than the multi-X system and enables a better 
reconstruction. 

Maximum entropy reconstruction method is studied to combat the under- 
sampling problem. It worked very well on the stacked cube test patterns, 
ov«i r. with simulated multi-X sampling of frequency axes and diagonals. It 
only offers a slight or no improvement over direct transform technique on 
the three test radiometer images. Since the algorithm has builtin measure- 
ment error and noise constraint, it may be a better reconstruction technique, 
»VKn the system is not undersampled. Further study Is needed to verify the 
universal worthiness of the method. 
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